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Préface 


Les lois de la physique sont des êtres bien étranges. Nous le* 
« Cotoyons » à chaque instant et pourtant elles nous semblent parfois 
aussi inaccessibles que les galaxies des confins de l'univers. 


Partez à la découverte de ces étoiles aux commandes de votre 
micro-ordinateur. Le temps de charger le programme et vous êtes à 
pied d'œuvre. Les concepteurs de votre vaisseau et de son logiciel 
ont tout fait pour que vous franchissiez sans problèmes la barrière épi- 
neuse de mathématiques qui a anéanti les tentatives de tant de vos 
prédécesseurs. 


Une fois cet obstacle passé, l'immense domaine des lois s'étend 
devant vous. Approchez de l'une d'elles : l'écran de votre moniteur 
vous en donne des images. C'est à travers de ces images que vous 
allez avoir la révélation de la loi : de même que vous imaginez la caros- 
serie d'une voiture à partir d'une série de photos prises sous des angles 
différents, de même vous allez vous familiariser avec la loi à travers 
les applications numériques qu'elle gouverne. 


Nous avons tout fait pour que votre mission réussisse mais nous 
vous laissons seul à la fin du voyage. Nous espérons que vous revien- 
drez ébloui de cette rencontre avec la loi. Désormais, vous saurez met- 
tre un visage derrière son nom, vous saurez la reconnaître sous son 
masque de frottements et d'imprécision des appareils de mesure lors- 
que vous la croiserez dans la vie de tous les jours. 


Après quelques unes de ces aventures, n'hésitez pas à venir nous 
rejoindre et concevoir des programmes de voyages vers d'autres lois. 


VII 


Sommaire 


MECANIQUE 
— Centre de gravité ............................. 1 
—: Loi: de: NeMWION: 4:55 era res es g 
— Mouvement d’un corps dans le champ de pesanteur 19 
— Effet Magnus................................. 21 
— Interactions gravitationnelles — Application à l'étude du 

mouvement des satellites terrestres .............. 28 

— Conservation de l'énergie ...................... 35 
— Plan: incliné .::::4.4.444424 urine eue 41 
+. DYNAMIQUE: : 5,5 22 Re ii es Ant 48 
— Composition de mouvements ................... 53 
— Le pendule simple ....................,........ 60 
— Le pendule de Foucault........................ 65 
— Oscillateurs couplés ........................... 72 

ELECTRICITE 
— Champ et potentiel créés par un système de charges 

ponctuelles .................................. 77 

— Circuit électrique ............................. 86 
— Circuits électriques : droites de charge ........... 94 
—+ | AéSOnNANGE:.: Linea man ous frinidhe ee ns 102 
— Filtrage RC simple ou double alternance .......... 107 

OPTIQUE 
— Lentilles minces .............................. 115 
ne LL APE RE ER EE ET RE RE ET 125 
— Dioptre plan ................................. 133 


CHIMIE 


— Propriétés d’une solution chimique ............... 137 
— Dosages volumétriques ........................ 143 
— Dissociation d’une solution chimique ............. 151 
— Couches électroniques ......................... 153 


PROGRAMMES GÉNÉRAUX 


— Equation de Van der Maals..................... 161 
— Méthode des moindres carrés — Application aux équa- 
tions non linéaires ............................ 165 


1 
Mécanique 


1. Centre de gravité 


On charge une caisse sur le plateau d’un camion. Suivant la 
charge et sa position, la répartition du poids sur les essieux est diffé- 
rente. Le camion risque même de basculer si le centre de gravité de 
l’ensemble se situe en arrière de l’essieu arrière. 


Ce programme vous montre de quelle façon les réactions du sol 
sur les essieux équilibrent le poids du camion et il vous familiarise avec 
la notion de centre de gravité. 


Pourquoi ce programme ? 


La notion de centre de gravité est importante car elle permet de 
simplifier considérablement la résolution des problèmes. Chaque tôle, 
chaque écrou du camion est soumis à la force d’attraction de la terre : 
son poids. Pour étudier le comportement du camion, il faudrait tenir 
compte de ces millions de petits poids ! Non, car lorsque vous étu- 
diez le mouvement global du camion (par exemple : va-t-il basculer ?) 
vous pouvez remplacer tous ces petits poids par une seule force, le 
poids du camion, appliquée au centre de gravité du camion. 


Utilisation du programme 
L'ordinateur vous indique la position du centre de gravité du 


camion vide (+), puis il vous demande des renseignements concer- 
nant la charge. Il vous indique alors le centre de gravité de la charge 


(+) et celui de l’ensemble camion + charge (X). Enfin il vous indique 
l'intensité des réactions du sol sur les essieux. 


Pour faire avancer le programme, taper — 


Les longueurs sont exprimées en mètres, les masses en kilogram- 
mes, les forces en Newton. Pour introduire une donnée, taper le nom- 
bre correspondant (sans indiquer l'unité), puis faites ENTER. 


Pour sortir du programme, faites ESC ESC. 


Le cœur du programme 


La caisse 


On suppose que le volume de la caisse est proportionnel à sa 
masse et que sa largeur est imposée. La surface du rectangle repré- 
sentant la caisse sur l'écran est donc proportionnelle à sa masse M” 

M” =h.c/1,6 (ligne 230) 


On suppose d’autre part que le centre de gravité de la caisse est 
au milieu de celle-ci. 


Le centre de gravité (fig. 1) 


Soit G’ le centre de gravité du camion vide de la masse M’ et 
soit G” celui de la caisse de masse M”. La position du centre de gra- 
vité G de l’ensemble est donnée par : 

(M' + M”).O0G = M’.0G” + M”’.0G” 

soit en projetant sur deux axes Ox et Oy: 

(M + M”). = M'y + My” 

(M' + M”).z = M'.z" + M'z” 

On retrouve ceci en lignes 270 et 280 avec 

M' = 5000 V' =" 292 Z' = 196 

y” = 180 (cabine) + d (espace) + c/2 (demi longueur de la caisse) 

z" = 146 (plateau) + h/2 (demi hauteur de la caisse). 





Fig. 1. — Position des centres de gravité 


Réactions du sol (fig. 2) 


Ai est la réaction du sol sur l’ensemble des deux roues, avant, 
R2 sur l’ensemble des deux roues arrière. 
Ces forces équilibrent le poids P du camion chargé 

donc : R1 + R2+ P = O 





Fig. 2. — Forces exercées sur le camion 


projetons sur un axe vertical vers le haut : 
R1 + R2-P=0 
avec P = (M” + 5000). 9,81 (ligne 320) 


donc : la somme des moments des forces par rapport à un point 
quelconque est nulle. On a choisi ce point sous les roues avant : 


R2 .(432-88) - P.(y-88) = 0 (ligne 310) 
Ces deux équations permettent de calculer R1 et R2. 


Une valeur négative de R1 signifierait que le sol tire l’essieu avant 
vers le bas. Ceci étant impossible en l'absence de boue gJluante, 
l’essieu avant quitte le sol : le camion bascule. 


Liste des variables 


C longueur de la caisse 

d distance caisse-cabine 

m masse de la caisse 

h hauteur de la caisse 

ÿ,Z coordonnées du centre de granité de l’ensemble 
r1,r2 réactions du sol sur les essieux 


Déroulement du programme 


10 présentation 
90,500 dessin du camion 
100 introduction des données 
220,500 dessin du camion 
240 dessin de la caisse 
270 calcul 
330,700 si le camion bascule 
410 fin 
illustrations 


L'illustration 3a a été obtenue avec les valeurs 15, 2, 20 000 ; l’illus- 
tration 3b avec 5, 15, 12, 500. 

Dans ce dernier cas, le centre de gravité de l’ensemble est alors 
en arrière de l’essieu arrière : le camion bascule. 






> 


+ + 
63879 181370 
: + 
le camion bascule Fig. 3a 
x 
Fig. 3b 
+ 


10 CLS 

20 LOCATE 12,8:FRINT'"Centre de gravite“ 

30 LOCATE 15,10:FRINT'"d'un camion" 

40 LOCATE 1,22:PRINT'"'pour faire avancer le program 
me,appuyer sur "CHRK#(243) 

30 REM 

69 REM 

70 q=INKEY(1):1F qt+0 THEN 79 

80 CLS 

90 GOSUB 500 

95 q=INKEY(1):1F qf>0 THEN 95 

100 FKRINT''masse du camion: 5000" 

119 FRINT: FRINT'longueur de la caisse 7 (entre 1! 
et 20)" 

120 INPUT c 

150 PFRINT:FRINT'"'intervalle entre la caisse et la c 
abine 7" 

140 PRINT'"(entre O0 et"20-c")" 

150 INFUT d 

160 FRINT:FRINT'masse de la caisse 7" 

170 PFRINT' (entre Ô et"2500*%c")" 

180 INFUT m 

190 c=20*c:d=20%*d 


200 REM 

219 CLS 

220 GOSUB 500 

230 h=1.6*xm/c 

240 PFLOT 180+d,146 

250 DRAUR €,0:DRAWR O,h: DRAWR -c,0:DRAWR O,-h 


255 TAG: MOVE 176+d+c/2,150+h/2:PRINT CHR#(228)::TA 
GOFF 
260 REM 
270 y=(232#5000+(180+d+(c/2) )#m) / (m+5000) 
280 z=(136#5000+(146+(h/2))#m) / (m+5000) 
290 TAG:MOVE y-4,2+42:PRINT"x "5: TAGOFF 
300 REM 
310 r2=(m+5000) #9. B1%(y-88) /344 

320 ri=((m+5000)*9.81)-r2 
330 IF r140 THEN GOTO 700 

340 LOCATE 6,21:PRINT CHR# (240) 
350 LOCATE 3,23:PRINT INT(r1) 

360 LOCATE 28,21:FRINT CHR# (240) 
370 LOCATE 25,23:FRINT INT (r2) 

390 REM 
400 q=INKEY(1):1F q#>0 THEN 400 

410 CLS:GOTO 100 
420 REM 
500 REM"camion" 
510 MOVE 180,144:DRAW 580,144:DRAWR O,-2 
520 DRAW 40,142:DRAWR O,-24:DRAWR 20,0 
530 MOVE 116,120:DRAW 404,120 
540 MOVE 70,144:DRAWR O,42:DRAWR 30,0 
550 DRAWR 20,40: DRAWR 60,0:DRAW 180,144 
560 MOVE 100,186: DRAWR 60,0:DRAWR 0,35 
570 DRAWR -40,0:DRAWR O,-35 

580 FOR a=0 TO 2*PI STEF O.1 
590 FLOT 88+22*C0S (a) ,118+22#SIN(a) 

600 FLOT 422+22*C0S (a) ,118+22*SIN(a) 

610 FLOT 88+26*C0S(a/2) ,118+26#SIN(a/2) 
620 NEXT a 

30 TAG: MOVE 228,140:PRINT CHR# (238) : : TAGOFF 


640 LOCATE 359,25:FRINT CHR#(247) 

650 RETURN 

700 REM 

710 LOCATE 1,1:PRINT'"'le camion bascule" 
720 Q=INKEY(1):IF qir0 THEN 720 

730 GOTO 100 


2. Loi de Newton 


Cette loi dit que deux corps s’attirent proportionnellement à leur 
masse et en raison inverse du carré de leur distance. C'est cette loi 
qui régit le mouvement des corps célestes. 


Mathématiquement, elle s'exprime de la manière suivante : 


Le corps de masse M, subit, de la part de M,, une force F diri- 
gée vers M, et égale à 


Fe kMMeS . 


F M, M, 

Inversement pour M;. 

Par ailleurs, le principe fondamental de la Dynamique dit que le 
produit de la masse M d’un corps par son accélération a est égal à 
la somme (vectorielle) des forces qui lui sont appliquées. 

Considérons alors deux corps de masses M, et M.. Ils se sou- 
mettent mutuellement à des forces attractives. (Il pourrait s'agir de deux 


charges électriques de signes opposés, dont on néglige la masse : 
les équations sont du même type). 


On peut écrire : 


_ 


Ma, = — Kk Ê M, 
r 
— F 
M, a» = — K 3 M, 
r 





Liste des variables 


X1 abscisse du mobile M, 
Y1 ordonnée du mobile M, 
X2 abscisse du mobile M, 
Y2 ordonnée du mobile M, 
V1 composante sur OX de la vitesse de M, 


w1 composante sur OY de la vitesse de M, 
v2 composante sur OX de la vitesse de M, 
W2 composante sur OY de la vitesse de M, 
Ai composante sur OX de l'accélération de M, 
B1 composante sur OY de l’accélération de M, 
A2 composante sur OX de l'accélération de M, 
B2 composante sur OY de l'accélération de M, 


Ces deux équations vectorielles conduisent aux relations 
suivantes : 


d°x: ER. (X — X2) dy; ET (V1 — 2) 
dt? (ii dt? F 
d°x2 D (X — X2) d*y2 SR (V1 — y2) 
dt? É dt? r° 


avec r = (Xi — X2) + (Yi — Ya). 
Elles s’intègrent de façon classique : 


Soit DT l'intervalle de calcul choisi ; la variation de vitesse d’un 
ensemble sur OX par exemple est égale à a, * DT. De même, la 
variation de l’abscisse x pendant le temps DT vaut v, * DT. La nour- 
velle coordonnée permet de déterminer r, donc l'accélération, et ainsi 
d'effectuer un nouveau pas de calcul. 


Nous proposons plusieurs variantes du programme de base, mon- 
trant les différentes manières dont on peut exploiter cette loi. 


— Afin d'obtenir des courbes intéressantes dans un premier 
temps, il a fallu choisir des valeurs numériques convenables ; pour 
cela, si l’on s’en tient aux phénomènes de gravitation universelle, nous 
posons : 

Ke. © G = 6,67 x 10-''MKSA 

(unité de r)°? 

— Pour les phénomènes de l’infiniment petit, m sera de l’ordre 
de 1,67 x 10? kg (masse du neutron et l’unité de distance sera de 
l'ordre de l’Angstrôm (10° m). 

Ces raisons nous ont amené à prendre K = 10 pour les applica- 
tions numériques. 


Le choix de DT est délicat : il faut prendre cette quantité assez 
petite, de façon à ce que les accroissements A1.DT, ., VA.DT, … sont 
petits devant V1, .., X1, ., mais pas trop pour ne pas rendre le pro- 


gramme trop lent. Un raffinement consisterait à asservir la valeur de 
DT à ces accroissements, c’est-à-dire par exemple à s'imposer des 
variations de X1, . constantes ou inférieures à une valeur déterminée. 


Programme n° 1 


I montre le cas de deux mobiles, lancés l’un vers l’autre, de mas- 
ses égales. On remarquera la courbure des trajectoires à proximité 
du point de rencontre. 


Programme n° 2 


ll permet d'étudier le cas où les forces entre mobiles sont toujours 
en 1/R°, mais répulsives (c'est le cas en électricité, pour deux char- 
ges de même signe). La courbure des trajectoires est moindre que 
précédemment, car nous avons pris des vitesses initiales plus gran- 
des. Le passage du programme 1 au programme 2 se fait simplement 
en changeant les signes des accélérations (lignes 50/60). 


Programme n° 3 


Il permet l'étude des chocs entre mobiles. Pour cela, les forces 
newtoniennes en 1/R? sont remplacées par des forces en 1/R’. Cela 
signifie que ces forces sont négligeables dans tous les cas, sauf quand 
les deux mobiles sont à proximité immédiate. On retrouve ici les phé- 
nomènes observés sur une table à coussin d'air ou sur un plateau de 
billard. 


Programme n° 4 


Il retourne à l'étude des forces newtoniennes en 1/R?, mais en 
envisageant cette fois-ci des mobiles de masses différentes. 


On peut écrire : 


me mm Ne sg Gm,° 

ma SG 2 Pr d'où, 4 = CU 
ii F 

a MM  , . Re Gm,* 

Mid = Gr 12°" d'où & = D 
Fe = 


Le coefficient de proportionnalité de a., soit Gm, est proportion- 
nel à m, et inversement. En introduisant ainsi un facteur 5 dans les 
expressions des composantes de a,, nous rendons compte d’une 
masse m, cinq fois plus grande que m.. Il est alors aisé de voir que 
la particule (1), plus légère donc, est nettement plus déviée que la par- 
ticule (2). 


Programme n° 5 


Il permet d'introduire l'influence d’un champ de gravitation uni- 
forme sur le phénomème d'interaction entre mobiles. Cela revient à 
incliner la table à coussin d’air. Nous avons gardé le rapport de mas- 
ses (5) du programme n° 4. 


Programme n° 6 


Il permet d'extraire de l'observation de ces phénomènes l’une des 
conséquences du principe fondamental de la Dynamique : le centre 
de gravité du système des deux mobiles se déplace d’un mouvement 
uniforme et est indifférent aux interactions entre ceux-ci. 


Programme n° 7 


Il permet de poursuivre l'étude précédente. Nous visualisons cette 
fois-ci le mouvement des deux mobiles dans le repère du “laboratoire”, 
comme nous l'avons fait jusqu'ici, mais également dans le repère du 
centre de gravité. Le mouvement (complexe) des deux mobiles dans 
le premier repère prend une forme simple dans le second : une nour- 
velle loi de la Mécanique en découle ; la quantité de mouvement glo- 
bale du système y est toujours nulle. 


10 REM EUR CE CE CE LES En sn CES LE Ci EE LE LES LES 
20 REM # attraction  *# 


30 REM + universellel * 
40 FEM # # 
50 K EM Hit RER HR ARE 
&Q REM 


70 MODE = : CLS 


GO k=10 3 dt =0.5 


90 

190 
110 
1z0 
130 
149 
150 
140 
170 
150 
170 
00 
210 
20 


250 


xXI=2O0 3 x£=<0:y1=100:y<#=300 


vis£siulis0,.4%s visttius=-0,S 


PÆSQR (CCAHE Cxi-- xs) ) TE) + (CARS (y 1-yE) ) 2) ) 


ali=-k#(x1-x#) /r°31aï=-al 
Hi=-k# (yi-yE) /r°Ssbz=-ti 
visvi+talædtiiut=ui+bhiæ#dt 
visvi+atdtims=nuz+thz#dt 
xi=xityiædtiyl=yl+tut+dt 
xé=x2+vEHadtiyi=yE+u2tdt 

IF x1 © © CR x1 % 437 THEN 
THEN 
THEN 
THEN 


IF xZ € OO OR x 
IF yE£ # © OR y£ 
IF yi < © OR yi 


ONE 


Ÿ 0j 
QU] 


on 


3 
CRC, 


x 
t 

« 
ri 


FLOT x1:y1:FLOT 
GÛTO 110 


STOP 
STOF 
STOF 
STOF 
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10 REM Hit RE EEE 

O0 REM # attraction *# 

30 KEM # universelle: # 

40 KREM # # 

MO REM tite 

£Q KREM 

7Q MODE Z s CLS 

HO k=10 2 dt =0.5 

JO x1=20 2 x2=20:yi1=100:y2=300 

100 vi=£tui=0.45s vis tiuts-O.5 

110 reSGRCCCAHES Cxl-xEe) ) 2) +6 (AH CyY1-YE) ) TE) ) 
120 al=k#(xi-xi) /r°Staï=-ai 

130 bi=kæ(yi-ye) /rtSs biz 

140 vi=vitaliædtiuli=ui+tliædt 

180 vievEt+alztdt sueur t+h td t 

160 xi=x1+virdtiyi=yli+ult+dt 

170 xZ=xE+veédti yE=yE+utdt 

180 IF x1 < ©Q OR x à THEN GOTO 240 
190 IF xZ © © OR xZ THEN GTI 240 
2OO IF yE © © OR y* THEN GOT 40 
Z10 IF yi € O OR y1 THEN GOTO 40 
220 FPLOT xi:y15FLOT 
#50 GOT 110 
#49 STOF 


D 
5j 0j 
ü 


@] 


oi mo 
a 


x 
Li 
< Ÿ 
là 
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10 REM Hot EE 
20 REM # attraction *# 
50 REM *# universelles # 
40 REM *# # 
SO REM attente 
£Q KENM 
79 MODE # +: CLS 
50 k=10 : 

90 xI=E0 3 x2= 
100 vi=itwl=0. 45: VE=Tsur=-0.S 

119 raeSGR(CCARHECxI-xE) )"2) + (CARE Cyi-y2Z)) TE) ) 
120 ai=k#(xi-xz) /rtSiaïs-al 

130 bi=ki# Cyi-ye) /r Es bE=-hi 

140 vi=vitalætdtimismit+tliædt 

150 visvetastdtsiutenus+hzt dt 

140 x1=x1+visdtiyi=yituiædt 

170 xE=x2+vVetdts yes y Edited t 

180 IF xi < O OR xi + 437 THEN GÜÛTO #40 

190 IF x2 & O UR x£ 5% 439 THEN GTO 240 

200 IF yE£ € © OR y# + 377 THEN GOTO 240 

210 IF yi < © OR y1i > 5977 THEN GOT 40 

220 PLOT x1:yisPFLOT xZ:yE 

#30 GQÛTO 110 

#40 ETOF 





EO: ÿ 1=100: ÿY£=500 


y£= 2#43.47:874 yi= 140.527135 


D 


ne UE 


Pen 


200 
£210 


20 


30 


Z40 
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REM 
REM 
REM 
REM 
REM 
REM 


idee he A EE EE 
+ attraction # 
# mniverselle4 + 
# +: 


RE 8 LE LE LE 2 28 LE Ge a LES LES si SA cd 


MÜDE 2 3 CLS 


H=10 dt =0.5 


=2O 5 x2=20:y1=100:y2=500 

Re ve=iiwgez-0,48 

FE=SCR (CCABS C(x1-x2) ) 2) + (ABS (y1-y2))*2)) 
ali=-Stk#(x1-x2) /r*3iai=-a1/s 

bi=-Stk# Cy1-y2) /r°%35s5b2=-b1/5 
visvitaiædtsul= rh 
vV£Z=vitarstdtius=uz+bhz#dt 


xi=xit+tvisdtiyi=yi+uwuirdt 
x£2=xE+vetdtiye=y£tuztrdt 
IF x1 < © OUR xi % 437 THEN GOTO Z40 
IF xZ2 < © OR x2 > 637 THEN GÛTO 240 
IF y? £ © OR y£ * 399 THEN GOTO 40 


PS 


IF yi < Q OUR yi + : 
PLOT xt:y1: FLOT x y: 
GOTO 110 

STÜF 


1 
Ru] 
ÿ 


THEN GÜTO 240 


ti 





Programme Attraction universelle 4 


(para 


10 
£O 
30 
40 
50 
40 
70 
50 
30 
100 
110 
120 
130 
140 
130 
140 
170 
150 
170 


#00 


mètre : 5 des lignes 120 et 130 remplacé par 50) 


REM Hit te 
REM # attraction *# 
REM # universelles * 
REM # # 
REM otttotntt 
REM 

MALE 
kH=10 : 


xi=#0 © 





u 
D 2) + (ABS (y1-ÿ2))°2)) 


rESER (CCARE (x i1-x2) 
al=-S+rkit (x1-x2) /r°3ari=-al/s 
bi=-Srkæ Cyi-yE) /rt3ibzi=-t1/5 
vi=vi+tairdtiui=mi+(bi1-0.02)+#dt 
ve=vitasw#dtiuz=us+(bz-0.02) #dt 


xi=xi+vitdtiylzylitiuilædt 
xE=x2+vEtdte yE£=yEtiusédt 
IF x1i < © OR xi 5% 4397 THEN GOTO 
IF x£ © O OR x2 > 637 THEN GOTO 
IF y£ < O OR y£ > 377 THEN GOTO 


bi 5 1 
+ + 
Co © 


ES 
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10 IF yi ? © OR yi x 3 


3 
220 PLOT x1:y1:FLOT x23 y: 


7 THEN GOT 40 


“530 GOTO 119 
“40 STOF 


19 
zO 
40 
40 
50 
40 
70 
&O 

70 

100 
110 
120 
130 
140 
150 
160 
170 
150 
170 
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REM tie 4e AE A 
FREM # attraction  # 
REM # universelles *# 
RENM # # 
REM + 

REM 

MÛDE # 3 LS 

E=10 3 dt =0.5 


XI=EHD 5 xZ=rO0ivyi1=100:yEZ=3500 


LS CRÉES LES CE LES DES DES LES LES DES LES SES LES LES LED 






Vviziiuli=0.45i VE=ius-OQ . 45 
r= 
al=-kH#(xi-x£) /r"35a2z-ai 
Hi=-k#(yi-y£) /r°SiEz=-h1 
visvitalætdtiui=mi+biædt 
visve+ai#dttuseus+b td t 
xi=xi+viædtiylizyitiunitdt 
xé=xe+vetadtiyéeyEstusedt 

IF x1 © © OR x1 % 437 THEN GOT «40 
IF xx © O UK x2 % 4379 THEN SOTG 40 


EUR (CCARS CxX 1 XIE) ) TE) + (CABS (y1-y£) ) 72) ) 


IF y£ © © OR y# x 377 THEN GÜTO X40 
IF yi © © OUR yi > 377 THEN GOTO 40 
FLOT x1:y1:FLOT xZ2 y 

PLOT Cxi+xe) /22 Cyi+y£) /E 

250 QT 110 

240 STOF 

250 STOP 









Dee 


LO  REM A MRRER AE AE TE AO MEME EE 
5 REM # attraction  # 
BEM + universelle] # 
REP + # 
Ne] [a] F4 FT M EE ES RS ES OS ES DS RS 2e RES CE CES © ES LE ES D 





79 MODE & 2 LE 

0 LOCATE 40: 4%FRINT centre de gravite® 
90 LOCATE EG: 20 | i 

100 Elo 4 dt 
110 x1=106 5 Ê 
120 vlie-0 ,4rmirlrs tr. 
150 PeSOR € CARS Crdex st) ) UE) + CARS Cyl-y) ) TE) 





“laboratoire! 
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ai=-k 
bi=-E 


ÆCxI-xE) /rTSrars=-al 
ÆCyi-ys) /r°%3s b£s=-hi 


vizvital#tdtiult=uittlædt 


VE=VE 
xi=xi 


+aztdtiui=we+t dt 
+viædtsy1=y1i+tiuiæadt 
xE=X EAN EH dt y y + dt 


IF xi Q OUR xi 437 THEN GOT 20 

IF x2 © © OR x 457 THEN GOÛT ZEO 

IF y£ © O OR y* 373 THEN GOT 280 

IF yi © © CR yi 5% 377 FHEN GOT £50 
FLOT  ZO0+x1-Cxi1+xE) /2: 3041 Cy1i+y) /E 
FLOT 2004 x2— (x 1442) /2 * 320+y2- (Yi +yx) /X 
PLOT xl: y1sFPLOT x y 

STE 130 

STOF 


STOF 


A 
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\ 
à. centre de gravite 
\ 


laboratoire 


3. Mouvement d’un corps 
dans le champ de pesanteur 


Ce problème est classique en Terminale. L'exemple cité est sou- 
vent celui de la bombe lâchée par un avion. La décomposition de la 
force (unique : c’est le poids de l’objet) sur deux axes orthogonaux 
OX et OY conduit, après intégration des équations différentielles linéai- 
res, d’une part à un mouvement rectiligne à vitesse constante sur OX, 
d'autre part à un mouvement rectiligne uniformément varié sur OY. 


Le programme propose de montrer, dans l’ordre, le mouvement 
de l’objet suivant OX, puis celui suivant OY, enfin la composition de 
ces deux mouvements qui donne l'aspect de la trajectoire pour un 
observateur terrestre. Dans le référentiel lié à l’avion, c'est la seconde 
observation qu'il faut prendre en compte. Dans un dernier temps, nous 
tentons de montrer l'influence de la résistance de l’air par l’intermé- 
diaire d’une force de frottement proportionnelle au carré de la vitesse 
de l’objet qui tombe. 

Pour rendre plus générale l’utilisation de ce programme, l’utilisa- 
teur doit introduire la vitesse initiale V (en mètres par seconde), l’angle 
de départ T (en degrés), les coordonnées de départ (X, Y) et éven- 
tuellement le coefficient de la force de frottement K (qui peut être com- 
pris entre — 0,001, force très faible, et —0,1, force très élevée), 


On retrouve facilement les cas classiques : 
— obus lancé d’un avion 


V = quelconque  X = 0 
T = 0 à -90° ; Y = 300 


lignes 80 à 110 
— obus lancé par un canon 


= quelconque : X = 0 
T = compris entre O0 et +909; Y 150 


< 
I 


10 REM HAE RAR EE MEET 


20 REM # + 
30 REM # chute litre *# 
40 REM *# # 
SO REM tt tt 
&40 REM 


70 MODE 2 : g=7.81 


50 
70 
100 
110 
120 
130 
140 
150 
140 
170 
180 
190 
Z00 
210 
220 
30 
£40 
250 
60 
270 


INFPUIT''Vitesse initiale: "s:vivo=v 
INPLIT'"'angle de depart:"st:TO=t 


INFLIT'abcisse de depart:s"sx:xo=x 
INFUT'"'ordonnee de departi";ysyo=y 
PLOT x:y:DRAWR 437-x:0 

PLOT x: y: DRAWR O0: 399-y 

FOR n=1i TO S:READ k 

V=VOoi x=x0oy=yos t= TO 

t=t*PF1/180: dt=0,.01 
VX=VHEOUS (t) 3 Vy=vV#SINCt) 

V= (vx "2+vy"2) "0.5 

Fx=kEkVEVXxS FY=KERVEVY 
ax=-fxiay=-fy-g 
vrk=vxtaxétdt:svy=vy+ay#dt 
x=xtvrxtdtiy=yt+vy#dt 

IF x>640 UÜR y3400 ÜR yÿy<o THEN 270 
PLOT x>y 

DATA 0:0.001:0.005:0.01:0Q.1 

GÜTO 180 

NEXT n 


PA Pi AE igitiale; ? 65 
e LA e à A 
s 

bise art 9 158 


rdonnee 
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itesse initiale:? 89 
ngle de depart:? 18 
bcisse de depart:? 8 
rdonnee de depart:? 289 


NT 


itesse initiale:? 98 

ngle de depart:? 0 
cisse de depart:? 8 

rdonnee de depart:? 308 








4. Mouvement d’un corps lancé 
Effet Magnus 


Les programmes suivants permettent d'étudier les trajectoires de 
corps lancés avec une vitesse initiale oblique dans le vide ou dans 
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l'air. Dans le cas particulier où le corps en mouvement est sphérique 
(cas d’une balle de tennis, de golf ou d’un ballon de football), on peut 
faire intervenir la force Magnus. L'effet Magnus étant dû à la différence 
de pression qui apparaît entre deux points diamétralement opposés 
d’une sphère lorsque celle-ci tourne. 


W > O balle liftée 





W < O balle coupée 





balle brossée 


Rappelons que la résistance de l’air est de la forme : 


— 


V 


Fa = — KSv2 À 
IV 


Pour une sphère S = x R? (section droite). 

Dans le cas d’une balle, la force Magnus est de la forme 

FM = AWA V 

K et À sont des coefficients qui dépendent de la balle et aussi du 


domaine de la vitesse. Ce sont des grandeurs déterminées 
expérimentalement. 


Dans le cas des forces Fa et FM, c'est la vitesse relative de 
l’air et de la balle qui intervient. Dans le cas où la vitesse du vent n’est 
pas négligeable devant celle de la balle, il faut calculer cette vitesse 
relative. 

Si le vent s'oppose au mouvement de la balle, il faut donc aug- 
menter Fa et FM à vitesse égale de la balle par rapport à l'observateur 
(VR = Vballe —- V vent). 


Appliquons la relation fondamentale de la Dynamique. 


My = Mg -KSV? _ AW 4 V 
VI 


Cette équation différentielle sera calculée numériquement dans 
les programmes. 
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Pour une balle de golf, par exemple 


KS 
M 


= 8.10 ° 


et 27À 2 210-° 


M 


Liste des variables 


V, VO et VF 


, WO, WF 


OZ >rx 


_ 
£ 


Déroulement du p 


Dans le program 


vitesses de la balle 

(on fait, dans une boucle, varier la vitesse V de 
VO à VF) 

coefficient égal à (KS/M) (Résistance de l’air) 
coefficient égal à 2x x A/M (effet Magnus) 
valeur de l’angle que fait la vitesse VO avec 
l'horizontale 

vitesses de rotation de la balle en tours/seconde 
accélération de la pesanteur 

Temps de l'observation (durée maximale 
prévue). 


rogramme 1 


me 1, on étudie les trajectoires d’une balle de 


golf liftée, l’angle de la vitesse initiale avec l'horizontale étant égal à 
5 degrés. On calcule 600 points pour 5 valeurs de W différentes (ligne 


145). 


lignes 10 à 125 
ligne 140 


ligne 200 


lignes 200 à 400 


lignes 410 à 560 


titre et entrées des données 

introduction des constantes 

pas du calcul (pris égal à 1/100 s) 

on fait varier la vitesse de rotation de la 
balle W de 0 à — 100 tours/seconde (balle 
liftée) avec une vitesse initiale constante 
de 80 m/s faisant un angle A avec l’hori- 
zontale égale à 5 degrés 

résolution numérique de l'équation diffé- 
rentielle du mouvement (projection sur 
l’axe OX et projection sur l’axe OY), les 
paramètres S et R permettent de tracer 
les courbes à l'échelle prévue par l’utilisa- 
teur (ligne 470) 

on trace les courbes sur l'écran et on 
imprime les différentes données. 
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Déroulement du programme 2 


Dans le programme 2, on étudie les trajectoires fictives d’une balle 
de golf liftée (drive) à laquelle on aurait communiqué une vitesse ini- 
tiale de — 300 tours/seconde, la vitesse initiale de celle-ci variant de 
10 à 80 m/s. 


L'angle de la vitesse avec l'horizontale est pris égal à 45°. 


19 REM EFFET MAGNUS 

20 REM 362666 

30 REM 

40 REM ENTREE DES DONNEES 

50 REM HERO EEE 

60 CLS 

70 INPUT "VALEUR DE V="34V 

80 INPUT "VALEUR DE ALFHA A=";A 

90 INFUT "VALEUR DE K="3;K 

100 INPUT "VALEUR DE L='";L 

110 INPUT "VALEUR DE WO ="; WO 

120 INPUT "VALEUR DE WF="; WF 

125 INPUT "TEMFS DE L'OBS.TM="3; TM 
130 CLS 

135 S=0:R=0 

140 G6=9.81:Z=A*P]1/180:DT=0.01 

142 LOCATE 12,12:PRINT "CALCUL EN COURS" 
145 DIM A(5,600):DIM B(5,600) 

150 REM ETUDE DES TRAJECTOIRES 

160 REM #26 DEEE DEEE DE DE MED DE DE HD 

165 N=0 

200 FOR W=WO TO WF STEP -25 

2905 NEN+1is:s LOCATE 3,15: FRINT N 

206 IF N>5 THEN GOTO 410 

210 X=0: Y=O: TX=V#COS(Z): TY=V#SIN(Z) 
215 M=0 

220 FOR T=0 TO TM STEF DT 

225 M=M+i 

226 IF M=600 THEN GOTO 400 

230 ATX=-K#SOR(TXATX+HTYATY)#TX+LEWETY 
240 ACTX=ATX#DT: TX=TX+ACTX:ACX=TX#DT 
250 X=X+ACX:A(N,M)=X 

259 IF ABS(X)>S THEN S=ABS(X) 

260 ATY=-KASORITXEATX+HTYEATY) #TY-LAWXTX-G 
270 ACTY=ATY#DT: TY=TY+ACTY: ACY=TY#DT 
280 Y=Y+ACY:BN,M)=Y 

285 IF ABS(Y)>R THEN R=ABS (Y) 

290 IF Y£0 THEN GOTO 400 

320 NEXT T 

400 NEXT W 
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410 
420 
430 
440 
445 
450 
460 
470 
480 
490 
200 
910 
920 
930 
940 
930 
960 


nie 


CLS 
MODE 2 


FLOT 0,0:DRAWR 640,0 
PLOT 0,0:DRAWR 0,400 


LOCATE 


3,2:FRINT 


FOR I=i TO N 
FOR J=1i TO M 
FLOT 4(1,J)/Sx600,B(1,J) /R#x300 


NEXT J 
NEXT I 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 


55,3:PRINT 
55,4:PRINT 
55,5:PRINT 
55,6: FRINT 
55,7:PRINT 
55,8:PRINT 
55,9:PRINT 


"Y'sLOCATE 78,24:FRINT "X" 


“V=U;V; "M/S" 
“A="; A3 "DEG." 
"KaUsK 

QC DETE 
“WO=" ; NO; "T/S" 
“WF="; WF 3 "T/S" 
“TM=" 3; TM 
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10 

20 

39 

40 

50 

60 

70 

75 

89 

90 

100 
110 
125 
130 
135 
140 
142 
145 
150 
160 
165 
290 
205 
206 
210 
215 
220 
225 
226 
230 
240 
250 
255 
280 
279 
280 
285 
290 
320 
490 
410 
420 
450 
440 
445 
450 
469 
470 
480 
490 
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REM EFFET MAGNUS 
REM 3636 2626626 6 EX 6 
REM 
REM ENTREE DES DONNEES 
REM 62626066 DE DE EH DE DEEE DEN DE 
CLS 

INPUT "VALEUR DE VO="3: VO 

INPUT "VALEUR DE VF="3sVF 

INPUT "VALEUR DE ALPHA A='"';A 
INPUT "VALEUR DE K='"3;K 

INPUT "VALEUR DE L=";L 

INPUT "VALEUR DE W='";W 

INPUT “TEMPS DE L'OBS.TM=";TM 
CLS 

S=0: R=O 
G=9.81:2Z=A*PI/180:DT=0.01 

LOCATE 12,12:PRINT "CALCUL EN COURS" 
DIM A(5,600) : DIM B(5,6400) 

REM ETUDE DES TRAJECTOIRES 

REM 26262626 D DD DE DD DEEE DE DE DE DEEE 

N=O 

FOR V=VO TO VF STEF 20 
N=N+1:LOCATE 3,15:FRINT N 

IF N>5 THEN GOTO 410 

X=0: Y=O: TX=V#COS(Z): TY=VYSIN(Z) 
M=0 

FOR Tz0O TO TM STEF DT 

M=M+i1 

1F M=600 THEN GOTO 400 
ATX=-K#SOR(TXATX+TYATY)ATXHLEWXTY 
ACTX=ATX#DT: TX=TX+ACTX:ACX=TX*DT 
X=X+ACX:A(N,M)=X 

1F ABS(X)>S THEN S=ABS(X) 
ATY=-KASQOR(ITXATXHTYATY) #TY—-LEUWXTX-G 
ACTY=ATY*#4DT: TY=TY+ACTY: ACY=TY*#DT 
Y=Y+ACY:B(N,M)=Y 

IF ABS(Y)2R THEN K=ABS(Y) 

1F Y<O THEN GOTO 400 

NEXT T 

NEXT V 

CLS 

MODE 2 

PLOT O0,0:DRAWR 640,0 

PLOT 0,0: DRAWR 0,400 

LOCATE 3,2:PRINT "Y'":LOCATE 78,24: FRINT "X" 
FOR I=i TO N 


FOR J=1 TO M 

PLOT A(I,J)/S#x600,B(I,J) /K#2:00 
NEXT J 

NEXT I 


900 LOCATE 55,2:PRINT "VO="; VO; "M/S" 
508 LOCATE SS,S:PRINT “VF="3; VF: "M/S" 
910 LOCATE SS,4:PRINT “A=";A;"DEG." 
929 LOCATE S5S,5:PRINT “"K='";5k 

930 LOCATE 55,6:PRINT "L=";l 

540 LOCATE 55,7:FRINT “W=";W;"T/S" 
990 LOCATE S55,9:PRINT "“TM="3 TM 
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5. Interactions gravitationnelles 
application à l’étude du mouvement 
des satellites terrestres 


La loi de gravitation exprime que deux masses M1 et M2 situées 
à une dinstance d l’une de l’autre s’attirent mutuellement avec une 
force d'intensité : 


F=cG M1M2 
d? 


où G est la constante de gravitation universelle (G = 6.67 x 
107" SI). 

Etant donnée la faible valeur de G, l'intensité F ne devient appré- 
ciable que lorsque les masses mises en jeu sont importantes. 


Ce programme permet de déterminer l'orbite parcourue par un 
satellite terrestre de masse m lancé à partir de l'altitude HO avec une 
vitesse initiale VO faisant un angle donné A, ni nul, ni plat avec la droite 
OM qui joint le centre de la Terre au centre de gravité du satellite. 





Remarque : 

La lune sera toujours considérée dans une position telle que la 
perturbation introduite par celle-ci sur l’orbite du satellite soit consi- 
dérée comme négligeable. 


Etude de l’orbite du satellite 


Expression de la force F qui agit sur le satellite : 


Fe - GP avecR =RT+H 
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Energie d'interaction de gravitation : W(R) = — G MAL 
R 
Les intégrales premières du mouvement traduisant la conserva- 


tion du moment cinétique et celle de l'énergie s’écrivent 
respectivement : 


MAR ny (J/m = constante des aires), 
dt 

_ QMTM 1m [Re à re Be | LE 
R 2 dt? dt 





En éliminant le temps entre elles, on obtient l'équation de la tra- 
jectoire du satellite en coordonnées polaires. 








dB de dR dR,, dB J: on 
R? CPR 2 2 CDS 2, 
Cat m°R? Ca (8? Ca m°R* Ge) 
d'où 
2 
_gMTm, J 1 fe ; à LE 
R 2mR° | R° dB 
Effectuons le changement de variable 
1 du 1 dR 
u = donc = — 
R dB R? dB 
il vient 
2 
-GMTmu + : pos + W| =E 
2m| dB 
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décrivons par rapport à B, 


2 2 
OL ENT EU 


dB Re lin 


La solution = 0 conduit à u indépendant de B, c'est-à-dire à 


des trajectoires circulaires ; la deuxième solution est celle de l’équa- 
tion différentielle : 
2 2 
CARE MTm 
dB? Hd 
On rappelle que l'orbite du satellite est une conique. La nature 


de cette conique dépend uniquement dela grandeur (et non de la direc- 
tion) de la vitesse initiale. 


Pour que le satellite puisse s'éloigner indéfiniment de la Terre, 
il faut que l’on ait E > O, soit, d’après les conditions initiales 


gi + À mvo? = 0 
(RT + HO) 2 
l'orbite ne pourra comporter une branche infinie que si 


Fycse..-Mr 


2 (RT + HO) 


La vitesse minimale à communiquer au satellite appelée vitesse 
de libération VL est donnée par 


re | 2GMT 
(RT + HO) 
— Si VO < VL, l'orbite est elliptique, 


— Si VO = VL, l'orbite est parabolique, 
— Si VO > VL, l'orbite est hyperbolique. 


Dans les deux derniers cas, l'orbite ne comprend que la branche 
de parabole ou d’hyperbole dont la concavité contient le foyer O, à 
la fois centre de la terre et centre d'attraction. 


Liste des variables 


G constante universelle de gravitation 
MT masse de la Terre (5.98E24 kg) 
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RT rayon terrestre (6 370 km) 

HO altitude initiale du satellite 

VO vitesse initiale ns 

A angle de VO et du rayon vecteur initial OM (en radian) 

B angle du rayon vecteur initial et du rayon vecteur OM à l’ins- 
tant t (varie de 0 à 2P1 au cours d’une révolution) 

C constante des aires (C = J/m) 
(on rappelle que C = VO (RT + HOjsinA 


U= :avecR = AT +H 


Déroulement du programme 


Le programme permet de tracer sur l'écran l'orbite du satellite 
et de faire apparaître à l'échelle de celle-ci la Terre ; dans le cas où 
les orbites sont elliptiques, on calcule la durée de révolution, l’altitude 
et la vitesse à l’apogée. 


lignes 10 à 60 
lignes 65 à 90 


ligne 180 


ligne 185 


ligne 186 


ligne 190 


ligne 210 
lignes 240 à 310 


lignes 320 à 350 


introduction 

introduction des conditions initiales nécessai- 
res à la résolution numérique de l'équation dif- 
férentielle de mouvement 

dans le cas d’une orbite parabolique ou hyper- 
bolique, on arrête le calcul lorsque la distance 
R = RT + H est supérieure à 150 fois le rayon 
terrestre 

dans le cas d’une orbite elliptique selon les 
valeurs de A et de VO, le satellite peut retour- 
ner sur Terre, on arrête donc le calcul lorsque 
R < RT 

on calcule, dan le cas des orbites elliptiques, 
la durée d’une révolution ou le temps que met 
le satellite à revenir sur Terre. 

calcul du paramètre M qui permet de tracer 
l'orbite à l’aide d’une échelle calculée 
automatiquement. 

résolution numérique de l'équation différen- 
tielle du mouvement 

tracé des axes et d’un cercle correspondant 
à la Terre à l'échelle de l'orbite du satellite. 
tracé de l'orbite du satellite. 
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10 
20 
30 
25 
40 
50 
60 
65 
79 
80 
90 
100 
110 
20 
130 
140 
150 
160 
170 
180 
185 
186 
187 
188 
190 
200 
210 
220 
225 
(W) 
2:30 
240 
245 
246 
230 
260 
270 
280 
290 
500 
310 
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325 


Le 
2 
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lignes 400 à 430 impression sur écran des différents 
paramètres. 


RKEM GRAVITATION 
REM HE EEE 
REM 
CLS 
REM INTROD.DES CONSTANTES 
REM EDEN DE DE DEN DE DEEE 
G=é6.67E-11:RT=62:70000:MT=S. 98E+24 
INFUT "PAS DU CALCUL DE=";DB 
INPUT "ALTITUDE INITIALE HO="; HO 
INFUT "VITESSE INITIALE VO="; VO 
INFUT "VALEUR DE ALFHA="3; A 
LOCATE 12,1S:PFRINT "CALCUL EN COURS" 
C= (RT+H0) #*VO#xSIN(A) 
É=GxMT/(C"2) 
L=2*FI/DB+1:DIM A(L):DIM H(L) 
U=1/ (RT+H0):TU=-(1/0C) #VOxCOS (A) 
M=O:N=0:F=0:Q=03: T=0 
FOR B=0Q TO 2*PFI1 STEF DE 
N=N+1:A(N)=1/U 
IF ACN)2>1SOKRT THEN A(N)=0:F=1 
IF A(N)E=RT THEN G=i 
IF Q=0 AND F=0 THEN H(N)=(1/U)"2/C*xDH 
IF G=1 AND F=0 THEN E(N)=0 
T=T+B (N) 
IF ABS(A(N)):M THEN M=ABS (A (N) ) : WEN 
LOCATE 1,17:FRINT N 
DU=E-U: AD=DU%#DE: TU=TU+AD: U=U+TU*DE 
NEXT E 
IF Q=0 OK Q=1 AND F=0 THEN K1=((A(W+1)-A(W))/E 


):KH2=C/A(W) : V=ESOR(K1KE1+KZ#+2) 
CLS 
MODE 2 


LOCATE 50,2S:PRINT INTC(M/19.9%0.001)+0.5S)3"EM 


FLOT 400,32: DRAWR 10,0 
FLOT 2320,0:DRKAWR 0,400 
FLOT 0,200:DKRKAWR 640,0 
FOR Z=0 TO 2*FI STEF 2#F1/:3:00 
X=RT/M#199xC0S(Z) +320 
Y=RT/M#199xSIN(Z) +200 
FLOT X,Y 

NEXT 7Z 

N=O 

FOR B=0 TO 2*FI1 STEF DE 
N=N+1 


327 IF A(N)=0O OR A(N)#O THEN GOTO 400 

320 IF A(N)£E=RT THEN GOTO 400 

340 FLOT 220+A(N) /Mx199xC0S (EH) ,200+A(N) /M#199xSIN( 

E) 

350 NEXT E 

400 LOCATE 1,19:FRINT "HO="3:HO%O0. 001: "KM": LOCATE 1 

s2OS PRINT "VO="SVOXxO. 001: "EM/S":LOCATE 1,21: FRINT 

"A='!"3 INT(A*180/F1+0.5) 3; "DEG. " 

405 IF A(N)=O OK A(N)£O THEN STOFP 

410 IF Q=0Q THEN LOCATE 1,1:FRINT "DUREE DE REVOLUT 

ION=";T/3600; "H" 

420 IF Q=1 AND F=0 THEN LOCATE 1,1:FRINT "DUREE KE 

TOUR="3;T/2:600; "H" 

4350 IF Q=0 OK Q=1 AND F=0 THEN LOCATE 1,2:FRINT "A 

LTITUDE AFOGEE=" 3 INT ((M-RT) #0. 001+0,.5):; "EM": LOCATE 
1,3:FRINT "VITESSE AFOGEE="3;Vx0,001,;"KHM/S" 

440 STOF 





= 6508 
de 96 DEC. 
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DE REVOLUTION= 1.49982783 HT" Rene 
ALTITUDE APOCEE= 373 KA Re, 
VITESSE APOGEE= 7. Ars US Tin, 





Fi RETOUR= 7 TN C] 
VI ge aPOLÉES 8° 8.876333287 M/S este 





DUREL RETOUR: 16.5231487 N 


DLTSSE = 0. DT &s 








6. Conservation de l’énergie 


Le wagon du grand-huit va-t-il réussir à monter jusqu'au sommet 
suivant ? Ce programme est destiné à montrer comment on peut résou- 
dre des problèmes de physique à partir du principe de conservation 
de l'énergie. Il étudie le mouvement d’un objet que vous lancez sur 
une route dont vous donnez le profil. 


Pourquoi ce programme ? 


Le principe de conservation de l'énergie est une « règle du jeu » 
que les physiciens tiennent absolument à respecter. 


Ils définissent pour cela de nouvelles formes d'énergie chaque 
fois qu'ils découvrent un nouveau domaine. Il faut dire que ce prin- 
cipe est très efficace pour résoudre de nombreux problèmes. 


Utilisation du programme 


Les données sont introduites sans unité, suivies de ENTER. 
Pour faire avancer le programme, taper — 
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Pour sortir du programme, faire ESC ESC puis RUN ou RUN 200 
ou MODE. 


L'objet a pour masse 1 kg. Les longueurs sont en mètres, les temps 
en secondes. Vous indiquez le profil de la route au moyen d’une équa- 
tion en ligne 1 000 et la loi de frottement en ligne 1100. Vous pouvez 
éviter de recopier les lignes 10 à 190 qui ne sont qu’un préambule. 

L'objet a, au départ, l'énergie E. L'ordinateur indique à chaque 
instant la valeur des énergies cinétique et potentielle sous la forme 
d’un pourcentage de E. Quand il y a frottement, de l’énergie mécani- 
que disparaît et la somme EC +EP est différente de 100 %. 


L'ordinateur augmente l’abscisse de l’objet de 4 mètres entre cha- 
que image et indique le temps écoulé depuis le départ. Il ne s'écoule 
donc généralement pas un temps fixe entre deux images. 


Les résultats peuvent être erronés, on observe par exemple des 
oscillations, lorsque la vitesse de l’objet est très faible. 


Ne donnez pas une vitesse nulle au départ de l’objet. 


Le cœur du programme 


Y est l’altitude de l’objet. E est l'énergie initiale. L'énergie cinéti- 
que est donnée par EC = (1/2).M.V?, ici : 0,5.V?. L'énergie potentielle 
est donnée par EP = M.g.Y, ici : 9,8.Y. En l'absence de frottement, 
EC + EP = E = constante. 


Lorsque le point avance de dX en abscisse, son altitude varie de 
dY. Ici, dX n’est pas infiniment petit mais égal à 4, on veillera donc 
à ce que le profil de la route n'ait pas d'accident inférieur à une dou- 
zaine de mètres en abscisse, sinon on corrigera la ligne 270. 


L'ordinateur calcule la nouvelle énergie potentielle du point à partir 
de la nouvelle valeur de Ÿ et en déduit la nouvelle valeur de l'énergie 
cinétique par EC + EP = E. Il en déduit alors la nouvelle valeur de 
la vitesse. 

S'il y a frottement, l'énergie de l'objet a diminué du travail de la 
force de frottement, soit 


F. V (dx)? + (dY}? 
Tout se passe bien tant que la différence E — EP n'est pas néga- 


tive, car cela aboutit à une énergie cinétique négative ce qui est impos- 
sible : l’objet fait alors demi-tour (lignes 560 et 550). 
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On a daté la position de l’objet ; ce travail, pénible lors d’une réso- 
lution par le calcul, est chose facile pour l’ordinateur. L'objet avance 
de la distance dL = V (ax)? + (dY)? à la vitesse V pendant l'inter- 
valle de temps dT = dl / V. La valeur de dT s’ajoute au temps précé- 
dent chaque fois que l’objet avance ou recule de dX = 4. 


Liste des variables 


X,Y 
YY 
V 
F 
É 
H 
AC 
S 


EC 
W 
dX 
dY 
dT 
Z 


abscisse et altitude de l’objet 

valeur de Ÿ précédente 

vitesse 

force de frottement 

énergie totale au départ 

= 100/E 

EC / E en %, AP = EP/E en % 

= 1: l’objet va de gauche à droite ; —-1, de droite 
à gauche 

énergie cinétique, EP énergie potentielle 

= E — EP 

progression de l’objet en abscisse 

variation d'altitude quand l’objet avance de dX 
variation de temps quand l’objet avance de dX 
compteur de boucle pour l'affichage. 


Déroulement du programme 


10 

200 

220 

230, 1000 
240 

290 


360 


430 
440 
450 
460, 1000 


préambule 

introduction de la vitesse initiale 

tracé de la route 

valeur de Y 

tracé du profil 

calcul de E, EC, EP initiaux, préparation de 
l'affichage 

début de la boucle 

calcul de EC, EP, effacement et affichage des 
valeurs 

effacement du mobile 

le mobile avance 

l’objet sort de l’écran, fin 

calcul de la nouvelle altitude 
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470 marquage de la nouvelle position 
490 calcul de l’augmentation de temps 
530, 1100 calcul de la nouvelle énergie 

La ligne 580 évite une division par zéro 

600 fin de la boucle 


Commentaire des illustrations 


On a utilisé les lignes 1000 et 1100 du programme ci-joint (la force 
de frottement est constante) avec au départ v = 1. 


Fig. 1 : l'énergie cinétique au départ est très faible, l’ordinateur 
l’arrondit à O0 %. L'énergie potentielle est grande. 

Fig. 2 : l’objet descend la pente, arrive en bas, son énergie poten- 
tielle est alors nulle, son énergie cinétique est grande. 

Fig. 3 : l'objet a perdu de l’énergie en raison des frottements. Son 
énergie mécanique est égale à 11 + 74 = 85 % de l’éner- 
gie initiale, aussi il n’a pas pu franchir la bosse. Après avoir 
fait demi-tour, (il a laissé une trace) il commence à 
redescendre. 

Fig. 4 : l'objet s’immobilise sous l’action des frottements sur la 
pente droite du vallon. Il lui reste encore un peu d’éner- 
gie mécanique. 


DANS LE PROGRAMME, LES LIGNES 1000 ET 1100 SONT DON- 
NEES A TITRE D'EXEMPLE. 


19 CLS:MODE 1 

20 LOCATE 6,8:FRINT'"Conservation de 1 ‘energie 
40 LOCATE 1,22:FRINT'"'pour faire avancer le program 
me, appuyer sur "CHR#(243) 

40 LOCATE 2Z79,25:FRINT CHR#(245) 

00 REM 

69 REM 

70 q=INKEY(1):IF qi:O THEN 70 

80 CLS 

90 FRINT:FRINT'"donner le profil du parcours," 
100 FRINT'y en fonction de x," 

110 FRINT'en lignes 1000 a 1090" 

120 FRINT"(O<x<600 O£y<300) 

30 PRINT:FRINT'"donner la force de frottement," 
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140 
150 
160 


PRINT'"f en fonction de v," 
PRINT''en lignes 1100 a 1190" 
FRINT:PFRINT'"'faire run 200 pour lancer le progr 


amme " 


170 
180 
190 
200 
210 
220 

30 
249 
230 
280 
270 
280 
299 
300 
310 
320 
339 
340 
350 
360 
370 
380 
390 


400 
410 


PRINT:LIST 1000-1090:LIST 1100-1190 
REM 

REM 

MODE 1:CLS: INFUT'vitesse initiale "sv 


CLS 

FOR X=0 TO 600 
GOSUB 1000 
PLOT 10+X,Y,1 
NEXT X 
WINDOW 1,40, 
x=Os:t=0:s=i: 
GOSUB 1000 
Yy=y: ec=0.5xv'"2:ep=9.8*yy :e=ec+ep 
h=100/e 

LOCATE 1,3: PRINT"ZEC=" 

LOCATE 532,3: FRINT"'XEF=" 

LOCATE 1,2:FRINT"t=" 

REM 

REM 

ec=0.5*v"2:ep=9.8*yy 
ac=INT(ec#h):ap=INT (ep#h) 

tt=0. OOLY#INT (1000%t€ ) 

FOR z=3 TO 13:LOCATE z,2:FRINT" ":NEXT z 


1 
dx=4 


LOCATE 3,2:PRINT tt 
FOR z=1 TO 3:LOCATE 4+z,3:FRINT" ":LOCATE 35+z 


SPRINT" MsNEXT z 


420 

30 
440 
450 
460 
470 
480 
490 
900 
910 
920 
330 
940 
9930 
960 
370 
980 
390 
600 


LOCATE S,3:PRINT ac:LOCATE 36,3:PRINT ap 
MOVE 10+x,yy+2:DRAWR 0,8,0 
x=x+dx*s 

IF x>600 OK x<O THEN GOTO 1200 
GOSUR 1000 

MOVE 10+%,y+2: DRAWR 0,8,1 

REM 

dy=y-yy 

dt=SOR (dx*2+dy"2) /v 

t=t+dt 

REM 

GOSUER 1100 
e=<e-f*S0OR(dx"2+dy"2) 

w=e-9,8+%y 

IF w<O THEN s=-s:60T0 360 
v=SOR (2*AES (W) ) 

IF v<i1#107-4 THEN v=i#10"-4 
YY=Y 

GOTO 360 
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610 REM 

990 REM 

1000 y=150+1350xC0S(0.012%x) 

1095 RETURN 

1096 REM 

1100 f=0.5 

1195 RETURN 

1196 REM 

1200 LOCATE 39,1:PRINT CHR#&(245) 
1210 Qq=INÉEY(1):1F q#>O THEN GOTO 1210 
1220 GOTO 200 


7EC= 11 FSEP= ‘74 


t= 218.796 
ZEC= 0 
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t= 0 
rEC-=- 0 FZEP= 99 


t= 23.954 
rEc= 9e rEP= 6 


7. Plan incliné 


Un objet placé sur un plan incliné est soumis à son poids à la réac- 
tion du plan et à une force F. Est-il en équilibre ? Sinon quel est son 
mouvement ? 

L'utilisateur indique à l’ordinateur des valeurs de tous les para- 
mètres. L'ordinateur dessine l’objet sur le plan, les vecteurs représen- 
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tent les forces et il indique les valeurs de ces forces. Eventuellement, 
l'ordinateur indique le sens de déplacement de l’objet. 


Pourquoi ce programme ? 


Un objet ponctuel soumis à plusieurs forces est en équilibre si 
la somme vectorielle de ces forces est nulie. Ce programme montre 
que le principe permet de résoudre des cas complexes comme celui-ci. 


D'autre part nous modélisons le phénomène du frottement sec 
en une ligne de calcul (300), est-ce conforme à la réalité ? En essayant 
le programme dans de nombreux cas, vous conviendrez que ce modèle 
très simple donne des résultats satisfaisants. 

On distinguera bien le principe qui régit l'équilibre que l’on ne peut 
remettre en question sans bouleverser la physique et l'équation qui 
représente le frottement qui n’est qu'une proposition intuitive, approxi- 
mative et qui peut être modifiée au besoin. 


Utilisation du programme 


Les masses sont en kg, les forces en N, les angles en degrés mais 
les données doivent être introduites sans unités. 


Taper ENTER après l'introduction d’une donnée. 
Taper — pour faire avancer le programme. 
Taper ESC ESC pour sortir du programme. 


Pour se familiariser avec le programme, on étudiera d’abord l’objet 
sur un plan horizontal (A = O0) et sans frottement (K = O). 


Le cœur du programme 


L'objet est soumis à son poids P, àlaforce F , à la réaction du 
sol que l’on a décomposé en N et FF (fig. 1). 
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Fig. 1. — Bilan des forces sur l'objet 


Supposons l’objet en équilibre : 
Pon F EN HE 0 
Projetons sur Oy : N = P cosA + F. cos (AF — À + 90°) 
sur Ox: FF + P.sinA = F. cos (AF — A) 
(lignes 260 et 270) 
À, AF, P, F, étant donnés, on en déduit N et FF. 


Une solution qui conduit à N < 0 est à rejeter, elle signifierait que 
le plan tire sur l’objet pour le maintenir en équilibre. Comme il ne peut 
pas le faire, l’objet décolle du plan. 


D'autre part, FF ne peut pas dépasser la valeur K.N (K : coeffi- 
cient de frottement) (ligne 300). 

Si la valeur de FF ne provient pas de l'équation (projection sur 
Ox), mais de FF = K. N, la somme des forces appliquées ne sera pas 
nulle et l’objet se mettra en mouvement en glissant sur le plan. 
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Liste des variables 


A 
AF 
F 
M 
P 
N 
FF 
K 
FM 


angle du plan avec l'horizontale 
angle de la force F avec l'horizontale 
valeur de la force 
masse de l'objet 
poids de l’objet 
composante normale de la réaction du plan 
composante tangentielle de la réaction du plan 
coefficient de frottement 
valeur de la plus grande force 
coefficient d'échelle 
— 1 l’objet glisse vers la gauche 
1 l’objet glisse vers la droite 
2 l’objet est immobile. 


Déroulement du programme 


10 

90 

250 

280, 800 
300 

360 


présentation 

introduction des données 
calcul 

si l’objet décolle 

l’objet glisse 

dessin des forces 


590, 1600rappel des données et résultats 


670 


fin 


680, 1600rappel 


700 


10 CLS:DEG 


retour à 90 


20 LOCATE 14,8:FRINT'"Plan incline" 
30 LOCATE 1,21:PRINT"pour faire avancer le program 


me ,appuyer 


sur "CHR#(243) 


40 LOCATE 39,25:PRINT CHR$# (2435) 


00 REM 
60 REM 


70 q=INKEY(1):1F q<>0 THEN 70 


89 CLS 


90 PRINT:FRINT'masse de l'objet (kg) 7" 
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100 INFUT m 

110 FRINT:PRINT'"angle du plan incline avec 1'‘horiz 
ontale" 

120 PRINT'"(entre © et 89 deg) 7?" 

130 INPUT a 

140 PRINT:PRINT'"'coefficient de frottement 7?" 
150 INFUT k 

160 PFRINT:PRINT'"'angle de la force f avec l'horizon 
tale" 

170 PRINT'"(entre -90 et 90 deg) 7?" 

180 INFUT af 

190 PRINT:PRINT'module de la force ?7 (en N) " 
200 INPUT +f 

210 LOCATE 39,25:PRINT CHR#(243) 

220 q=INFHEY(L):IF qgi>O THEN 220 

230 CLS 

240 REM 

230 p=9.8*m 

260 n=p*xC0S(a)+f*xC0S(af-a+90) 

270 ff=f*CO0S(af-a)-pxSIN(a) 

280 IF n<O THEN GOTO 800 

290 c=2 

300 IF ABS(fF)>k#m THEN CESGN(ff):ff=k#nxc 
310 fm=MAX (p,f,n,ff) 

340 ech=1i60/fm 

350 REM 

360 MOVE 320-180+*C0S (a) ,200-180%SIN(a) 

370 DRAWR 360#C0S (a) ,360XxSIN(a) 

380 TAG:MOVE 316,204: PFRINT CHR# (2:38); 

3970 ORIGIN 320,200 

400 DRAW ech+*f*C0S(af) ,ech#*fxSIN(af) 

410 MOVER 10,10: FRINT'"'F'"'3 

415 ORIGIN 320,200 

420 DRAW ech#nxC0S (a+90) ,ech#nxSIN(a+90) 

450 MOVER 10,10: FRINT'N'; 

440 IF ff=0 THEN GOTO 500 

445 ORIGIN 310,210 

459 DRAW ech*ff#C0S (a+180) ,ech#*ff#*SIN(a+180) 


460 MOVER -c#10,-c#10:PRINT'FF"'3 
500 ORIGIN 320,200 

910 DRAUWR O0,-ech*#p:MOVER 0,-10:FRINT"'P"3 
S20 TAGOFF:ORIGIN O,0 

530 REM 

990 GOSUER 1600 

640 KEM 

650 LOCATE 59,25:PRINT CHR#(245) 
6860 q=INKEY(1):IF q<>0 THEN 660 
6870 CLS 

680 FRINT'essai precedent" 
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6970 GOSUB 1600 

700 GOTO 90 

710 REM 

720 REM 

800 PRINT:PRINT"1 "objet s'envole" 

810 LOCATE 39,2S:FRINT CHR#(243) 

820 q=INKEY(1):1F q<>0 THEN 820 

830 REM 

1590 REM 

1600 FRINT'"m="m:PRINT'"a="as PRINT"k="k:PRINT"'af="af 
sPRINT"'f="+# 

1610 PRINT: PRINT"'p="INT (p):PRINT'#="INT (#) 

1620 FRINT'"'Nn="INT (nn): PRINT"'#fF="INT(f#$) 

1630 IF c=-1 THEN PRINT"I'objet glisse":PRINT" ver 
sS la gauche" 

1640 IF c=1i THEN PRINT'"1'objet glisse":PRINT'vers 
la droite" 

1650 RETURN 


m= à? 
k= 0 N 
af- 9 

£f= 100 
= 98 

= 100 
nz= 98 

£ff=- 6 : 

l’objet lisse 

vers la roite 

P 
ES 

m= 10 
a= @ 
k= 96.5 
af = 45 
f= 100 F 

= 98 
?= 199 
Ffapise lisse à 

Vers ln droite 2 

D, À FF 
P 
ES 
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8. Dynamique 


Votre astronef, dont les moteurs sont en panne, réussira-t-il à 
échapper à l'attraction de la planète Maudite ? Quelle relation y a-t-il 
entre la force appliquée à un point matériel et son mouvement ? C’est 
ce que vous montre ce programme où vous pouvez à tout instant modi- 
fier la force. 


Pourquoi ce programme ? 


On utilise ici de la façon la plus simple qui soit le principe fonda- 
mental de la dynamique, clé de voûte de toute la mécanique. 


Utilisation du programme 


Les masses sont exprimées en kg, les distances en mètres, les 
temps en secondes, les forces de Newton, mais les données doivent 
être introduites sans unité. 

Après introduction d’une donnée, taper ENTER 

Pour faire avancer le programme, taper — 

Pour sortir du programme, faire ESC ESC, puis RUN ou MODE 1. 

Différents types de force peuvent être imposés au point ; soit des 
forces pré-programmées (cas 2, 3, 4), soit une force déterminée par 
l'utilisateur (cas 1). 

Après avoir choisi le cas, taper ENTER. 

Cas 2 : la force est constante de module 20 N et dirigée vers le 
bas. On pourra simuler ainsi une chute libre. 


Cas 3 : la force est dirigée vers le centre de l'écran, son module 
est constant égal à 100 N. La ligne 625 est destinée à éviter une divi- 
sion par zéro. 


Cas 4 : la force est centrale et inversement proportionnelle au carré 
de la distance du point au centre de l'écran. On pourra simuler ainsi 
le mouvement d’un objet sous l’action d’une force de gravitation. 


Dans les cas 2, 3 et 4, on fera avancer le temps d’une seconde 
en appuyant sur la touche. 


Cas 1 : tous les dix tours de boucle, l'utilisateur redéfinit la force 
en donnant d’abord son module puis l’angle (en degrés) de sa direc- 
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tion avec l'horizontale (taper ENTER après l'introduction de chaque 
donnée). Essayez donc de faire suivre au point un trajet fixé à l'avance, 
ce n’est pas toujours facile. 


Le cœur du programme 


Le point évolue sous l’action de la force F. Cette évolution est 
déterminée la relation 


_ 


F = M. A (ligne 320) 

L'ordinateur détermine les composantes AX et AY de l’accéléra- 
tion et en déduit l'augmentation de vitesse par 

dV = À .dt 


dt a été choisi, non pas infiniment petit, mais égal à 0,1 s. L'ordi- 
nateur calcul la nouvelle vitesse égale à V + dV par ses composantes 
VX + AX.dtet VY + AY. dt (ligne 330). 


De la même façon, l'ordinateur détermine la nouvelle position du 
point par: X + VX. dtet Y + VY. dt (ligne 340). 


Dans le cas (1), l'ordinateur effectue dix fois cette boucle de cal- 
cul avant de demander à l'utilisateur de redéfinir la force appliquée 
au point. 


Liste des variables 


t temps 

Z compteur de boucle 

X,Y position du point 

dd distance du point au centre de l'écran 

XX différence entre l’abscisse du point et l’abscisse du point 
et l’abscisse du centre de l'écran 

V module de la vitesse 

an angle de la vitesse avec l'horizontale 


VX,VY composantes de la vitesse 
ax,ay composantes de l'accélération 


C type de force 

f module de la force 

S angle de la force avec l'horizontale 
fx, fy composantes de la force 
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Déroulement du programme 


10 présentation 

90 introduction des données 

240 début de la bouche 

280 sic = 1, on va en 400 : introduction de la force 
300 dessin de la force 

310 sous boucle de z = 1 à 10 


312 vers 500, 600 ou 700 : calcul de la force 
sic = 1, c’est toujours la force introduite en 400 
315 calcul du mouvement 
350 fin de la sous-boucle 
360 fin de la boucle 


llustrations 


Fig. 1 : l'utilisateur choisit la force (taper ENTER après chaque 
valeur). Départ 10, 390, vitesse 40, angle 0, masse 10. 
Fig. 2 : force constante (taper — pour faire avancer le point). 
Départ 10,390, vitesse 20, angle 0, masse 10. 
Fig. 3 : force centrale de module constant. 
Départ 320, 400, vitesse 16, angle 0, masse 50. 
Fig. 4 : force centrale en 1/d°. 
Départ 320, 400, vitesse 6, angle 0, masse 50. 
Plus l'objet est près du centre de l'écran, plus la force est 
grande. 


10 CLS:DEG: MODE 1 

29 LOCATE 12,8: PRINT"Dynami que" 

30 LOCATE 1,22:PRINT'"'pour faire avancer le program 
me,appuyer sur "CHR$(243) 

40 LOCATE 39,25:FRINT CHR$ (243) 

59 REM 

60 REM 

70 q=INKEY(1):1F q<>0 THEN 70 

80 CLS 

90 PRINT:PRINT'"choix de la force" 

100 FRINT:FRINT"i quelconque","2 constante" 

110 PRINT:PRINT"3S centrale","4 centrale en 1/d°2" 
120 FRINT: INPUT c 

130 FRINT:PRINT'"position initiale 7?" 
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140 
150 
160 
170 
180 
185 
190 
200 
210 
220 
230 
235 
240 
2950 
255 
260 
270 
280 
S00 
310 
312 
515 
320 
330 
340 
3930 
360 
370 
400 
430 
440 
450 
460 
900 
910 
920 
600 
610 
625 
630 
640 
650 
700 
710 
720 
725 
730 
740 
750 


INPUT "'x=" 5% 

INPUT" y="3y 
PRINT:PRINT'"'vitesse initale 7?" 
INPUT"v=" 5; v 

INPUT "angle=";an 
vx=v#CO0S (an): vy=v#SIN (an) 
PRINT: INPUT 'masse="; m 

REM 

LOCATE 39,25:PRINT CHR$(245) 
q=INKEY(1):1F q<>0 THEN 220 
CLS: TAG 

IF c=i THEN WINDOW 1,6,1,24 
MOVE x—4,y+4: PRINT CHRK$ (203); 
t=t+i 

1F c=i THEN GOTO 280 

MOVE 620,10:PRINT CHR$(243);3 
q=INKEY(1):1F q<>0 THEN 270 
1F c=i THEN GOSUEB 400 

MOVE x,yY:DRAUWR fx,fy 

FOR z=i TO 10 

IF c<>1 THEN ON (c-1) GOSUR 500,600,700 
fx=f+COS(S):fy=fxSIN(S) 
ax=fx/ms:ay=fy/m 
VX=VX+ax #0. 1: vy=vy+ay#O. 1 
KX=X+VXHO. 1: yY=y+Vvy#0. 1 

NEXT z 

GOTO 240 

REM 

REM'"quelconque" 

TAGOFF: LOCATE 1,5+2#t: INPUT f 
LOCATE 1,6+2*%t:INPUT s:TAG 
fx=f#C0OS(S):fy=fxSIN(S) 
RETURN 

REM'"constante" 

f-20:s5s--90 

RETURN 

REM'"centrale" 

f=100 

XX=Xx-320:IF xx=0 THEN xx=0.01 
S=ATN((y—-200) /xx) 

IF x2320 THEN s5=5+180 

RETURN 

REM'"centrale en i/d*2" 
dd=(x-520)"2+(y-200) 2 
f=(1*10"6)/dd 

KXX=X—S20:1F xx=0 THEN xx=0.01 
Ss=ATN((y-200) /xx) 

IF x>320 THEN s=s+180 

RETURN 
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Fig. 4 


9. Composition de mouvements 


Un autobus parcourt les rues tortueuses de la ville, un passager 
se déplace d’un siège à l’autre dans l’autobus. Quel est le mouvement 
du passager par rapport au sol ? C'est la composition du mouvement 
de l’autobus par rapport au sol et du mouvement du passager par rap- 
port à l’autobus. Ce programme d’animation vous propose quelques 
exemples de composition de mouvements simples. 


Pourquoi ce programme ? 


N’avez-vous pas observé qu’une bille lâchée sur le plancher de 
l’autobus qui prend des virages a un comportement curieux ? Qu'il 
est difficile de suivre une trajectoire simple lorsque vous marchez sur 
le plateau d’un manège qui tourne ? 


On s'intéresse au mouvement par rapport au sol parce que les 
lois de la physique permettent d'expliquer plus facilement le mouve- 
ment d’un objet par rapport à un repère fixe. 


Utilisation du programme 


Vous choisissez d’abord le mouvement de l’objet par rapport au 
rectangle (mouvement relatif), puis le mouvement du rectangle par rap- 
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port au sol (mouvement d'entraînement), enfin vous observez le mou- 
vement de l’objet par rapport au sol. 


Pour faire avancer le programme (quand apparaît par exemple 
le mot ‘‘fin’’), taper — 


Pour choisir un mouvement, taper le chiffre correspondant puis 
ENTER. 


Pour sortir du programme, faire ESC ESC. 


Remarquons que ce programme permet également d'étudier la 
loi de composition des vitesses. La vitesse moyenne entre deux mar- 
ques est égale à la distance entre ces points divisée par l’unité de temps 
de l'horloge. Mesurez les vitesses absolue, d'entraînement, et rela- 
tive pour la même valeur de h ; dans le cas où le rectangle ne tourne 
pas vous pourrez vérifier que 


— 


Vo= V,+ V, 


Le cœur du programme 


Connaissant la position du rectangle par rapport au sol (sommet 
A en xa, ya, le côté AB fait un angle a avec l'horizontale) et la position 
de l’objet par rapport au rectangle (x,, y,), il faut calculer la position 
de l’objet par rapport au sol (x,y). Ceci est réalisé en lignes 4020,4030 
comme le montre la figure suivante (fig. 1). 





! 
x x 


Fig. 1. — Composition des mouvements 
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Liste des variables 


M=1 lors du mouvement absolu 

h horloge 

r choix du mouvement relatif 

d choix du mouvement absolu 

a, Xa, ya position du rectangle 

xr, yr position du point dans le rectangle 
X, Yÿ position du point 


Déroulement du programme 


10 présentation 

90 choix du mouvement relatif 
150 boucle de h = 1 à 25 

3000 calcul de la position relative 
4000 position du point 

210 choix du mouvement d'entraînement 
280 boucle de h=1 à 25 

1000 position du rectangle 

2000 dessin du rectangle 

350 mouvement absolu 

370 boucle de h = 1 à 25 

1000 position du rectangle 

3000 calcul de la position relative 
4000 position du point 

1470 le rectangle est effacé 

999 fin 


Quelques exemples d’utilisation 


La combinaison de deux mouvements rectilignes uniformes con- 
duit à un mouvement rectiligne uniforme (fig. 2). 

La combinaison d’un mouvement rectiligne uniforme et d’un mou- 
vement rectiligne uniformément accéléré conduit à un mouvement 
parabolique (fig. 3a, 3b, 30). 

La combinaison d’un mouvement rectiligne uniforme et d’un mour- 
vement circulaire conduit à un mouvement cycloïdal (fig. 4a, 4b, 4c). 
Ce mouvement est analogue à celui du point d’une roue d’une voi- 
ture se déplaçant à vitesse constante. 
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10 
20 
30 
49 
me , 
90 
69 
70 
89 
90 
100 
110 
120 
130 
140 
1350 
195 
160 
179 
180 
190 
195 
200 
210 


220 
230 
240 
250 
280 
270 
280 
285 
290 
300 
310 
315 
320 
350 
340 
330 
360 
370 
380 
390 
400 
405 
410 
429 
430 
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CLS: INK O,0:INK 1,26 

LOCATE 14,8:PRINT'"COMPOSITION" 

LOCATE 13,12:PRINT"DE MOUVEMENTS" 

LOCATE 1,22:FRINT'"pour faire avancer le program 


appuyer sur "CHR#(243) 

REM 

REM 

q=INKEY(1)3:1F q<>0 THEN 70 
CLS 


PRINT'"Choix du mouvement relatif" 
m=0:xa=0:ya=0:a=0 :e=i 
PRINT:PRINT'"is:immobile","2:uniforme" 
FRINT:PRINT'"3:accelere","4:circulaire" 
PRINT: INPUT r 
CLS: GOSUB 2000 
FOR h=i TO 25 
FOR z=0 TO SOO:NEXT z 
LOCATE 35,1: PRINT'"h="h 
ON r GOSUB 3100,3200,3300,3400 
NEXT h 
LOCATE 1,1:FRINT"#fin" 

REM 
q=INKEY(1):1F q<>0O THEN 200 
CLS:FRINT"Choix du mouvement d'entrainement" 


m=0: xr=0:yÿyr=0 

FRINT:PRINT'"is:uniforme horizontal lent" 
PRINT:PRINT'"2suniforme horizontal rapide" 
PRINT :PRINT'"'3suniforme vertical" 
PRINT:PRINT'"'4:circulaire" 

FRINT:INPUT d 

CLS:FOR h=1i TO 25 

FOR z=0 TO SOO:NEXT z 

LOCATE 35,1:PRINT'"'h="h 


ON d GOSUE 1100,1200,1300,1400 
NEXT h 

LOCATE 1,1 :FRINT"fin" 

REM 

REM 


q=INKEY(1):1F q<>0 THEN 340 
CLS: PRINT' mouvement absolu" 
m=i 


FOR h=1 TO 25 

LOCATE 35,1:FRINT"h="h 

ON d GOSUB 1100,1200,1300,1400 
NEXT h 

LOCATE 1,3:PRINT'"#in" 

REM 

REM 


q=INKEY(1):1F q<>0 THEN 430 


440 CLS:GOTO 90 

450 REM 

999 END 

1000 KEM 

1910 KEM'entrainement" 

1100 a=0:ya=120 

1110 xa=7#%h +10 

1120 GOTO 1435 

1200 a=0:ya=120 

1210 xa=14%*h 

1220 GOTO 1435 

1300 a=0:xa=170 

1310 ya=8*xh+10 

1320 GOTO 1435 

1400 a=0.25%h 

1410 xa=340-170#C0S (a+0.49) 
1420 ya=200-170#SIN(a+0O.49) 
1435 e=1:GOSUB 2000 

1440 IF m=0 THEN RETURN 

1450 ON r GOSUB 5100,5200,3300,:400 
1460 FOR z=0 TO SOO:NEXT z 
1470 e=0:GOSUB 2000 

1480 RETURN 

2000 REM 

2010 REM'"rectangle"“ 

2020 ORIGIN »xa,ya 

2030 s=SIN(a):c=C0S(a) 

2040 DRAWR 300+%c,500%s,e 

2050 DRAWR —-160%s,160%c,e 
2060 DRAWR —-300%c,-1:00+%s,e 
2070 DRAWR 160%s,-160%c,e 
2080 RETURN 

3000 REM 

3010 REM'relatif" 

3100 xr=30:yr=30:GOTO 49000 
3200 xr=12*xh:yr=100: GOTO 4000 
3300 xr=0.48*%hx#h:yr=100::GOTD 4000 
3400 xr=80+60*xC0S (0.252#h) : yr=80+60xSIN(O.252#h):6G 
OTO 4000 

4000 REM 

4010 REM'"marquage du point" 
4015 ORIGIN 0,0 

4020 x=xa+xr*xC0OS (a) -yrxSIN(a) 
4030 y=ya+xrxSIN(a) +yrxC0S(a) 
4040 FLOT »x,y 

4050 RETURN 
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mouvement absolu h= es 
fin 


Fig. 2 
h= es 

Fig. 3a 
h= 25 

Fig. 3b 


nouvement absolu 


fin 


mouvement absolu 


fin 


Fig. 3c 
h= 25 
Fig. 4a 
h= 2e5 
Fig. 4b 
h= 25 
Fig. 4c 
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10. Le pendule simple 


Il est à la base de l’étude des mouvements sinusoïdaux. Rappe- 
lons l’équation différentielle qui régit son mouvement : 
2 
| LE + gsin0 = 0 
dt? 
Si l’on se limite à l'étude des ‘‘petits’”’ mouvements, sin 8 # 8 et 
donc l'équation se simplifie : 
2 
| g + g0=0 
dt 


dont on connaît la solution analytique : 
8 = 60, COS (wt + oo) 


pour laquelle : 


u= Mar = 25) 
T g 


Nous proposons de simuler une telle expérience, et de comparer 
les résultats obtenus à ceux fournis par l'intégration directe de l’équa- 
tion différentielle non simplifiée. Le pendule sera lâché dans une cer- 
taine position, sans vitesse initiale. 


Programme 1 


ll permet de visualiser le mouvement de deux pendules, l’un réel, 
obéissant à l'équation différentielle de départ, l’autre théorique, satis- 
faisant à l'équation simplifiée. Les dessins (1) à (5) montrent divers 
aspects de l'écran à des instants consécutifs. Il s’agit là d’une simu- 
lation du phénomène. 


Programme 2 


Il permet de représenter les élongations angulaires des deux pen- 
dules précédents en fonction du temps. On constate que, si les deux 
pendules sont bien synchrones pour les petits angles initiaux : (6), ils 
ne le sont plus dès que ces mêmes angles ont de grandes valeurs : 
(7) et (8). 
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10 

20 

39 

49 

90 

60 

70 

80 

dro 
90 

100 
110 
120 
130 
140 
130 
160 
170 
180 
190 
200 
210 
220 
230 
240 
230 
260 
270 
280 
290 
500 
:RE 
3510 
RET 
3520 
KET 
3350 
RET 


REM HE EREUEUEEN 
REM *# + 
REM * pendule * 
REM * * 
REM HEURE 
REM 
INFUT'"'longueur du pendule:"3;1 
FRINT'position initiale en degres:":INFUT"(+ a 
ite,- a gauche)";to 
g=9.81 
CLS: #=320:y=550 
t=2*PI#(1/g)"0.S:ti=to 
dt=t/100:t1=t0 
tO=2#PI#tO0/3560 
ti1=2*xFI#t1/360 
a=g/l1:p0=0:pi=0:GOSUR 260 
t=0 
FOR n = Ô TO 478 
dO=-axto:di=-axSIN(t1i) 
tet+dt:LOCATE 2,23:PRINT "tstst 
po=pôo+do*dt:Pi=Fi+dixat 
GOSUR 310:GOSUB 330 


tO=tO+poxdt:ti=ti+PFixdt 
GOSUR 300:GOSUR 3:20 
NEXT n 

STOF 


LOCATE 2,20: FRINT'"1= "sls'"m." 
LOCATE 2,21:FRINT'periode theorique=";t;"s." 
LOCATE 2,22: PFRINT'"'angle de depart=";ti;"degres 


RETURN 

FLOT x,y:DRAWRK 1#19OXxSIN(tO) ,-1#190xCO0S (t0) ,1 
TURN 

FLOT x,y:DRAWR 1#19OXxSIN(tO) ,-1#190xC0S(t0) ,0: 
URN 

FPLOT x,y:DRAWR L#i9OxSIN(EI1) ,-1#190xC0S(t1) ,1: 
URN 

PLOT x,y:DRAWR Ll#19OXxSIN(E1) ,-1#190xC0S(t1) ,0: 
URN 
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Fig. 1 
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periode theori 
degres, 


angle de F 
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1= Lan. 
periode teoriquse 2.08686668 s. : 
angle CT epar 29 degres. Fig. 2 

= La. 

periode theoriquez 2: 2.08686668 s. 

angle de SR degres. 

t= 4.5 

Fig. 3 


62 


=: 
période th theori na ny lis s. 

e de depar egres £ 
16 53977 17187 Fig. 4 
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10 REM HER 

20 REM +% * 

30 REM *% pendule * 

40 REM +%# * 

SO REM HHXRUMUMUONHEUE EX 

60 KEM 

79 INFUT'"longueur du pendule:";1 
80 FRINT'position initiale en degres:": INPUT" (+ a 
droite,- a gauche)";to 

90 g=9.81 

100 MODE 2? 

110 CLS: x=320: y=2:50 
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120 t=2*xFl#(1/g)"0.S:ti=to 

150 dt=t/100:t1=tO 

140 tO=2*FIxt0/760 

150 ti=2#PI#t1/:60 

160 a=g/l1l:p0=0:pi=0:GOSUER 5:00 

170 t=0 

180 FLOT 0,200: DRAW 659,200 

190 FLOT 0,30: DRAW 0,:80 

200 FOR n=0 TO S :FLOT 5,200+(FI1/2-nxF1/4)#100:DRA 
WR -5,0 : NEXT 

210 FOR n = 0 TO 639 

220 dO=-axto:di=-axSIN(t1) 

220 t=t+dt:LOCATE 2,23:FRINT "te"st 

240 pO=po+do*xdt:F1i=Fi+dixdt 

230 tO=tO+pOrxdtitisti+Fixdt 

2609 FLOT n,200 + TO#100 

270 FLOT n,200 + Ti1x100 

280 NEXT n 

290 STOF 

300 LOCATE 2,20:FRINT"I1= “s1l3;"m." 

310 LOCATE 2,21:FRINT'"periode theorique=";t;"s." 
220 LOCATE 2,22:FRINT'"'angle de depart=";ti:"degres 


330 RETURN 





1= 


La. 
periode theorique= 2.808686668 s. 
angle de ASE = 2 degres. 
t= 12.83882 


Fig. 6 
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= La. 

perieat pti à vez, “ 0686668 5. 
angle d degres. 

t= 12. TES 


Fig. 7 





= M. LA AV RNA JAN VA 
periode de depart 69 à 2.08686668 s. 
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degres. 
18.8 


Fig. 8 


11. Le pendule de Foucault 


L'expérience du pendule de Foucault est l’une des plus célèbres 
du dix-neuvième siècle. 


Foucault attribua le phénomène observé à la rotation de la Terre 
autour de l’axe des pôles. 


Le plan d'oscillation d’un pendule simple abandonné sans vitesse 
initiale tourne lentement autour de la verticale du point d'attache du fil. 
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Mise en équations du phénomène 


On se propose d'étudier l'influence de la rotation de la Terre sur 


le mouvement, rapporté à un référentiel lié à la Terre, d’un pendule 
pesant supposé réduit à une particule M de masse m suspendue à 
un fil de longueur | fixé en F. 





X mg 
Z est la latitude du lieu considéré. 


Z est la latitude du lieu considéré. 
Dans ce référentiel (non galiléen), la masse m est soumise : 
— à son poids — mgk (ce qui revient à négliger la force d'inertie 


d'entraînement par rapport à la force d'attraction terrestre) ; 
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— àlatension T dufil; 

— à la force de Coriolis (force d'inertie complémentaire) : 
fc = —- 2mA,V ; 

— À est le vecteur rotation de la Terre. 


Appliquons la relation fondamentale de la Dynamique. 
m = mg + T- 2mA v 
Soit en projetant sur les axes : 


2 
Pie EX sa cos 2 Sn 2%) 
dt? m dt dt 


2 

Passe, Ni az x 

dt? m | dt 

2 
Le = — Q + : (1 — “à + 2Acosz %* 
dt m dt 
Si le mouvement est supposé plan, on prend z = 0, 
d?z dz 

- = Oet = (0. 

dt? dt 
On a donc 0 = - g + L + 2Acosz 

m dt 


Le dernier terme étant négligeable devant g, on peut écrire T= mg. 
Le système d’équation devient : 


2 
dx _ = À + + 2Asinz Ÿ 
2 
PV or 9 re SASnz 2 


Les caractéristiques du pendule monté par Foucault au Panthéon 
à Paris étaient les suivantes : | = 67m, période d'oscillation 
(T = 2PI V lg = 16,421 s, g = 9,81 m/s? et Z = 48951’. 
Remarque : Pour la Terre À = 0,73.10-“rad/s. 


Foucault énonça dans les termes suivants la loi qui relie la durée 
du plan de révolution du plan d'oscillation à la latitude du lieu : 


‘‘A mesure que l’on approche de l'équateur, le plan de l'horizon 
prend sur l’axe de la Terre une position de plus en plus inclinée, et 
la verticale, au lieu de tourner sur elle-même comme au pôle, décrit 
un cône de plus en plus ouvert ; il en résulte un ralentissement dans 
le mouvement du plan d’oscillation, mouvement qui annule à l’équa- 
teur pour changer de sens dans l’autre hémisphère”. 


Le programme suivant permet la résolution numérique des équa- 
tions du mouvement du pendule et calcule la déviation au cours de 
chaque oscillation. 


Le mouvement est simulé à la latitude de Paris, en un lieu de même 
latitude en valeur absolue dans l’hémisphère sud et enfin à l’équa- 
teur où l’on observe aucune déviation. 
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Liste des variables 


A vitesse angulaire de la Terre (73E-4 rad/s) 

G accélération de la pesanteur à la latitude de 
Paris (9.81 m/s°) 

Z latitude de Paris 48.85 (exprimée en degré déci- 
mal) à l'équateur Z = 0 

L longueur du pendule (67 m) 


Déroulement du programme 


lignes 10 à 50 titre et introduction des données 

ligne 70 pas de calcul (dT = 1/100 5) 

lignes 90 et 100 tracé des axes 

lignes 110 à 200 résolution numérique des équations différentiel- 
les du mouvement et tracé sur l’écran du mou- 
vement de l'extrémité du pendule sur le plan 
horizontal. 

ligne 300 l'ordinateur imprime sur l'écran la déviation 
observée au cours d’une oscillation. 


Remarque : 

Nous retrouvons bien les 2.7 mm observés par Foucault au cours 
de son expérience au Panthéon. La rotation du plan d’oscillation à la 
latitude de Paris est d'environ 12 degrés par heure. 


10 REM PENDULE DE FOUCAULT 

20 REM 66e DE DEEE ED DE DEEE NE 

30 REM 

40 REM MOUVEMENT DANS L'HEMISFHERE NORD LATITUDE D 
E FARIS 48 51° 

50 A=0.000075:6=-9.81:2=48.85*FI1/180:L=607 
60 CLS 

70 DT=0.01 

80 MODE 2 

85 LOCATE 32,1:FRINT "FENDULE DE FOUCAUT" 
90 FLOT 329,0: DRAWKR 0,400 

100 FLOT 0,200: DRAWR 640,0 

110 REM CONDITIONS INITIALES 

120 REM #66 EEE DEEE EEE UE NE 

130 X=3:Y=0: TX=0: TY=0 

140 REM EQUATION DE LA TRAJECTOIRE 
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1950 REM He EEE DE EME DEEE DEEE ME DEEE DEEE 

160 FOR F=0 TO 16.421 STEP DT 

170 ATX=-G/LaX+2#AXSINCZ)#TY: ACTX=ATX#DT: TX=TX+ACT 
X2ACX=TXY#DT: X=X+ACX 

180 ATY=-G/LExY-2#AXSINC(Z)XTX: ACTY=ATYY#DT: TY=TY+ACT 
Y2ACY=TYADT: Y=Y+ACY 

190 PLOT X#100+2:20,Y#20000+200 

200 NEXT F 

S00 LOCATE 2,3:PRINT "DEVIATION="3:;Y#1000: "MM' 


PENDULE FOUCAUT 


DEVIATION=-2.70772619 MM 





10 REM PENDULE DE FOUCAULT 
20 REM HR EDEENEMEUEDEDEDEN DEEE DE 


S0 REM 

40 REM MOUVEMENT DANS L'HEMISFHERE SUD LATITUDE - 
48 51° 

90 ÀA=0.000073:6=9.81:27=-48.8S*F1/180:L=67 

60 CLS 

70 DT=0.01 

80 MODE 2 


85 LOCATE S2,1:PRINT "PENDULE DE FOUCAUT" 
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90 FLOT 320,0:DRAWR 0,400 

100 FLOT 9,200:DRAWR 640,0 

110 REM CONDITIONS INITIALES 

120 REM 366626 DD DE DEEE EN 

1530 X=S: Y=O0: TX=0: TY=0 

140 KEM EQUATION DE LA TRAJECTOIRE 

150 REM 3636 DE EMEHE DE DE DE DEN DE DE EDEN DEN 

160 FOR F=0 TO 16.421 STEP DT 

170 ATX=-G/La#X+2#AXxSINC(Z)#TY:ACTX=ATXEDT: TX=TX+ACT 
X=:ACX=TX#DT: X=X+ACX 

180 ATY=-G/L#xY-2RAXSINÇ(Z)#TX: ACTY=ATYEDT: TY=TY+ACT 
Y2:ACY=TY*#DT: Y=Y+ACY 

190 PLOT X#100+2:20,Y#20000+200 

200 NEXT F 

300 LOCATE 2,3:FRINT "DEVIATION="3Y#1000: "MM" 


PENDULE FOUCAUT 


DEVIATION= 2.78772619 MM 





19 REM FENDULE DE FOUCAULT 

20 REM #6 EEE EE DEEE DEEE NE DEEE 

30 REM 

40 REM MOUVEMENT À L'EQUATEUR LATITUDE 0 
50 A=0.000073:6=9.81:7=0:L=67 


70 


60 CLS 

79 DT=0.01 

80 MODE 2 

85 LOCATE 52,1:FRINT "PENDULE DE FOUCAUT" 

90 PLOT 320,0: DRAWR 0,400 

100 PLOT 0,200:DRAWR 640,0 

119 REM CONDITIONS INITIALES 

120 REM #36 62666 De DEEE DE DE DEEE EDEN 

130 X=3:Y=O0: TX=0: TY=0 

140 REM EQUATION DE LA TRAJECTOIRE 

150 REM 36666 MEME DE DEN NE DE DE DE DE DE DE DEEE DE DE DEN 

160 FOR F=0 TO 16.421 STEF DT 

170 ATX=-G/L#X+2#AXSIN(Z)#TY:ACTX=ATX#DT: TX=TX+ACT 
X:ACX=TX#DT: X=X+ACX 

180 ATY=-G/LAxY-2#A4SIN(Z)#TX: ACTY=ATY*#DT: TY=TY+ACT 
Y:ACY=TY*X4DT: Y=Y+ACY 

190 PLOT X#100+:20,Y#20000+200 

200 NEXT F 

300 LOCATE 2,5:PFRINT "DEVIATION="3;Y#1000; "MM" 


PENDULE FOUCAUT 


DEVIATION= @ MA 
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12. Oscillateurs couplés 


La production d’oscillations forcées dans un système oscillant met 
en évidence des forces dont les modules sont des fonctions sinusoi- 
dales du temps. 


Ces forces sont produites par un système oscillant appelé ‘‘exci- 
tateur’’ que l’on fait agir sur le système à étudier appelé ‘‘résonateur”’. 


Si la transmission des forces se fait entre les deux systèmes par 
une liaison élastique, on parle alors d’oscillateurs couplés. 


O1 
(R1) + 
Aî 
(R) 
A2 


(R2) 


O2 





Mise en équations 


Envisageons le cas où les deux oscillateurs sont des ressorts R1 
et R2 liés par un troisième ressort R. Les constantes de raideur sont 
respectivement K1, K2 et K. Les ressorts R1 et R sont liés par une 
masse M1 en A1. Les ressorts R2 et R sont liés par une masse M2 
en A2. 

— appelons | la distance O, O, 


— appelons at, a2 et a les longueurs respectives des ressorts 
R1, R2 et R hors tension. 


Soit Z1 et Z2 les cotes de A1 et A2 en mouvement et Z10 et Z20 
à l'équilibre. 
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Les équations d'équilibre s’écrivent : 
—M1g + K1(l —- ai —- Z10) - K(Z10 - Z20 -— a) = O, 
— M2g - K2(Z20 -— a2) + K(Z10 - Z20 - a) = 0 

Si les forces de frottement fluide sont supposées négligeables, 


c'est-à-dire si l’on considère les ressorts comme des oscillateurs non 
amortis, les équations du mouvement s’écrivent : d71 


— Mig + K1(1 —- Z1 - at) - K(Z1 - Z2 - a) = M1 














dt? 
2 
— M2g -— K2(Z2 - a2) + K(Z1 - Z2 - a) = Me TEE 
Posons  X1 = Z1 - Z10 
et X2 = Z2 — Z20 
On peut écrire : 
2 
M1 ax = — (Ki + K)X1 + KX2 
dt 
2 
M2 en = — (K2 + K)X2 + KX1 


Les exemples présentés traitent les cas concernant un couplage 
faible entre deux oscillateurs identiques (K1 = K2) liés par un res- 
sort de constante de raideur K. (M1 = M2) 

C'est le cas par exemple de deux pendules identiques faiblement 
couplés. 


Remarque : 

L'utilislateur pourra étudier à loisir ses exemples personnes. Nous 
avons aussi introduit dans le programme une force de frottement fluide 
de la force B dx/dt (ligne 270 et ligne 310). 


Liste des variables 


K, K1, K2 constantes de raideur des ressorts R, R1 et R2 

M1, M2 valeurs des masses M1 et M2 

D durée de l'observation (maximale prévue) 

B1 et B2 coefficients correspondant à des forces de frot- 
tement fluide 

XO et YO valeurs de Z10 et Z20 à l'instant initial 


TXO et TYO valeurs de dZ1/dt et dZ2/dt à l'instant initial. 
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Déroulement du programme 


10 
29 
30 
40 
90 
99 
60 
70 
80 
90 
100 
110 
150 
160 
170 
180 
190 
200 
210 
215 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
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lignes 10 à 200 titre et entrées des données (ne pas oublier B1 
et B2 si vous désirez observer le mouvement 
avec un frottement fluide) 

ligne 210 on efface l’écran 

ligne 230 conditions initiales 

lignes 240 à 360 résolution numérique des deux équations dif- 
férentielles : l’ordinateur calcule 600 points 
pour chaque équation (ligne 220). Les variables 
S pour X et R pour Y permettent le calcul d’une 
échelle automatique pour tracer les courbes. 

lignes 400 à 600 on trace les axes, les courbes X = f (T) et 
Y=f(T) et on imprime les données. 


REM OSCILLATEURS COUPLES 
REM HERMIONE DE DE DE DE DEN NE 
REM 
REM ENTREE DES DONNEES 
REM 6H DEDE DE DEN DEN 
CLS 
INPUT "VALEUR DE K=";K 
INPUT "VALEUR DE K1i="3K1 
INPUT "VALEUR DE Mi=";mMi 
INPUT "VALEUR DE K2=";K2 
INPUT "VALEUR DE M2="3; M2 
INPUT "DUREE OBSERVATION D=";D 
REM CONDITIONS INITIALES 
REM HE MED DE MEME DE DE DE DE DEEE 
INFUT "VALEUR DE XO0="3;X0 
INPUT "VALEUR DE TXO="3;TXO 
INFUT "VALEUR DE YO="3; YO 
INPUT "VALEUR DE TYO="3:TYO 
CLS : DT=D/600 
LOCATE 12,12:PRINT "CALCUL EN COURS" 
DIM A(600):DIM B(600) 
X=X0: S=ABS (X0) : Y=YO: R=ABS (YO) s:A(0)=X0:B(0)=YO 
FOR T=0 TO D STEP DT 
N=N+1:LOCATE S,15:FRINT N 
IF N=600 THEN GOTO 355 
ATX=(-(K1+K)#X+K#Y) /M1+HLXTX 
ACTX=ATX#DT: TX=TX+ACTX 
X=X+TX#DT:A(N)=X 
1F ABS(X)>S THEN S=ABS(X) 
ATY=(-(K2+K) #Y+K2#X) /M2+B24TY 


320 
339 
340 
330 
355 
360 
400 
410 
420 
430 
440 
450 
900 
9510 
920 
530 
540 
9930 
9560 
970 
2580 
399 
600 


.V 


ACTY=ATY#DT : TY=TY+ACTY 
Y=Y+TY#DT:B(N)=Y 

IF ABS(Y)>R THEN R=ABS (Y) 
NEXT T 

IF S>=R THEN F=$S 

IF R>S THEN F=R 

MODE 2 

PLOT 0,0:DRAWR 0,400 

PLOT 0,200:DRAWR 640,0 

FOR 1=0 TO N 

PLOT I,A(1)/F#199+200:PLOT 1,B(1)/F#199+200 
NEXT I 

LOCATE 15,17:PRINT "K="35K 
LOCATE 15,18:PRINT "Ki=";K1 
LOCATE 15,19:FRINT "K2=":K2 
LOUATE 15,20:PRINT “"Mi="3; M1 
LOCATE 15,21:PRINT "M2="3;M2 
LOCATE 15,22:FRINT “D=";D 
LOCATE 15,23:PRINT "DT=";DT 
LOCATE 35,17:PRINT "XO="3: XO 
LOCATE 35,18:PRINT "TXO="; TXO 
LOCATE 35,19:PRINT “YO="3 YO 
LOCATE 35,20:PRINT "TYO="; TYO 








ri 





Ê PO ge : . 
Se CR 
” CNE vd . . = 8° % so. 2: Te +! 
; ME Là NN  F-'e : JU Le ea ea Eos 
Nè= 10 = St D, nr fie Æ KA tes. MO 
-.b= 9% ce 4 Nes ; .! .! : : ee : 
DT: 815. 4 HR & 
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2 
Electricité 


1. Champ et potentiel créés par un système 
de charges ponctuelles 


I. Ce programme intéresse aussi bien la mécanique que l’élec- 
tricité. En effet, les équations fournissant les composantes du champ 
électrique produit par une charge ponctuelle sont de même nature que 
celles fournissant les composantes du champ gravitationnel produit 
par une masse ponctuelle. Le traitement de ce genre de problème est 
donc le même dans les deux cas : cependant, ce n’est qu’en électri- 
cité qu’il existe des charges négatives ; en mécanique, toutes les mas- 
ses sont positives. 

Ce programme propose donc d’esquisser les lignes de champ et 
les équipotentielles du champ produit par un système de charges (ou 
de masses) ponctuelles placées sur l’écran. Seules sont esquissées 
en chaque point du plan (l'écran) les tangentes à ces courbes. 

Rappelons que : 

— la force produite par une charge électrique est proportionnelle 
au champ électrique qu’elle crée. Les lignes de champ sont tangen- 
tes en tout point au champ électrique ; 

— une équipotentielle est une ligne perpendiculaire en chaque 
point à la ligne de champ qui y passe. 

Il. Le programme est simple : une fois définies les positions et 
les valeurs des charges sur l’écran, il faut calculer la distance du point 
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courant à chacune des charges, évaluer les composantes sur OX et 
OY des champs électriques produits par les charges en ce point, en 
déduire les composantes du champ résultant, calculer le module du 
champ résultant, puis la tangente de la ligne de champ en ce même 
point (d’où la longueur des calculs). On trace alors un petit segment 
de droite, dont la pente est donnée par la valeur de la tangente ainsi 
calculée. Un certain nombre de points sont ainsi examinés sur l’écran. 
Le tracé des équipotentielles se fait de la même manière, dans la 
mesure où l’on sait que, par définition, en chaque point, une équipo- 
tentielle est perpendiculaire à la ligne de champ. 

Il. Le programme ne présente pas de difficultés majeures (jusqu’à 
la ligne 240). A partir de là intervient un problème lié à l'esthétique 
du dessin. Nous nous sommes imposé de tracer des petits segments 
(symbolisant des petits morceaux de ligne de champ) de longueur cons- 
tante égale à 9 pixels. En effet, pour une pente donnée, si l’on impose 
AX, il est possible que, suivant la valeur de la tangente, YŸ devienne 
très grand puisque AY = AX.tga (dessin 1) s'imposant, on a donc : 


A? = (AXŸ + (AY), on connaît B = - = tga 
A 


D'où, pour un point de coordonnées X et Y, les composantes du 
petit segment que l'on y trace : 


2 
ie 22 AY = B.AX 
B°+1 

On voit qu'il n’y a plus de problème même si B tend vers l'infini. 

En outre, pour le tracé des équipotentielles, le calcul reste le 
même, sauf lignes 260 à 270. En effet (dessin 2), on peut tracer un 
segment de droite perpendiculaire au précédent enutilisant comme 
composantes : 


AX' = AY et AY' = — AX 


Ay 
B = tga 





Fig. 1 
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(Ay)' 





Liste des variables 


X1, X2, Y1, Y2 positions des deux charges 


A2 

Q1 et Q2 
Yet X 
R1 et R2 


carré de la longueur d’un segment (= 9 x 9ici) 
valeurs des charges électriques 

position du point courant 

distances du point courant aux deux charges 
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EX, EY 
FX, FY 
EX, EY 
E 
B 


DX, DY 


composantes du champ produit en X, Ÿ par la 
première charge 

composantes du champ produit en X, Y par la 
seconde charge 

composantes du champ résultant 

module du champ résultant 

tangente de l'angle entre l'horizontale et le 
champ résultant 

composantes (AX, AY) du petit segment. 


Détail du programme 


80 
90 


170 - 180 
190 - 200 


210 - 220 
230 
240 


250 - 260 
270 


ti 


} 
") REM # 


position des deux charges 

valeur des deux charges et longueur des petits 
segments 

définissent les points qui vont être examinés 
évaluation de la distance du point courant à la 
première charge ; calcul du champ produit par 
celle-ci 

idem pour la seconde charge 

calcul des coordonnées du champ total 
calcul de la tangente du champ par rapport à 
l'horizontale 

calcul de AX et AY (voir plus haut) 
positionnement du point et tracé du petit seg- 
ment de droite 


Un OREM  ctdete e e Ae e AE EE AE AE  MEE 


4: 


30 REM # champ electrique # 


40 KREM # 


+ 


M0) REM HA de AA AE AE AE A AE AE 


49 REM 
70 MODE € 


RO x1i=HOO0+x2=400: y 1er yE#=E00 





90 as 


100 FRINT #£& 
110 PRINT #5 
120 PRINT #5 : 
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un reseau 


2. Circuit électrique 
. pour vous présenter une méthode de résolution par approxi- 


Soulevez le coin d’une nappe, quelle forme prend-elle ? Chauf- 
mations successives. Cette méthode est applicable immédiatement 


fez le bord d’une plaque, quelle est la température en ses différents 


points ? Appliquez une tension électrique entre deux points d’ 
du réseau ? Ces problèmes sont voisins, c’est le dernier exemple cité 


de résistances électriques, quel est le potentiel en chacun des nœuds 
que nous avons choisi de traiter ici. 


Pourquoi ce programme ? 
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à des réseaux très étendus. Si nous nous sommes contentés ici d'un 
réseau à six nœuds, c’est pour que la résolution soit rapide, facile- 
ment vérifiable ‘‘à la main”, pour pouvoir dessiner facilement le réseau 
et indiquer les valeurs des résistances et des potentiels. Mais le pro- 
gramme ne serait pas plus long à écrire pour un réseau de six cents 
nœuds. Enfin elle présente l’avantage de ne comporter qu’une ou deux 
lignes de calcul. 


Cette méthode pose le problème du compromis rapidité-précision : 
si vous désirez des résultats dix fois plus précis, les calculs seront dix 
fois plus longs. 


Utilisation du programme 


Vous donnerez les valeurs des résistances en Ohms, sans intro- 
duire l'unité, par un nombre entier compris entre 1 et 999. Vous pou- 
vez également envisager des résistances nulles ou des circuits ouverts. 
Le programme assimile les court-circuits à des résistances de un 
micro-Ohm. 

La tension appliquée entre les bornes du circuit est de 100 Volts. 
Les résultats sont donnés à un volt près. Pour obtenir une plus grande 
précision, modifiez les lignes 380 et 520. 


Le cœur du programme 


Deux relations très simples sont utilisées (fig. 1) : 
La loi d'Ohm : U,-U=R,., 
ou en posant G,;=1/R, 


1, =G,.(U, -U) 

La loi des nœuds : 1, +1,+1,+1,=0O 

soit : 

G;.(U, — U) + G,.(U, — U) + G3.(U3 — U) + G;.(U, — U) = O 
d’où : 


G..U, + G:.U; + G;.U; + Gs.U4 = U.(G; + G: + G; + Gi) 
Dans le programme, on a noté : G=G., +G,+G;,+G,(ligne 360) 
d'où : U={(G.;.U, + G.U, + G;.U; + G,.U,)/G (ligne 370) 
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Fig. 1 


Au départ, tous les potentiels sont fixés arbitrairement à 50 V. Puis 
on applique les potentiels O et 100 V aux bornes de gauche. On se 
propose de calculer le potentiel de chaque nœud au moyen de la rela- 
tion précédente, mais on ne peut pas les calculer simultanément. Le 
calcul du potentiel en un point fait intervenir les potentiels voisins qui 
sont erronés. Les premiers résultats obtenus sont donc erronés. 


Après avoir fait un tour de réseau, on en refait un autre avec des 
valeurs un peu moins erronées que les valeurs arbitraires de départ... 
et ainsi de suite jusqu’à ce que deux tours de réseau fournissent des 
résultats pratiquement identiques (ligne 380). 


Voici ce que cela donne sur un exemple simple (fig. 2) : 


premier calcul : G = 1 + + = 1,25 


U = (100 ! + 50 | 1,25 = 60 
4 1 


100 100 100 100. 100 100 100 


50 50 60 44 37,6 35 34 33,3 
10 44 4 
50 30 22 18,8 17,5 17 16,6 


_ e ue 
0 0 0 0 0 0 0 


8 


——— = © —— | 
initialement on applique 100 V pratiquement en réalité 
identiques 
Fig. 2 
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l'est très agréable de n'avoir qu'une ligne de calcul, malheureu- 
sement tous les nœuds du réseau ne sont pas dans la même situa- 
tion. Il y a des points marginaux (fig. 3) et plus le réseau est petit, plus 
la proportion de points marginaux est importante. 





O nœud marginal 
X nœud normal (4 voisins) 


Circuit électrique Fig. 3 


Nous avons donc été amenés à créer des nœuds fantômes et des 
résistances fantômes. Ces résistances sont infinies (ligne 260), afin 
que le potentiel des nœuds fantômes n'aient aucun effet sur le réseau 
réel. La figure 4 vous montre la notation utilisée. 





Fig. 4 


| 
l 
! 
à 3,7 & 57 


En italiques : indices des conductances 
En caractères droits : indices des nœuds 
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Liste des variables 


A < 


1<CONTXN 
Il 


Il 
O 


compteur de boucle 
indices 

résistance 

conductance 

si la conductance est nulle 
nouveau potentiel 

ancien potentiel 

si le calcul est terminé 


Déroulement du programme 


10 présentation 

90 dessin du circuit 

120 introduction et affichage des valeurs des résistances 

260 conductances fantômes 

290 initialisation des potentiels 

340 calcul des potentiels 

380 le calcul est-il terminé ? 

430 les potentiels calculés deviennent ‘‘anciens”” 

680 en route pour un nouveau tour de piste 

500 affichage des résultats 

Illustrations 

Fig. 5 : toutes les résistances sont identiques. 

Fig. 6 : un cas très simple... 

Fig. 7 : le potentiel des nœuds supérieurs est divisé par deux à 
chaque cellule 

Fig. 8 : le potentiel des nœuds supérieurs est diminué de vingt 
volts à chaque cellule 

Fig. 9 : les nœuds de droite sont isolés du circuit. Ils sont restés 
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au potentiel arbitraire de 50 volts. 


10 CLS 

20 LOCATE 12,6: FRINT"'CIRCOUIT ELECTRIQUE" 

350 LOCATE 10,12:FRINT'"court-circuit taper: -" 

40 LOCATE 10,14:FRINT'"'circeuit ouvert taper: /" 

90 LOCATE 1i,22:FRINT'"pour faire avancer le program 
me, taper'"CHR$(2435) 

60 LOCATE 39,25: PRINT CHR#(2435) 

70 q=INKEY(1):1F qf>0 THEN 70 

80 CLS: TAG 

90 MOVE 48,54: DRAWR 568,0:MOVE 48,328:DRAWR 568,0 
100 FOR z=0 TO 2:MOVE 232+192%#z,54: DRAWR 0,274: NEX 
T z 

110 DIM g(9,9):DIM u(9,9):DIM v(9,9) 

120 FOR X=2 TO 6 STEF 2 

130 FOK Y=3 TO 5 

140 MOVE -64+96%X,576-96%Y:PRINT"?"; 

150 MOVE O0,400:INFUT r+# 

160 IF r$="/" THEN G(X,Y)=0:C=1: GOTO 200 

170 IF R$="-" THEN G(X,Y)=1000000:GOTO 225 

180 C=0:G(X,Y)=1/VAL (RS) 

200 FOR Z=0 TO 2 

210 IF Y<:4 THEN MOVE —-64+9bxx+16%2z,747-137#y:FRIN 
T CHR#(143-15#C); ELSE MOVE 22+96%x,202-16%z2:FRINT 
CHR#(143-15%0C) ; 

220 NEXT Z 

225 MOVE -64+96#%X,576—-96XY:FRINT r#: 

230 NEXT Y 

240 NEXT X 

245 FOR z=0 TO 624 STEF 16:MOVE z,400:FRINT" "32:NE 
XT z 

250 REM 

260 FOR X=3 TO 7 STEF 2:6(X,2)=0:G(X,6)=0:NEXT X 
270 G(8,23)=0:G6(8,5)=0 

280 REM 

290 FOR X=1 TO 9:FOR Y=1 TO 9:U(X,Y)=50:V(X,Y)=50: 
NEXT Y:NEXT X 

300 U(1,3)=100:U(1,5)=0: MOVE 16,384: FRINT"100V'"3;:M 
OVE 16,232:PRINT"O V'": 

310 REM 

330 T=i 

340 FOR X=3 TO 7 STEF 2 

350 FOR Y=3 TO S STEF 2 

260 G=G(X-1,Y)+G(X-1,Y+1)+G(X-1,Y-1)+G(X+1,Y) 

370 UCX,Y)=(U(X-2,Y)#G(X-1,Y)+U(CX,Y+2)#G(X-1,Y+1) + 
U(X,Y-2)#G(X-1,Y-1)+U(X+2,Y)#G(X+1,Y))/G 

380 IF ABS(U(X,Y)-V(X,Y))5:0.2 THEN T=0 

390 NEXT Y 

400 NEXT X 

410 REM 

420 IF T=i THEN GOTO 500 
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450 FOR X=3 TO 7 STEP 
440 FOR Y=3 TO 5 STEF 
450 V(X,Y)=U(X,Y) 

460 NEXT Y 

470 NEXT X 

480 GOTO 330 

490 REM 

500 FOR X=3 TO 7 STEF 2 

510 FOR Y=3 TO 5 STEF 2 

520 MOVE -96+96%x,902-172%y:FRINT STR# CINT (u (x ,y)) 
)+"v"z 

530 NEXT Y 

540 NEXT X 

550 MOVE 16,620:FRINT CHR (247); 


SES) 


S60 q=INFEY(1):1F qf:0 THEN 560 
365 CALL &AOGO00 
97/70 END 


la ligne 565 est destinée à l'imprimante 


100U 63v S4v Sev 





e U 3 7v Au A49v 
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186U 939v 39 3S93v 





e uv av Ov 493v 
1868U SO v e5v 13v 
16 16 10 
e6 20 16 
e u ; av - Ov | Gv 
1008U 86üv 6Ov 40v 
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106U 7Svu sov SOv 


16060 


3608 





e vu Au sou SsOu 


3. Circuits électriques : droites de charge 


Ce programme propose l'étude de circuits électriques simples, 
contenant un générateur et un récepteur. Celui-ci peut être, au choix 
de l’utilisateur : un résistor, une batterie ou une diode Zener. 

Le problème suivant est posé : les éléments du générateur (force 
électromotrice et résistance interne) étant connus, ainsi que ceux du 
récepteur, quel est le point de fonctionnement de l’ensemble ? Le pro- 
gramme teste alors les résultats introduits au clavier, vérifie leur vali- 
dité (à + 10 % près) et en cas d'erreur, fournit les bonnes valeurs de 
U et de I. 

Puis, à titre de vérification, les caractéristiques des deux dipôles 
sont tracées ; leur intersection fournit le point de fonctionnement de 
l’ensemble. 

Toutes les quantités numériques proposées sont absolument aléa- 
toires, ce qui assure des problèmes perpétuellement nouveaux. 


Programme 


Il se compose de plusieurs parties : 
e le tableau de présentation (80 - 170), 
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e le choix aléatoire des éléments (190 - 260), 
+ le calcul des éléments (280 - 800) et de la solution U, |: 
- générateur (280 - 370), 
- résistor (420 - 470), 
- diode Zener (480 - 550), 
- batterie (560 - 680). 
+ l'introduction de la solution (690 - 710) et sa vérification (720 
- 730), 
+ l'évaluation de son exactitude et l'affichage de la solution pro- 
posée par l'ordinateur (770 - 790), 
e le dessin des axes sur l'écran et leurs graduations (820 - 930), 
+ le dessin de la caractéristique du générateur (940), 
le dessin de la caractéristique du résistor (960 - 1010), 
le dessin de la caractéristique de la diode Zener (1020 - 1080), 


e le dessin de la caractéristique de la batterie (1090 - 1150) et 
sa prolongation côté | <0 si par hasard V> E (1170 - 1260). 


Quelques remarques concernant les calculs 


e Choix aléatoire des éléments ; on a ; 


10 <E < 50 E « vz < Æ 
3 3 
100 < R < 1000 10 < RZ < 20 


F<R < À +8/3"RND 10 < V < 50 


100 < RP < 1000 
Toutes ces valeurs sont bien sûr modifiables. A ce propos, il est 
prévu une sécurité pour VZ > E (ligne 520). 

e Les échelles sont telles que : 

- sur OX, axe des courants | : 88 pixels, 0,1 ampère ; 

- sur OY, axe des tensions U : 60 pixels, 10 volts. 
e Une droite de pente R se construit alors de la manière suivante : 
R = dU _ dV6 

di dX/800 


95 


Donc, pour dX=1 pixel, il faut avoir 


SR 


880 


D'où le tracé (lignes 980, 1050...). 
e Dans le cas où V < E, le point d’intersection des caractéristi- 


dY (en pixels) 


ques est à gauche de l’axe des ordonnées (! < 0) ; celles-ci sont pro- 
longées en conséquence (lignes 1170 à 1260). 


Définition des variables introduites dans le texte et dans le 


programme. 


10 
20 
50 
40 
90 
60 
79 
80 
90 





La 


R; 
| | 
batterie , | 
dipôle générateur resistor diode batterie 


zener 
SR —_——————— 


dipôles récepteurs 


REM #66 ee DD DE DEEE DEEE DEEE EE EEE 


REM *% % 
REM * point de fonctionnement +* 
REM # * 
REM HE DE DE DEMI DEDE DEEE DEN 
REM 
CLS 


FRINT:FRINT'"'point de fonctionnement" 
FRINT:FRINT'"-generateur " 


100 FRINT:FRINT'"il1l est defini par sa resistance in 
terne r et sa f.e.m e“ 

110 FRINT:FRINT'"-recepteur" 

120 PFRINT:FRINT'vous avez le choix entre:" 


1270 


FRINT:FRINT'resistor ,de resistance r....,..4(1)" 


140 FRINT'"diode zener ,de tension vz,r... (2)" 
150 FRINT'batterie,f.e.m e,resistance r.(3)" 
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160 
170 
180 
1970 
200 
210 
220 
250 
240 
230 
260 
270 
280 
290 
300 
510 
320 
S350 
340 
550 
360 
370 
380 
390 
400 
410 
20 
450 
440 
450 
460 
LL 
470 
480 
490 
300 
9310 
920 
929 
930 
940 
990 
960 
370 
9380 
93970 
600 
610 
620 


50 


FRINT:FRINT'quel recepteur choisissez vous?" 
INFUT a 

CLS 

REM choix des elements 
FRINT''appuyez 4 fois sur la touche a" 
FOR n=i TO 4 

x=RND (1) 

a$f=INKEYS 

IF a$f<> "a" AND af <> "a" THEN 220 
a (n)=x 

NEXT n 

CLS 

REM generateur 

e=a(1)*x50 

IF e < 10 THEN e=e+i0 

e=INT(e) 

r=a (2) #1000 

IF Fr < 100 THEN r=r+100 

r=INT(r) 

PFRINT:FRINT'"generateur 

FRINT:FRINT" f.e.m='";es"v." 
FRINT'resistance interne="s;r3:"ohm." 
REM recepteurs 

FRINT:FRINT "recepteur " 

ON a GOSUB 420,480,560,670 

GOTO 690 

REM resistor 

ri=r/2:+8%a(3) /3irisINT (r1) 
u=rixe/(r+r1) 

i=e/(r+r1i) 

FRINT:FRINT'resistor de resistance :";r1s"ohm. 


RETURN 

REM diode zener 
vzze/3texa(z)/3:vz=INT(vz) 

rz =10+10%#a(4):rz=INT(rz) 
u=vzt+rzk(e-vz)/(r+rz) 
iz(e-vz)/(r+rz)siF 14 O©O THEN i=zo:uze 
FRINT:FRINT'"diode zener" 
FRINT'tension de seuils:"svzs"v," 
FRINT'"resistance:":;rzs"ohm. 
RETURN 

REM batterie 

v=a (35) #50 

IF v <10 THEN v=v+i0 

V=INT (v) 

rp=a(35)*1000 

IF rp#100 THEN rp=rp+100 

rp=INT (rp) 

i=(e-v)/(r+rp) 
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640 u=e* rp/(r+rp)+vxr/(r+rp) 

650 FRINT:FRINT'"batterie" 

660 FRINT" f.e.ms: “"svs"v." 

670 FRINT" resistance interne :";rp3s'"ohm." 

680 RETURN 

690 PRINT:FRINT'"'donnez moi votre solution" 

700 FRINT:INFUT'"'u=";us 

710 INPUT'"i=";is:FRINT 

720 IF us<1i.1x#u AND usr0.9xu THEN GOTO 730 ELSE G0Q 
TO 750 

7350 IF ABS (is)<1.1*AERS(i) AND ARS(is):0.9*ARS(i) T 
HEN GOTO 740 ELSE GOTO 750 

740 PKINT'"'bonne reponse":6G0OTO 800 

750 FRINT'mauvaise reponse" 

760 PRINT 

770 PRINT'"il fallait repondre: '":PRINT 

780 FRINT'u="353:FRINT USING'"##.#"sus:FRINT"V." 

790 PRINT'"'i="3;s5:FRINT USING'"'##H.##"35:15:FRINT" a." 
800 z$=INKEYS:1F z#$="" THEN 800 

810 MODE 2 

820 PLOT 40,40:DRAWR 0,320 

830 FLOT 40,40: DRAWR 560,0 

840 PFLOT 40,320: DRAWR —-4,-4: PLOT 560,40 : DRAWR -—4, 
4 

850 FLOT 40,320 :DRAWR 4,-4:FLOT 560,40 :DRAWR -4, 
—4 

860 REM graduation o©oyÿy : 1 graduation/ 10 volts . 
870 FOR j = 1 TO 4 : FLOT 28 , 60 x j + 40 : DRKAWR 
4,0 : NEXT j 

880 LOCATE 7,5 : FRINT"'U" 

890 LOCATE 7,8:FRINT"40Q V'" 

900 REM graduation ox : 1 graduation/ 0.1 ampere . 
910 FOR j=1 TO 5 : FLOT 42 + 880 * O.1 % 3j ,40 : D 
RAWR O,-4 : NEXT j 

920 LOCATE 69,22:FRINT'"I" 

950 LOCATE 60,24: PRINT"O0.S5 A " 

940 PLOT 40 , 6*e +40:DRAW 40+880%e/r,40 

950 ON a GOTO 960,1020,1090 

960 REM resistor 

970 FLOT 40,40 

980 x=x+i1 5: y= y+3#r1/440 

990 IF x > 560 OR y > 320 THEN 1270 

1000 DRAWR %»,y 

1010 GOTO 980 

1020 KEM diode zener 

1030 x=0:y=0 

1040 FLOT 40 , 40+6%vz 

1090 x=x+1 : y = y+3xrz/440 

1060 IF x > 560 OR y 520 THEN 1270 

1070 DRAUR %x,y 
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1080 GOTO 1050 

1090 KREM batterie 

1100 x=0 s: y=0 

1110 FLOT 40 , 40 + 6%v 

1120 x=x+1: y=y+3#rp/440 

1120 IF x % 560 OK y > 220 THEN 1160 
1140 DRAWR *x,y 

1150 GOTO 1120 

1160 IF e > v THEN GOTO 1270 

1170 x=02:y=O: FLOT 40 , 40 + 6%v 

1180 x=x-1: y=y-5S+rp/440 

1190 IF 40+x < © OR y+éxv + ©Q THEN 1220 
1200 DRAWR —-1 , —-3xrp/440 

1210 GOTO 1180 

X=0:y=0 : FLOT 40,40+é6+%e 

É K=x—1:y=y+5i#r/440 

1240 IF 40+x < © OK 40+y 5% 399 THEN 1270 
1290 DRAUR -1,3%#r/440 

1260 GOTO 1250 

1270 STOF 





point de fonctionnement 

-generateur 

il est defini par sa resistance interne r et sa f.e.n e 
-recepteur 

vous avez le choix entre: 

resistor de resistance r...... (1) 

diode zener,de tension vz,r...(8) 

batterie,f.e.n e,resistance r.(3) 

quil recepteur choisissez vous? 


generateur 


f.e.n= 19 v 
resistance Taternez 178 ohn. 


recepteur 
resistor de resistance : 68 oha. 
donnez noi votre solution 


u=? 4 
i=2 8.88 


nauvaise reponse 
il fallait repondre: 


u= 4.8v. 
i= 9.68 a. 


99 





point de fonctionnement 

-generateur 

il est defini par sa resistance interne r et sa f.e.n e 
-recepteur 

vous avez le choix entre: 

resistor,de resistance r..... .(1) 

diode sr de tension vz,r...(2) 

batterie,f.e.n e, resistance r.(3) 


gs! recepteur choisissez vous? 


generateur 


f.e.n= 41 v. 
resistance interne 664 oha. 


recepteur 

ension de sul ou v, 
resistance: 19 oha 
donnez noi votre solution 
u=® 85 

i=? 6.823 


bonne reponse 
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point de fonctionnement 
-generateur 

il est defini par sa resistance interne r et sa f.e.n e 
-recepteur 

vous avez Le choix entre: 
resistor,de resistance r...... { 


diode zener ,de tension vz,r... 


1) 
8) 
batterie,f.e.n e,resistance r.(3) 


gus! recepteur choisissez vous? 


generateur 


f.e.n= 3% v 
resistance Taternez 478 oha. 


recepteur 
batterie 

f.e,n: y. 

resistance interne : 587 oha. 
donnez noi votre solution 


u=? 29.5 
i=? 8.081 


bonne reponse 


point de fonctionnenent 

-generateur 

il est defini par sa resistance interne r et sa f.e.n e 
-recepteur 

vous avez le choix entre: 

resistor,de resistance r......(1) 

diode zener ,de tension vz,r..,(2) 

batterie,f.e.n e,resistance r.(3) 


gs! recepteur choisissez vous? 


8.54 
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generateur 


f.e.n= 16 v. 
resistance interne- 544 ohn. 


recepteur 
batterie 

f.e.n: V. 

resistance interne : 514 ohn. 


donnez noi votre solution 


u=? 28.6 
iz2 8.0885 


bonne reponse 





4. Résonance 


On étudie la résonance en tension aux bornes du condensateur 
d’un circuit R L C soumis à une tension sinusoïdale de un volt. L'ordi- 
nateur trace les courbes donnant la tension aux bornes du condensa- 
teur en fonction de la fréquence, et cela pour diverses valeurs de la 
résistance, ainsi que le lieu des points correspondant à la fréquence 
propre du circuit. Après quoi l’utilisateur donne les valeurs des élé- 
ments du circuit et la fréquence de l'excitation. L'ordinateur place alors 
le point correspondant sur le graphique. 


Pourquoi ce programme ? 


Pour la forme : la représentation d’une famille de courbes en pers- 
pective est toujours attractive et parle mieux à l'imagination. 
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Pour le fond : la résonance est un phénomène très courant en 
mécanique, en électricité : pourquoi les vitres du vieil autobus 
tremblent-elles quand le moteur tourne au ralenti ? Pourquoi la bou- 
teille dans laquelle vous soufflez émet-elle un son ? Pourquoi pouvez- 
vous sélectionner avec votre poste un émetteur de radio ? Pourquoi 
faites-vous déborder le café de la tasse quand vous agjitez la cuillère 
à une certaine vitesse ?... Parce qu'il y a résonance. 


Nous étudions ici un système électrique, mais les résultats sont 
valables pour un système mécanique avec la correspondance 


suivante : 
Fe cé | 


Fig. 1 


Utilisation du programme 


Les données sont introduites en Ohm, Henry, microFarad, Hertz, 
mais sans indiquer l’unité, suivies de ENTER. 


Pour faire avancer le programme, taper — 
Pour sortir du programme, faites ESC ESC puis RUN ou MODE 1 


Pour voir sur l’écran le point correspondant au circuit que vous 
proposez, il faut que la fréquence soit inférieure à deux fois la fréquence 
propre Fo et que la résistance soit comprise entre 3,75 et 0,023 fois 
la résistance caractéristique Ro du circuit. 


Le cœur du programme 


L'ordinateur trace les courbes en coordonnées réduites, ainsi 
quelle que soit la fréquence de résonance, la courbe représentant le 
circuit se situe à la même place. 


La valeur remarquable de la fréquence est la fréquence de réso- 
nance en courant Fo, qui correspond à la pulsation «, telle que : 


L.CuË = 1 
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La valeur remarquable de la résistance est la résistance critique 


Ro telle que : 


R, = 2.F. V L/C 


La fréquence réduite est donc F / Fo = w/wo = fr 
La résistance réduite est donc R / Ro = rr 
Soit | l'intensité complexe dans le circuit, Zc l'impédance com- 


plexe du condensateur et Z l’impédance complexe totale du circuit : 


Ue = 2.1 Us =Zc.l Us = Ue.Zc/Z 

Z=R + j(Lw-1/Cw) et Zc = -— j/Cw 

En multipliant haut et bas par jCw, on obtient : 

Us = Ue. 1 / (RCw -— LCu? + 1) 

Us = Ue .1(2j:rr.fr + 1 — fr?) 

En module : 

Us = Ue.1/F. V 4rr2.fr2 + (1 — fr.) (ligne 1010) 





Liste des variables 


n numéro de la courbe 

rr résistance réduite 

fr fréquence réduite 

L inductance de la bobine 

C capacité du condensateur 

Fo fréquence de résonance en courant 
Ro résistance critique 

X,Y coordonnées du point 

u tension aux bornes du condensateur 
a,b coordonnées de la case texte effacée 


Déroulement du programme 
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10 présentation 

82 préparation du graphique 

95 boucle des 20 courbes 

1000 calcul des coordonnées du point 
1020 et 1030 effet de perspective 

180 tracé de la courbe 


1000 calcul des coordonnées du point 


250 introduction des données 

280 calcul de Fo et Ro 

340 calcul de fr et rr 

1000 calcul des coordonnées du point 
360 affichage des résultats 

400 effacement des résultats 


Commentaire des illustrations 


Dans le cas de la figure 2, la fréquence imposée est égale à la 
fréquence propre du circuit mais l'amortissement est élevé si bien que 
la tension aux bornes du condensateur est inférieure à la tension appli- 
quée. |l en est de même pour la figure 3 bien que l’amortissement soit 
faible parce que la fréquence est trop élevée. Par contre, dans la figure 
4, la fréquence est voisine de la fréquence propre, l'amortissement 
est faible, il Y a surtension aux bornes du condensateur. 


10 CLS:MODE 1 

20 LOCATE 14,8:PRINT'"resonance" 

30 LOCATE 1,22:PRINT'"pour faire avancer le program 
me, appuyer sur "CHR#(245) 

40 LOCATE 59,25:PRINT CHR$(243) 

930 REM 

60 REM 

70 q=INKEY(1):1F q<>0 THEN 70 

80 CLS 

82 MOVE 14,94: DRAW 14,34: DRAW 494,-26: DRAW 550,60 


84 WINDOW 1,6,1,5 

86 LOCATE 2,1:FRINT'"'u'":LOCATE 2,2:FRINT CHR$# (209): 
LOCATE 1i,3:PRINT'"'r"CHRS# (47) CHR& (208) "#f" 
90 rr=3.75 

95 FOR n=0 TO 19 

100 fr=0:GOSUB 1000 

110 MOVE x%x,y 

120 FOR fr=0 TO 2 STEP 0.05 

130 GOSUB 1000:DRAW *%x,y 

140 NEXT fr 

150 rr=rr/1.22 

160 NEXT n 

170 REM 

180 rr=53.75:fr=1i:GOSUE 1000 

190 FOR n=0 TO 59 
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200 MOVE x»%,y 

210 GOSUB 1000:DRAW x%x,y 

220 rr=rr/1i.07 

2350 NEXT n 

240 REM 

250 WINDOW 25,40,1,11 

260 LOCATE 25,1: INPFUT'"L="31 

270 LOCATE 25,2: INPUT"C="3c 

280 c=c#10"-6:fo=i/ (2#PIXSQR(1%#c)) 
290 ro=2*SGR(1/c) 

300 LOCATE 25,4:PRINT'"Fo="fo 

310 LOCATE 25,5: INPUT"f="sf 

320 LOCATE 25,7: PRINT'"KRo="ro 

330 LOCATE 25,8: INFUT'"'r="sr 

340 fr=f/fosrr=r/ro 

350 GOSUE 1000 

360 LOCATE 25,10: PRINT'u="u 

370 LOCATE 38,12:PRINT CHR$(243) 
380 TAG: MOVE x-4,y+4:PRINT CHR#& (203); 
390 Qq=INKEY(1):1F q<>0 THEN 390 
400 MOVE x-4,y+4: PRINT" “3: TAGOFF 
410 FOR a=25 TO 40 

420 FOR b=i TO ii 

430 LOCATE a,b:PRINT" “ 

440 NEXT b:NEXT a 

450 GOTO 240 

470 REM 

990 REM 

1000 REM 

1010 u=1/SGR(4#fr#frérraenmmr+(i-frafr) 2) 
1020 x=-20+(90-15#rr)+240%fr 

1030 y=((60xu) +(—-30%fr)+(90—-15#rr) ) 
1040 RETURN 






u L=? 1 
a 1 £ C=?2 06.1 
Fos S923-#9e1e1 
£f=? 503.3 
Roz= 6324.6555632 
r=? 6324 
se 6.500036077 
RER 
Fig. 2 
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Fig. 3 
u L=? 1 
cg C=? 06.1 

Fo= 3.29 
ÉF23 se3 2121 
Ro= 6324.555 
r=? 600 8 
ÿ= 3.80679979 


5. Filtrage RC 
Simple ou double alternance 


A partir d’une tension sinusoïdale Ve issue par exemple du secon- 
daire d’un transformateur dont le primaire est alimenté par le secteur, 
on désire obtenir en sortie une tension Vs unidirectionnelle, on dit ‘‘con- 
tinue”, si possible avec une faible ondulation. 
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Ve Vs 





La loi des nœuds en A donne 


i= in + à Less 
R dt 
1er cas : Si Ve > Vs, la diode conduit. Alors Vs = Ve 
2e cas : Si Ve < Vs, la diode est bloquée. i = o 
Vs , eds _5 dVs = - —!_ Vsdt 
R dt RC 


Utilisation du programme 


Entrer les données R,C,FQ. Tracer les courbes Ve, Vs. Mettre ainsi 
en évidence le rôle des entrées sur la forme de Vs. 
Représentation graphique : on veut représenter deux périodes de 
Ve. ?T doit correspondre au 640® point. Si l’on donne à t des accrois- 
sements dt, chaque accroissement fera passer d’un point à un autre. 
Donc dt.640 = 2T d'où dt = 2e ou encore dt = ee puis- 
{ 640 640.FQ 
que FQ = ca Ce calcul est fait dans les premières instructions de 


la ligne 110. 


Déroulement du programme 


lignes 10 à 70 entrées des données. A l'instant t = 0 où l'on 
met sous tension le redresseur, la tension 
d'entrée Ve n'est pas nécessairement nulle. 
Autrement dit, Ve = Vsin (Wt + PH) avec PH: 
phase à l'origine. 
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ligne 80 


ligne 100 


ligne 140 
ligne 160 


ligne 170 


Au moment où la diode cesse de conduire 
(instant t,, point A) le condensateur se décharge 
dans R suivant une loi exponentielle et Vs est 
représentée par la portion AB. Au point À, la 
pente de la sinusoïde et la pente de l’exponen- 
tielle sont égales. Au contraire, dans la partie AB, 
la tangente à la sinusoïde a une pente plus 
grande que celle de l’exponentielle. On se ser- 
vira de cette remarque dans le programme (lignes 
200 à 230). 

3 tableaux de 640 valeurs À, B, C 


À Ve, Vs 


B 


Ÿ 


Ÿ-------- 








b 


on forme 640 valeurs de t puisque t = 2.dt.F et 
que F varie de 1 à 640. On calcule 640 valeurs 
de 180.sin (2 FQt.PH) que l’on place dans le 
tableau A. On pourra ainsi les extraire du tableau 
pour dessiner la sinusoïde sur deux périodes. 
initialisation de Vs, dVs, i à zéro. 

On calcule l'accroissement D = A (F+1) —- A 
(F) qui représente dVe. 


dVs = Vs 2 . Le condensateur se décharge 


dans R. La diminution de tension qui en résulte 

est telle que : 

j = C (= dVs) ou Vs = =: SG dvs 
dt R dt 


d'où dVS = 212 Vsdt 
RC 
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ligne 180 on suppose que la diode conduit d’où la condi- 
tion i > 0 et que l'augmentation de Ve est supé- 
rieure à la diminution de Vs. Alors Vs = A (F) 
= Ve. C'est la partie OA qui sera décrite, con- 
fondue avec Ve. 


ligne 190 à partir de t,, la tension aux bornes de C baisse 
et Vs (nouvelle valeur) = Vs (ancienne valeur) 
+ dVs. 

ligne 200 on compare Vs nouvelle valeur à Ve 
SiVe > Vs ; vs = A (F) 

lignes 210 on veut tracer i (t) dans la diode d’où le calcul de 
. _ Vs cdVs 
\= — + — 

R dt 
ligne 230 si i > O0 faire i = O 
ligne 240 tableau de 640 valeurs de Vs = B (F) 


tableau de 640 valeurs de i = C (F) 
Le reste est le tracé proprement dit des courbes. 
On pourra essayer FQ = 1000 R = 500 


T , 
ee C = 10E -6, PH = 3 par exemple, ou d’autres 


valeurs. En particulier, faire varier uniquement 
C (ou R) en croissant. La courbe tend vers l’ho- 
rizontale. 


Remarques : 
Dans le cas du redressement monalternance, le lecteur pourra 
rentrer ligne 120 : À (u) = 180sin (27. FQt + PH). 


Dans le cas d’un redressemert bialternance, réécrire la ligne 120 
avec : À (u) = ABS (180sin (27. FQt + PH) ). 


Le lecteur pourra aussi faire apparaître les temps t, et t, (inter- 
valle de temps correspondant à la décharge du condensateur dans 
la résistance R). 


10 REM FILTRAGE RC 

20 REM #36 EEE 

30 REM 

35 CLS 

40 INPUT "VALEUR DE LA FREQUENCE FQG=";FQ 
50 INPUT "VALEUR DE LA RESISTANCE R="3;:R 
60 INPUT "VALEUR DE LA CAPACITE C=";c 

70 INPUT "VALEUR DU DEPHASAGE PH="3;PH 


80 DIM A(640):DIM B(640):DIM C(640) 
90 LOCATE 12,12:FRINT "CALCUL EN COURS" 


100 
110 
120 
1:30 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
330 
360 
370 
380 
390 
400 
410 
420 
430 
449 
450 
460 


FOR F=1 TO 640 

DT=1/ (FO*#640) : T=2#DT%F 
A(F)=ABS(180XSIN (2XxPI#XT4FQ+FH) ) 
NEXT F 

VS=0: DVS=0: 1=0 

FOR F=i TO 640 
D=A(F+1)-A(F):VE=A(F) 
DVS=-VS#DT / (R#C) 

IF ABS(D)>=DVS AND I1>0 THEN VS=A(F):GOTO 200 
VS=VS+DVS 

IF VE>VS THEN VS=A(F) 

I= (VS/R)+(D#C/DT) 

IF I>M THEN M=I 

IF 1<O OR VS>VE THEN 1=0 
B(F)=VS:C(F)=1 

IF F=639 THEN GOTO 270 

NEXT F 

CLS 

MODE 2 

FLOT 0,200:DRAWR 640,0 

FLOT 0,0:DRAWR 0,400 

LOCATE 1,1:PRINT "VE, VS, 1" 
LOCATE 30,1:PRINT "FILTRAGE RC“ 
LOCATE 80,12:PRINT “T" 

FOR F=i TO 640 

VE=A(F):PLOT F,A(F) +200 

NEXT F 

FOR F=1 TO 640 

VS=B(F):FLOT F,VS+200 

NEXT F 

FOR F=1 TO 640 

I=C(F):PLOT F,1/M#180 

NEXT F 

LOCATE 60,20:FPRINT “FG=";FQ; "Hz" 
LOCATE 60,21:PRINT "R="35R3"OHM" 
LOCATE 60,22:FRINT "C=";C; "FAR. * 
LOCATE 60,23:PRINT “PH="3; PH; "RAD." 
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VE, VS, 1 FILTRACE RC 








\  R- 
\ È \C= 0.008801 FAR. * 
A: “ 1%= e . 4 
; à BN= 1.57 RAD. \ 
ns Li “ “ 
, , , 
a ————_——— Le 





3 
Optique 


1. Lentilles minces 


Les instruments d'optique classiques permettent d'obtenir des ima- 
ges par réfraction. Ils sont constitués de lentilles. Une lentille est un 
milieu tel que le verre, transparent à la lumière, limité par deux calot- 
tes sphériques. On définit son axe optique. Nous nous limiterons à 
l'étude des lentilles minces. 


Programme 1 


Il permet de générer l’image d’un objet, et de visualiser la cons- 
truction géométrique qui y conduit. Il faut fournir au programme, par 
référence au bord gauche de l'écran 


- la position horizontale de l’objet, de hauteur 70 pixels, appelée 
P; 

- la position horizontale de la lentille, de hauteur 140 pixels (x 2), 
appelée F ; 

- la distance focale de la lentille, appelée DF. 


Toutes ces grandeurs seront exprimées en pixels et pourront donc 
être imaginées en mètres, centimètres, etc. 

Rappelons qu'il suffit de tracer deux rayons particuliers pour déter- 
miner la position de l’image A’ B'): 

- l’un, partant du point B (voir dessin) et parallèle à l’axe optique, 
passe après traversée de la lentille par le foyer F’ ; 


- l’autre, partant de B également, traversant la lentille en son cen- 
tre, n’est pas dévié. 
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Il a fallu tester la position de l’objet en ligne 230. En effet, suivant 
la valeur de celle-ci, l’image peut être devant ou derrière l’objet. Un 
‘drapeau’ | vaut donc + 1 ou — 1 suivant le cas, ce qui détermine 
la direction dans laquelle on trace le rayon. Celui-ci est construit à l’aide 
de la tangente de l’angle qu’elle fait avec l'horizontale, soit 70/DF à 
droite de la lentille, 70 (F —-P) à gauche. 


Plutôt que de donner une longue explication du programme, nous 
préférons donner un dessin sur lequel est porté, pour chaque segment 
de droite, le numéro de la ligne de programme qui le construit. 


Programme 2 


Ilréalise la même construction dans le cas de lentilles divergentes. 


Programme 3 


Il ne génère pas de construction géométrique. Il se contente 
d'exploiter les relations mathématiques classiques sur les lentilles. On 
peut ainsi déterminer, connaissant la position et la grandeur de l’objet, 
la position et la grandeur de l’image. 

Les exemples proposés traitent les différents cas les plus 
classiques. 


Programme 1 : Lentille convergente 
1et2 objet à plus de deux fois la distance focale F 
3 objet à une distance comprise entre 2 F et F de la lentille 


4 objet à la distance F de la lentille 
5 et 6 objet à une distance inférieure à F de la lentille 


Programme 2 : lentille divergente 


7 objet à une distance supérieure de F de la lentille 
8 objet à une distance inférieure à F de la lentille 
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Programme 3 : cas général 


9 et 10 identique à 1 et2 


11 identique à 3 
12 identique à 5 
13 identique à 6 
14 identique à 7 
15 identique à 8 


Remarques sur la construction des rayons : 
1° Un rayon parallèle à l’axe émerge de la lentille en direction 
du foyer. La tangente de l’angle qu'il fait avec l'horizontale est égale à 
longueur de l’objet 
distance focale 


soit 70/DF d’après nos notations. Un déplacement horizontal de 1 pixel 
est donc lié à un déplacement vertical de 70/DF pixels (ligne 240, pro- 
gramme 1). 


2° De même, pour un rayon passant par le centre optique, la tan- 
gente de l’angle qu'il fait avec l’horizontale est égale à 
longueur de l'objet _ 70 
distance objet/lentille F-P 
d'où la construction (ligne 240). 


3° L'intersection des deux rayons est trouvée lorsque le second 
rayon tracé rencontre un point de couleur (ligne 330). L'image est alors 
tracée. 
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R : résistance (Q) R — F F : coefficient de frottement visqueux (N/(m/s)) 
L : inductance (H) L — M M : masse (kg) 
C : capacité (F) C — «as K : raideur du ressort (N/m) 
K 
u., u: tensions (v) u. — x, x., x : déplacements (m) 
u — x 


Correspondance entre système électrique et système mécanique. 


120 à 170 





Définitions des grandeurs utilisées dans le programme : 
- en gras : grandeurs géométriques, 
- en maigre : numéros des lignes. 


10 REM 36e EEE 

20 REM * # 

50 REM * lentilles * 

40 REM + * 

SO REM HD EENE 

60 KEM 

79 INFUT'position de la lentille":f 
80 INFUT'"distance focale":df 

90 INFUT'"position de l'objet":;p 
100 MODE 2:CLS 

110 FLOT O0,180:DKAW 500,180 

120 FLOT f,40:DRAW f-6,46 
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130 PLOT f,180:DRAW f,40 

140 FLOT #f,40 :DRKAW f+6,46 

150 FLOT f,180 :DKAW f,:20 

160 FLOT f,320:DRAW f-6,314 

170 FLOT f,:20 :DRAW f+6,:14 

180 x=f-df:y=180:GOSUER 380:G6G0SUE 420 

190 x=f+df:y=180 : GOSUEH 380 

200 FLOT p,180 :DKAW p,180+70 

210 FLOT p,180+70 :DRKAW f,180+70 

220 PLOT f,180+70:DRAW f+df,180 

230 IF p * f-df THEN i=-1 ELSE i=1 

240 a=ati:b=b+70%i/df 

230 IF a+df+f>:600 OR 180-b:35350 THEN 290 

260 IF atdf+f <0 OR 180-b <O THEN 290 

270 FLOT atdf+f,1890-b : GOTO 240 

280 REM 

290 a=0:b=0 

300 1F i:0 THEN FLOT p,i8980+70 ELSE FLOT f,180 

310 a=a+i/2:b=b+70 *i/2/(f-p) 

320 IF a+p£O OR 189 +70-b+40 OK a+p:600 OK 180+70-b 
+350 THEN 350 

5350 IF TEST(a+p,180+70-b)=1 THEN FLOT a+p,180+ 70-— 

b:DRAW a+p,180 

340 GOTO 310 

5950 IF i:O THEN FLOT p,180+70: DRAW a+p,180+70-b 


360 IF 1£0 THEN FLOT f,180:DRAW a+p,180+70-b 
370 STOF 

3680 REM croix 

390 FLOT x-4,y-4: DRAW x+4,y+4 

400 FLOT x+4,y-4: DRAW x-4,y+4 

4109 RETURN 

420 LOCATE (f/8)-2,20:FRINT + 

450 LOCATE (p/8)-2,20:FRINT p 

440 LOCATE 1,22:FRINT'distance focale=":df 
450 RETURN 






distance focale= 48 
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, \ 
. L : D 


distance focale= 48 L \ 





distance focale= 48 L 





distance focale= 48 
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10 
20 
50 
40 
50 
60 
70 
80 
90 
100 
110 
120 
1:0 
140 
150 


160 





distance focale= 78 L 





distance focale= 79 L 


REM HERO ME DEEE NE 
REM + * 
KEM % lentilles * 
REM * divergentes* 
REM *% # 
REM HER EEE EX 
REM 
INFUT'position de la lentille";f 
INFUT'"'distance focale":df 
INFUT'position de 1'‘objet";p 
MODE 2:CLS 
FLOT 0,180: DRAW 500,180 
FLOT #f,40:DRAW f-6,:4 
FLOT f,180:DKAW #f,40 
FLOT #f,40 :DRAW f+6,74 
FLOT f,180 :DKAW #f,:20 


121 


170 
180 
190 
200 
210 
220 
230 
240 
2950 
260 
270 
280 
290 
300 
310 


FLOT f,220: DRAW f-6,2:26 

FLOT f$,320 :DRAW f+6,:326 
x=f-df:y=180:GOSUR 3560:GOSUR 400 
x=f+df:y=180 : GOSUR 5460 

FLOT p,180 : DRAW p,180+70 

PLOT p,180+70 :DRAW f,180+70 
a=at+i/2:b=b-70/d+f/2 

IF a-df+f:600 OR 180-b::50 THEN 280 
IF a-df+f <0 OR 180-b 0 THEN 280 
FLOT a-df+f,180-b : GOTO 220 

REM 

a=i:b=3535/(f-p):FLOT p,180+70 
a=at+1/2:b=b+70 /2/(f-p) 

IF a+p?600 OK 180+70-b 350 THEN 3350 
IF TEST(a+t+p,180+70-b)=1 THEN FLOT a+p,180+ 70- 


b:DRAW a+p,180 


320 


330 
340 
35350 
360 
370 
380 
390 
400 
410 
420 
430 
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GOTO 290 

FLOT p,180+70: DRAW a+p,180+70-b 
STOP 

STOF 

REM croix 

FLOT x-4,y-4:DRAW x+4,y+4 

PLOT x+4,y-—-4:DRAW x—-4,y+4 
RETURN 

LOCATE (f/8)-2,20:PFRINT + 
LOCATE (p/8)-2,20:FRINT p 
LOCATE 1,22:FRINT'distance focale=";:df 
RETURN 





distance focalez 7% | D 


10 
20 
30 
40 
50 
60 
70 
80 
90 
100 
110 
120 
(2) 
130 
140 
60 
150 
—{: 
160 
170 
180 
190 


220 





distance focale= 7% Il LS 


REM EEE 


REM * * 
REM # lentilles * 
REM + * 
REM HE EEE NE 
REM 
CLS 
INPUT “distance focale :"3;f 
FRINT distance de l'objet a la lentille :" 
INPUT " ( — a gauche , + a droite )":;p 
INPUT " grandeur de l'objet :"3:0 
PRINT " lentille convergente (1) , divergente 
“ 
INPUT a$ 
IF at ="i1" THEN PRINT'lentille convergente ": 
TO 170 
IF a$ ="2" THEN FRINT'lentille divergente":f= 
GOTO 170 
GOTO 140 
pi=i/(1/f+1/p) 
i=pi*o/p 


PRINT'"'la grandeur de l'image est :"::FRINT USI 
NG "#4. HH" si 

200 FRINT'"distance de l'image a la lentille :"; 
210 FRINT USING "##.##":pli 


STOF 
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distance focale : 40 

distance de l’objet a la lentille : distance de 
l'otjet a la lentille :-100 

grandeur de l’objet *: 70 

lentille convergente (1) : divergente (£) 
lentille convergente 
la grandeur de l'image est :-44.67 
distance de l’image a la lentille 146.67 


distance focale * 40 

distance de l'objet a la lentille :-:00 
grandeur de l'objet : 70 

lentille convergente (1) : divergente (2) 
lentille convergente 

la grandeur de l°image est :-17.:0 
distance de l'image a la lentille 550.00 


distance focale : 40 

distance de l'objet la lentille :-460 
grandeur de l'objet 7Q 

lentille convergente (1) : divergente (=) 
lentille convergente 

la grandeur de l'image est :%-140.00 
distance de l'image a la lentille :%120.00 


= à 


distance focale : 70 
distance de l'objet a la lentille :-40 
grandeur de l’objet * 70 
lentille convergente (1) ;: divergente (x) 

lentille “onvergente 

la grandeur de l'image est 1143. 

distance de l”'image a la lentill 


| 


33 
ë 


CR ECS 


4 


distance focale : 70 

distance de l’objet a la lentille *-33 
grandeur de l’objet : 70 

lentille convergente {1) : divergente (ZX) 
lentille convergente 

la grandeur de 1°image est :140.00 
distance de l'image a la lentille :%-70.00 
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distance focale : 70 

distance de l’objet a la lentille :-100 
grandeur de l'objet =: 70 

lentille convergente (1) : divergente (2%) 
lentille divergente 

la grandeur de l'image est : 25,87 
distance de l'image a la lentille :%-41.13 


distance focale : 70 

distance de l’objet a la lentille :-50 
grandeur de l'objet : 70 

lentille convergente (1) : divergente (2%) 
lentille divergente 
la grandeur de l’image est : 40.83 
distance de l'image à la lentille :%-29.17 


2. Lentille 


Quel est le chemin parcouru par la lumière quand elle traverse 
une lentille ? Lentille d'objectif photographique ou d’un lecteur de dis- 
que laser ou simple loupe... Quelle image donne la lentille d’un objet 
ponctuel ? Ce programme vous montre la marche des rayons lumi- 
neux et vous montre dans quelles conditions on peut espérer une 
image. C’est l'utilisateur qui détermine les caractéristiques de la len- 
tille et la position de l’objet. 


Pourquoi ce programme ? 


L'ordinateur permet ici d'aborder un domaine de la physique de 
façon entièrement nouvelle : c’est un compagnon qui vous libère des 
calculs mais vous laisse la responsabilité de la décision finale. 

Les cours de physique traitent en détail les lentilles minces qui 


sont peu utilisées. L'étude des lentilles épaisses pose en effet deux 
problèmes : 


— la seule loi optique utilisée est la loi de réfraction de Descar- 
tes, mais les problèmes géométriques conduisent à des calculs 
complexes, 
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— l'image que la lentille fournit de l’objet n’est pas nette, elle est 
mal définie. Aussi, il est impossible de la localiser avec précision, donc 
impossible de donner une formule pour calculer sa position. 

L'ordinateur vient aisément à bout de la première difficulté. Pour 
la seconde, il fournit une figure sur laquelle il vous laisser décider si 
l’image obtenue est suffisamment nette pour vos besoins et appré- 
cier sa position. 


Utilisation du programme 


La lentille est placée au centre de l’écran, dans un milieu d'indice 
1, le diamètre de la lentille doit être égal ou supérieur à la hauteur 
de l’écran. L'objet peut se situer en dehors de l’écran mais il doit être 
réel. Si l’objet est à l’intérieur de l'écran, il est possible qu’il faille donner 
au verre un indice très élevé pour obtenir une image dans les limites 
de l'écran. 

Les rayons réfléchis par les faces de la lentille ne sont pas tra- 
cés. Dans la réalité ils peuvent donner naissance à des images 
parasites. 

Dans le cas de rayons extrèmement inclinés sur l'axe, l’ordina- 
teur peut fournir des images fausses car il n’envisage pas deux ren- 
contres successives du même rayon avec la même face de la lentille. 


Pour faire avancer le programme, taper — 
Pour sortir du programme, faire ESC ESC. 


Le cœur du programme 


Rappelons tout d’abord la loi de Descartes (fig. 1) : 






milieu 2 
indice n; 








milieu 1 
indice n, 


\ 


Loi de Descartes 
n,. Sini, = M. Sin k 
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La figure 2 montre l'entrée dans la lentille et les notations utili- 
sées dans le programme. 


On a utilisé les relations 
| = E + C (ligne 505) 
Ci =R-E (ligne 520) 


avec tanC = c,tanE = te, tan! = ti, tan R = tr, sin R = sr 


ES 
_ 
_ 
_ 
—— 
_ 





centre du 
diopte d'entrée 







axe du système 


Fig. 2. - Le rayon entre dans la lentille 


La figure 3 montre la sortie de la lentille, on a utilisé les relations 
| = S — C1 (ligne 610) et C2 = S — R (ligne 630). 





Fig. 3. — Le rayon sort de la lentille 
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Liste des variables 


“6 3 


s 
V 

e 

de,ds 

re,rs 

X;Yÿ 

XO,yO 

a 

pe, ps 
xe,ye,XS,ys 
c 

ze di 

m,n 

te, ti,i,sr,tr 
ci 

ZS 

mm,nn 

c2 


en deuxième lettre désigne l'entrée de la lentille 
en deuxième lettre désigne la sortie de la lentille 
indice du verre 

épaisseur de la lentille 

type de dioptre 

rayon du dioptre 

position du point sur le rayon lumineux 
position de l’objet 

angle variable pour le dessin de la lentille 
position des faces de la lentille 

position du point pour tracer la lentille 

pente du rayon avant de pénétrer dans la lentille 
distance du point au centre du dioptre d'entrée 
position du point où le rayon entre dans la lentille 
voir figures 2 et 3 

pente du rayon pendant la traversée de la lentille 
distance du point au centre du dioptre de sortie 
position du point où le rayon sort de la lentille 
pente du rayon après la sortie de la lentille 


Déroulement du programme 
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10 

100 
250 
330 
420 
430 
450 
470 
490 
500 
930 
550 
570 
590 
600 
640 
650 


présentation 

introduction des données 

dessin du premier dioptre 

dessin du second dioptre 

départ du rayon 

avant d'arriver à la lentille 

est-on arrivé à la lentille ? 

x avance d’une unité 

entrée dans la lentille 

calcul de ci 

propagation à l’intérieur de la lentille 
est-on arrivé à la sortie de la lentille ? 
x avance d’une unité 

sortie de la lentille 

calcul de c2 

après la lentille 

arrivée en bout d'écran 


660 x avance d’une unité 


700 choix d’un nouveau rayon, retour à 430 
710 le dernier rayon vient d’être tracé 
800 conclusion 


Commentaire des illustrations 


Fig. 4 : l’objet et l’image sont en dehors de l'écran. 

Fig. 5 : obtenue avec dans l’ordre : 3, 200, 1, 400, 1, 400, 0, ©. 
L'indice du verre étant plus élevé que précédemment, on 
peut observer l’image et l’objet sur l'écran car la lentille 
est plus convergente. L'image n’est pas ponctuelle car 
les rayons à la sortie de la lentille ne se recoupent pas 
au même point. 

Fig. 6 : même lentille, mais l’objet est en 0, 100. 

Fig. 7 : avec 2, 100, 1, 600, — 1, 240, O, 100 
La lentille est divergente. 

Fig. 8 : avec 1,5, 60, —-1, 700, — 1, 700, 150, 100. 

Dans les deux derniers cas, vous pouvez prolonger vers 
la gauche les rayons sortant de la lentille pour rechercher 
l’image virtuelle de l’objet. 


10 CLS 

20 LOCATE 14,8:PRINT'"Lentilles" 

30 LOCATE 14,10: FRINT'"epaisses" 

40 LOCATE 1,22:FRINT'"pour faire avancer le program 
me, appuyer sur'"CHR#(245) 

930 LOCATE 35,25:FRINT CHR#(243) 

60 REM 

70 REM 

80 q=INÉEY(1):1F q<>0 THEN 80 

100 FRINT:INPFUT'indice du verre="3;v 
110 FRINT: INFUT'"epaisseur de la lentille'";e 
120 FRINT:FRINT'"premier dioptre" 

130 INPFUT'"'concave:-1,convexe:l';de 
140 INPUT'rayon (3200) ='"'sre 

150 FRINT:FRINT'"deuxieme dioptre" 
160 INFUT'"concave:-1,convexe:l";ds 
170 INFUT'"'rayon (3200) ="3rs 

189 PRINT:FRINT'"position de l'objet" 
190 INPUT'"'x="3x»x0o: INFUT'"'y=": yo 

200 LOCATE 539,25:PFRINT CHR#(243) 

210 q=INKEY(1):1F qg£>O THEN 210 
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Z20 
240 
250 
269 
270 
280 
290 
300 
310 
320 
325 
33O 
540 
350 
360 
370 
380 
385 
390 
400 


410 
420 
420 
440 


450 
4860 
470 
480 
490 
500 
905 
510 
520 
920 
940 
990 
960 
970 
580 
590 
600 
610 
620 
625 
630 
640 
68350 
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REM 

CLS 

MOVE 0O,200:DRAW 640,200 
pe=320-e/2:ps=320+te/2 

a=0 

xe=petde*re+*(i-COS(a)) 
ye=re*SIN(a) 

IF ye>200 THEN GOTO 330 

PLOT xe,200+ÿe:FLOT xe,200-ye 
a=a+0.005:GO0OT0O 280 

RENM 

a=0 

*xs=ps-ds#rs*(i-C0S(a)) 
YS=rsS#SIN(a) 

IF ys>200 THEN GOTO 390 

PLOT xs,200+ys:PLOT xs,200-ys 
a=a+0.005:GO0OTO 340 

REM 

IF xs<xe THEN CLS:GOTO 100 
MOVE »xe,599:DRAW xs,:599 : MOVE »xe,0:DRAW xs,0 


REM 

x=xo:c=(160-y0o) / (pe-xo) 

Y=yYo+cCXx(x—-x0) 

1F ABS(y)<200 THEN FLOT x“,200+y ELSE GOTO 700 


ze=de*((x-pe-dexre)"2+y"2) 

IF ze“dex(re”2) THEN GOTO 490 
x=x+1:GO0TO 430 

REM 

ME=X 3: N=Y 

te=n/(petde*re-m) 

ti=(te+c)/(i-te*#c) :i=ATN(ti) 
sr=SIN(i)/vstr=sr/S0OKR(i1-sr"2 
ci=(tr-te)/(i+trrte) 

y=n+cix(x-m) 

IF ABS(y)<200 THEN PLOT x,200+y ELSE GOTO 700 
zs=ds*#((x-ps+ds+rs)*2+y7"2) 

IF zs>ds#(rs*"2) THEN GOTO 590 
x=x+1:GO0OTO 530 

REM 

mMm=}x:NN=Yÿ 

ts=nn/(mm-ps+ds*rs) 
ti=(ts-ci)/(i1+ts*xcil):i=ATN(ti) 
sr=SIN(i)#vV:1F ABS(sr)2>1 THEN GOTO 700 
tr=sr/SOR(i-sr"2) 

c2=(ts-tr)/(irttræts) 

ÿY=nn+c2#%(x-mm) 

IF ABS(y)<200 THEN FLOT x,200+y ELSE GOTO 700 


660 x=x+i1:IF x<640 THEN GOTOQ 640 

700 c=c-0.2 

710 IF c<(-160-ÿo) / (pe-xo) THEN GOTO 800 
720 x=xo:GOTO 430 

790 REM 

800 TAG:MOVE 620,20: PRINT CHR# (243) ; : TAGOFF 
810 q=INKEY(1):1F qg<:0 THEN 810 

820 CLS 

830 FRINT:PFRINT'meme lentille? "CHK#(243)",sinon E 
SC ESC" 

840 qg=INKEY(1):1IF qi>O THEN 840 

850 GOTO 180 

860 END 


Ar et avancer le programme, appuyer 
sur 


+ 

indice du verre=?7 1.5 
epaisseur de la lentille? 2060 
premier diopire 
concave:-i,convexe:1? 1 
rayon (2206) =? 406 
deuxieme dioptre 19 1 

e : on 
22004" 252065°23°486 
position de l’objet 


x=?2 
u=? © + 


Fig. 4 
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3. Dioptre plan 


Ce programme permet de visualiser le trajet d’un rayon lumineux 
lorsque celui-ci rencontre un dioptre plan. Nous ne donnerons pas la 
démonstration mathématique complexe qui conduit à l'algorithme 
employé. Seul compte le résultat ! 


Ce programme est largement modifiable et peut donc être adapté 
à la mise en évidence de phénomènes particuliers. 


La position du dioptre plan est imposée ligne 110 (X= 200) ; la 
position du point lumineux intervient ligne 150 (X = 100) et ligne 160 
(où le nombre 100 représente la distance entre ce point et le dioptre). 
Test l’angle de ‘‘départ”” du rayon, en radians. Les valeurs prises par 
T sont imposées ligne 140 et 290. 


10 REM HR fee A ME AE EE EE EE 


£Q KREM # # 
30 KEM *# dioptre plan # 
49 REM # # 
SO REM tt tete eee 
&Q REM 


70 MODE Æ# 5 CL: 
&Q INPUT 'nis"snd 
90 INPUT Mnésetsns 
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100 CLS: a0=ni/n2 

110 FLOT Z00:0:LRAW 700,377 

120 LOCATE 1,1: FRINT'ni="sni 

130 LOCATE 4%: 1: FRINT'"'nz="snZz 

140 FOR t=-FI/4 TO(FI/4+FI/32)STEP F1/14 
170 PLOT 100: 200%s=-Û0ÛE (t) 

1640 DRAW 00: ZO00+100#TAN(t) 

170 cenETE -nitZ +(ni+#s) "2 

150 IF s<O THEN a=-1:b=2#C0S(t):GOTO 00 
170 a=a0 th=s+(a0-(c/ (n2+s) T2) 70.5) 

“OO deSIN(t)/(b-akrs) 
X=2003y=2O00+100#TAN(Ct) 
x=Xx+5SON (CC): y=y+d 

IF x 5657 OR x<O OR y2377 OR y£0 THEN 250 
240 DRAW xy:G0TO 270 

250 IF céO THEN 270 

LEO x=x-1:ÿ=y-d 

270 IF x 3637 OR x<O OR y>397 OR y < © THEN 2970 
280 DRAW x:y : GÜTA Z60 

270 NEXT t 

300 STOF 





ni= 1 nez 1.5 
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4 
Chimie 


1. Propriétés d’une solution chimique 


Ce programme propose de calculer le pH d’une solution aqueuse, 
connaissant le pKa du couple acide-base correspondant, ainsi que sa 
molarité. 


Un couple acide-base, représenté par la réaction : 
A + H,0 a B + H,0* 
est défini par son Ka : 


ka - (B) (H:0*) 
(A) 
on a aussi pKa = -— log Ka. 


Pour une solution d'acide ou de base, on peut écrire (pour plus 
de commodité, on pourra prendre l’exemple de l'acide acétique) : 

CH;CO;H + H,0 = CH;,CO, + H,0 + 

On a alors (A) = (CH;CO.H) et (B) = (CH;,CO,°) 

(B) + (OH°) = (H30*) 

(A) + (B) = C (molarité initiale) 

(H3O*).(OH=) = 1074 


ka - (B) (H:0:) 
(A) 
On en déduit donc : 
10:77 
OH°) = — B) = (H,30*) —- (OH- 
pes (B) = (H:0*) - (OH-) 
107 
B) = (H:0*) - 
ne 
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(A) = C - (8) = C - (H,O*) - 107 




















; (H30*) 
(H0*) - es . (H50*) 
d'où Ka = ss 
C2 oO 
(H3O*) 
Cette équation relie le pH = — log (H:0*) ) la molarité de la 


solution. Elle dépend de Ka. 
il est donc possible de l’utiliser des quatre manières suivantes : 


— Tracer des courbes universelles (graphique de Flood) donnant 
le pH en fonction de C, avec comme paramètre le pKa. 

— Calculer le pH d’une solution connaissant le pKa et la mola- 
rité C (3). 

— Calculer la molarité C, connaissant le pKa et le pH de la solu- 
tion (1). 

— Calculer le pKa du couple, connaissant le pH de la solution 
et sa concentration (2). 


Ces trois dernières possibilités seront regroupées dans le 
programme. 


Liste des variables 


P pH 

H (H30*) 

K pKa 

KA Ka 

C concentration 
PK  pKa 

PA pKa 


Détail du programme 


Il permet d'étudier une solution particulière : 


80 - 170 Premier cas:  pKa et pH connus. C est obtenue par 
application directe de la relation de 
base. 
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19 
20 
30 
40 
90 
680 
70 
80 
90 
1900 
110 
(1) 
120 
(2) 
1:50 
3) 
140 
150 
160 
170 
180 
199 
200 
210 
220 
230 
240 
2350 
260 


200 - 290 Deuxième cas : 


450 - 620 Troisième cas : 


C et pH connus. Partant de pKa = 0,1, 
toutes les valeurs de pKa sont exami- 
nées jusqu’à ce que la formule de 
base donne une valeur de C presque 
identique à celle qui est donnée. 
pKa et C connus. Partant de pH = 1/2 
pKa, toutes les valeurs de C données 
par la formule sont examinées. Le cas 
d’une base est traité comme celui d’un 
acide en utilisant son pKb. La valeur 
finale du pH est alors 14 — (pH 
trouvé). 


630 - 640 Calcul de la concentration Ke = 10-74. 
630 Formule de base, donnant la concentration C en fonc- 
tion de (H*) et de Ka. 


REM HE END DEE DE DEEE 
REM *# 

REM *% SOLUTION CHIMIQ 
REM +% 

REM MH DE DE DM DEEE DEEE 
REM 

KHE=10"(-14) 


CLS 
FRINT 
FRINT " faites votre 
FRINT " connaissant 
PRINT " connaissant 
FRINT "“ connaissant 
FRINT 
PRINT'tapez 1,2 où 3 
INPUT n 
FRINT'choix numero: 
IF n=i OR n=2 OK n=5 
ON n GOSUB 200,300,4 
FRINT:PRINT:PRINT 
INPUT "“ph=":pth=10"( 
FRINT:PFRINT:FRINT 
INFUT'"pka=";k:ka=10" 
IF p>7 THEN h=ke/h:k 
GOSUB 6350:FRINT:FRIN 
FRINT" la concentr 


XX 
+ 
UE * 
+ 
EX 


choix 
pka et pH , vous voulez c 


c et pH ,; vous voulez pKa 


pkHa et c, vous voulez pH {( 


u 


sn 
THEN 190 ELSE 150 
60 

—p) 

(—K) 

a=ke/ka 


T:c=f 
ation est c="; 
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270 
280 
290 
300 
310 
520 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
450 
440 
450 
460 
470 
480 
490 
500 
910 


520 
5:30 
540 
590 
560 
970 
980 
990 
600 
6410 
629 

T0 
640 


FRINT USING'"#4.##H#H##":Ccs:FRINT" mole/litre"“ 
ZF=INÉEYS#S:1IF z$="" THEN GOTO 280 

GOTO 10 

PRINT:PFRINT:FRINT 

INFUT "“concentration=";c 

PRINT: INFUT'"ph="3;p:h=10"(-p) 

IF p>7 THEN h=ke/h 

pk=O. 1 

pk=pk+0.1:ka=10"(-pk):GOSUER 6350 

IF fc THEN 2370 ELSE 3550 

pk=pk-0,.1 

pk=pk+0.01:ka=10"(-pk):GOSUBR 650 

IF fc THEN 400 ELSE 380 

IF p >7 THEN ka=ke/ka 

FRINT:FRINT'le pka est:"; 

FRINT INT(100%(-LOG(kHa))/2.3)/100 
ZH=INKEY#:1F z$="" THEN GOTO 430 

GOTO 10 

PRINT: FRINT:FRINT 
INFUT'"pka=";pa:ka=10"(-pa) 

PRINT: INFUT'"concentration=";c 
FRINT:FRINT'"'solution acide ou basique?" 
INFUT af 

IF a$t="a"THEN FRINT'solution acide" 

IF at="Db'"'THEN ka=ke/ka:FRINT'"'solution basique“ 


p=0.Sx*pa:IF a$="b" THEN p=0.5*(1i4-pa) 
p=p+0.1:h=10"(-p):GOSUR 630 
IF f# € THEN 5350 ELSE 5250 
p=p-0.1 
p=p+0.01:h=10"(-p):GOSUE 6:30 
IF f #4 € THEN 580 ELSE 560 
IF at="b" THEN p=1i4-p 
FRINT:FRINT 
PRINT" le pH est ";INT(p#100)/100 
ZHF=INKEY#:IF z#="" THEN GOTO 4350 
GOTO 10 
f=(h-(ke/h))#((h/ka)+1) 
RETURN 


faites votre choix 


Cormaissant pla et pH : vous voulez © * (1) 
Connaissant &« et pH : vous voulez pKHa : (2) 
Connaissant péa et © : vous voulez pH 5 (5) 


tapez 1: où 3 


choix mumeroï 
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concentrations 0.01 
phe= 11 
le pka est: 10,0% 


faites votre choix 

LConnaissant pKa et pH : vous voulez © (1) 
Connaissant c et pH :; vous voulez pa : (2) 
Connaissant pKa et © : vous voulez pH (3) 
tapez 1:# ou 3 

choix numeroï 3 

phaz= 10 


faites votre choix 

Cornnaissant pka et pH >; vous voulez & #(1) 
Connaissant &c et pH : vous voulez pKa : (2) 
Connaissant pEa et & : vous voulez pH : (3) 
tapez 1: ou 3 

choix numeroï 3 

pka= 10 


Solution acide où hasique * 
solution acide 


le pH est 4.97 


faites votre choix 

Connaissant pka et pH : vous voulez = :(1) 
Connaissant © et pH : vous voulez pa : (£) 
Cornaissant pla et © : vous voulez pH :(3) 
tapez 1:Z ou 3 

choix numeroëï 3 

pha= # 


solution acide où basique 
solution Hhasique 


le pH est ‘7.54 
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faites votre choix 

Connaissant pa et pH : vous voulez © (1) 
Lonnaissant = et pH : vous voulez pKa 5 (2) 
Lonnaissant pa et © »: vous voulez pH #:(3) 
tapez 1:2 ou 3 

choix numeroï 1 


ph= 4 


pkha= 6 
la concentration est c=0.0101 mol/litre 


faites votre “choix 

Connaissant pa et pH >; vous voulez € #(1) 
Connaissant © et pH :; vous voulez pKa : (2) 
Connaissant pKa et © : vous voulez pH # (3) 
tapez 1: ou 3 

choix numeros 1 


ph= 10 


pha= 7 
la concentration est c=0.1001 mol/litre 


faites votre choix 


Lonraissant pa et pH : vous voulez € (1) 
LConnaissant c et pH : vous voulez pKa : (2) 
Connaissant pkÉa et € : vous voulez pH (3) 


— 


tapez 1:# ou 3 
choix nmumeroï %Z 


concentration= 0.01 
ph= 4 
le pka est: € 
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2. Dosages volumétriques 


Introduction 


Le dosage volumétrique d’une solution s'effectue en ajoutant m 
à un volume déterminé de cette solution un réactif convenable. 


La fin de la réaction (équivalence) est le plus souvent marquée 
par le changement de coloration d’un indicateur coloré. 


Le titre d’une solution est donné par sa molarité (nombre de moles 
présentes dans un litre de solution) ou par sa normalité (nombre de 
valences-grammes par litre de solution). 


Acidimétrie et alcalimétrie 


— Une valence-gramme d'acide est équivalente à un ion-gramme 
H,0*. 

— Une valence-gramme de base est équivalente à un ion-gramme 
OH. 

— La normalité d’une solution acide est exprimée par le nombre 
d’ions-grammes H;0* présents dans un litre de solution. 


Les valences-grammes des principaux acides sont : 


HC1 : 7 HO: : H,PO, : CH;COOH etc. 


— La normalité d’une solution basique est exprimée par le nom- 
bre d’ions-grammes OH” présents dans un libre de solution. 


Les valences-grammes des principales bases sont : 


NaOH : KOH : NH,OH : , 8e (OH), etc. 


Soit v, et v,, N, et N, les volumes et les normalités des solutions 
d'acide et de base qui ont réagi. A l’équivalence, c'est-à-dire pour le 
même nombre de valences-grammes d'acide et de base correspon- 
dant à la réaction complète. 


H:30* + OH —- 2H,0, on peut écrire : N;v, = Nov 
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Equilibres ioniques 


Les problèmes d’équilibres ioniques reviennent fréquemment à 
des calculs de concentrations ioniques. Le plus souvent, on calcule 
la concentration en ions H,0* qui se traduit par le PH (PH = -— log 
(H30*)). 

Les difficultés du calcul résident dans l’établissement correct des 
relations dont on peut éventuellement disposer. 


Intéressons-nous au cas des équilibres ioniques homogènes (la 
solution aqueuse est l'unique phase). 


1° Equilibres présents en solution 


Une solution aqueuse (solvant) contient une grande proportion de 
molécules d’eau dont on peut considérer la concentration comme 
constante. 


Parmi les ions apportés par le ou les solutés figurent toujours des 
ions H,0* et OH” en équilibre avec les molécules d’eau. Cet équili- 
bre peut se traduire par l'équation suivante : 

2H,0 —H;0* + OH” 

De plus, nous classerons les ions présents en solution en deux 
groupes : les ions ‘“‘forts’”’ et les ions ‘‘faibles””. 


— Union est dit ‘fort’ lorsqu'il est sans action sur les autres ions 
ou molécules de la solution ; sa concentration est celle qui correspond 
à son introduction en solution. 


— Un ion est dit ‘‘faible”’ s’il déplace l'équilibre précédent en 
créant un nouvel équilibre. 


Ex. : L'ion acétate CH,COO ” agit sur les ions H,0* présents en 
solution selon l'équilibre. 


CH;COOH + H,0==CH,COO- + H,0* 


En présence de soude (NaOH) par exemple, les espèces présen- 
tes en solution sont : 


Na‘ ; H30* ; CH,COO- ; OH et CH;COOH 


2° Relations entre les concentrations ioniques et molaires 


Ces relations utilisent les données : (volumes, concentrations, etc.). 
Elles peuvent être classées en deux groupes : 
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— Relations linéaires : 
° Expression de la neutralité électrique. 


On écrit que, dans un litre de solution, la somme des charges élec- 
triques (exprimées en électrons-grammes ou faradays) portées sur les 
cations est égale à celle portée par les anions. De cette manière, ce 
sont les concentrations ioniques qui interviennent dans la relation. 


(Attention : certains ions peuvent éventuellement porter plusieurs 
charges). 


Ex. : Pour une solution de sulfate ferrique (Fe,(SO,);) : 
3(Fe°*) +(H,0*) = 2 (S0,2-) + (OH°) 
+ D'autres relations linéaires sont obtenues en exprimant que la 


solution contient ce qui a été introduit (Loi de conservation de la 
matière). 

Ex. : Pour une solution d'acide acétique de concentration 
molaireC, 

(CH;COOH) + (CH;:COO') = C, 

— Relations non linéaires : 


Ce sont les relations obtenues par application de la loi-d’action 
de masse aux équilibres, y compris évidemment l'équilibre ionique de 
l’eau. 


(H:0*) (OH) = Ke (avec K, = 107“). 
Ex. : Toujours pour l'acide acétique 
CH,COOH + H,0 = CH,COO- + H,0* 
On peut écrire : 

(CH;COO!) (H:0*) 


KE L- avec K, = 1.85.10 ° 
(CH,COOH) 


Neutralisation d’un monoacide par une base 
1° Considérons l’acide HA de concentration molaire C, en 
solution 
— Appliquons la loi d'action de masse à l’équilibre 
HA + HO=Z= A7 + H,0* 
Remarque : Pour simplifier l'écriture, nous écrirons dorénavant : 
HA==AT + H* 


145 


On peut écrire : 
K,= AH) (si l'acide est fort (HA) — 0, donc K; — «) 
(HA) 

La loi de conservation de la matière entraîne : 
(HA) + (A°) = C; 
de plus (H*) (OH-) = K,  (K, = 107‘) 
Ecrivons l’électroneutralité de la solution acide : 
(H*) = (A) + (OH°) 
On montre facilement que : 

: K H+ 
C, = (> Ke) «1 + M), 

(H*) K; 

Remarque : (H*) = 10°" 


2° Considérons la base BOH de concentration molaire C, en 


solution 
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— Appliquons la loi d'action de masse à l’équilibre 
BOHZ=B* + OH” 
On peut écrire : 
2 = (ES OPE) (si l'acide est fort (BOH — O0 donc K; — ©) 
(BOH) 
La loi de conservation de la matière entraîne 
(BOH) + (B*) = C, 
Ecrivons l’électroneutralité de la solution basique : 
(B°) + (H*) = (OH°) 
On montre facilement que : 
Ke 


Ke 
C; = (1 + — (H* 
2 na 





3° Considérons le mélange BOH + HA 
Ecrivons toujours les équations : 
K, =CA) (HT) . L - (B°) (OH) 


(HA) ? (BOH) 
(H*) (OH°) = K, et  (H*) + (B*) = (OH°) + (A°) 


Supposons que l’on soit parti d’un volume v, de solution acide 
et que l’on ajoute un volume v de solution basique. 


La loi de conservation de la matière conduit aux deux nouvelles 


équations 


(HA) + (A7) = 


Go. et (BOH) + (B*) = -CŸ. 
Vo+V Vo+v 


On montre donc facilement que : 











C; ___Ke (H*) 
PCOCE 

CG  __ Ke , gs) 

__K___ (H*) 

KH”) 


Liste des variables 


K1 
K2 
p1 
p2 
K 
VO 
R 
C1 
C2 
Vi 


V2 
B 


constante d’acidité 

constante de basicité 

PH de la solution acide de départ 

PH de la solution à la fin de la réaction 
produit ionique de l’eau 

volume de la solution acide initiale 
facteur d'échelle pris à 1 au départ 
concentration de la solution acide pur 
concentration de la solution de base pure 
volume de la solution correspondant à la neutralisation 
totale 

volume de solution basique ajoutée 

PH du mélange 


Déroulement du programme 


lignes 10 à 100 entrées des données 


Remarque : le volume VO est pris égal à 
50 cms. Le facteur d'échelle (ligne 50) R est 
pris égal à 1 au départ (l’utilisateur peut éven- 
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10 
20 
30 
40 
Te) 
99 
60 
70 
80 
90 
100 
110 
120 
1:20 
149 
150 
160 
170 
180 
190 
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lignes 110 à 150 


lignes 160 à 210 
ligne 230 


lignes 250 à 340 


ligne 270 


ligne 400 


ligne 410 


ligne 460 


tuellement le changer selon le problème 
étudié). 

on calcule les concentrations C1 et C2 de la 
solution acide et de la solution basique 
tracé des axes de coordonnées 

on fait varier le PH (B au clavier) de P1 jusqu’à 
P2 

on calcule pour chaque valeur du PH le 
volume de base V ajouté correspondant 

la courbe PH = f (V) est tracée en continu 
l'ordinateur imprime le volume V2 correspon- 
dant à la neutralisation totale (V2 = V1) 
nous trouvons un point d'arrêt ; l'ordinateur 
interroge l'utilisateur : ‘Avez-vous observé la 
neutralisation totale ? » 

entrée d’une variable chaîne A$ correspon- 
dant à “OUI” ou “NON” 

si A$ = “OUI”, l'ordinateur imprime : ‘‘Expé- 
rience terminée” 

si A$ = ‘‘NON’’, l'ordinateur imprime ‘‘Nou- 
veau facteur d'échelle ?”’ L'utilisateur peut 
rentrer un nouveau facteur d'échelle R et l’on 
recommence. 


REM COURHRE FH=F (V) 
REM HE EE EEE 


REM 


REM INTRODUCTION DES DONNEES 
REM HE DEEOEHEDEDEDEDE HE DEN DEN DEN DEN DEN 


CLS 


INFUT "CONSTANTE D'ACIDITE Ei=";Kk1 
INFUT "CONSTANTE DE HASICITE KH2=";r2 
INFUT "FH DE DEFART Fi=";Fi 

INFUT "FH D'ARRIVEE F2="; F2 


CLS 


EH=1E-14:VO=S0:R=1 
REM Ci=CONCENT.DE L'ACIDE 
Ci=(10"(-F1)-K/10%C-F1))#(1+10% (F1) /K1) 
REM CZ=CONCENT.DE LA HASE 
2=(i1+KH/ (ER#107(C-F2)))#(H/10%(-F2) -10%(-F2) ) 


Vi=CixVvo/c2z 
MODE 2% 


FLOT O,0:DRAWR 0,400 
FLOT 6,0:DRAWR 640,0 


200 LOCATE 2,3:FRINT "FH" 
210 LOCATE 78,24:FRINT "V" 
220 T=0:G=0 
2:50 FOR B=F1 TO F2 STEF 0.05 
240 H=10"*(-H) 
2950 V=VO#(C1/(1+10%(-B) /KH1)+KH/107%(-B)-10*(-B) ) / (C2 
J U+E/ (H2#10%(-HB)))-EH/10%*(-H)+107(-HB)) 
260 IF G=1i THEN GOTO 300 
270 IF Vr=Vi THEN G=1:LOCATE 50,10: FRINT "TL NTR'": 
LOCATE SO,11: PRINT "V2="3V1:LOCATE 50,12: FRINT "FH 
=":B 
300 IF V/K:=3:00 THEN GOTO 400 
310 IF T=0 THEN FLOT V/R,24*H:T=1:GOTO0 3230 
320 DRAWR V/K-AV,24#xE-AE 
350 AV=INT(V/R+0,.5) : AB=INT(24*xEH+0.5) 
340 NEXT E 
400  LOCATE 1,1:FRINT "AVEZ VOUS OBSERVE LA NEUTRA 
LISATION TOTALE?" 
410 INFUT "A$="3;A# 
20 IF A$="OQUI" THEN GOTO 500 
450 IF A+="NON" THEN INFUT "NOUVEAU FACTEUR D'ECHE 
LLE R=";K:GOTO 170 
900 LOCATE 50,18:FRINT "VOLUME D'ACIDE VO="Vo 
910 LOCATE 50,19:PRINT "Hi="s;ki 
920 LOCATE 50,20:FRINT "E2="sk2 
930 LOCATE S50,21:FRINT "FH DE DEFART="FI 
940 LOCATE S50,22:PRINT "FH ARRIVEE=":F2 
990 LOCATE S0,23:FRINT "FACTEUR D'ECHELLE="3;R 


AVEZ VOUS OBSERVE LA NEUTRALISATION TOTALE? 
nie! OuI 
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+ à à OBSERVE LA NEUTRALISATION TOTALE? 








L 
ve= 939.76311 
PH= 7. 


VOLUME D'ACIDE VE= 58 
1 
Pi DE DÉPART 2.7 

U ARRIVÉE- 11 


FACIEUR D'ECHELLE= 1 


$ OBSERVE LA NEUTRALISATION TOTALE? 


P 
FACTEUR D’ ÉCHÉLLE- 1 





AVEZ VOUS OBSERVE LA MEUTRALISATION TOTALE? 
LEA QUI 





150 


3. Dissociation d’une solution chimique 


Considérons un acide AH, susceptible de se dissocier en solution : 
AH + H,0 A7 +H,0* 


Si l’on appelle : 


H la concentration en ions H30* 
KA la constante de dissociation 
A le degré de dissociation 

AO la concentration de la solution 


On démontre les relations suivantes : 
KA 
KA+H 


(2) AO = [H- ; IA  (K = 107‘) 


()  A= 


Pour une solution donnée, si H et KA sont connus (le plus sou- 
vent par l'intermédiaire de pH et pKa), on peut grâce à (1) en déduire 
À, puis avec (2) en déduire A0. La solution étant alors parfaitement 
définie, il est possible d’en déduire la courbe A = f (pH). Le programme 
exploite donc directement ces deux relations. 


Dans un premier temps, le programme évalue AO à partir du pH 
de la solution ainsi que de son pKa. 


Après appui sur une touche, le programme trace la courbe 
A=f (pH) pour la solution donnée. 


10 FEM Re CES ES LES DES LES LES LES Lu LE SES SES CES LES LES En Eu SES SES CE LES SES LES LES SES 6 


20 KEM # # 
30 REM # degre de dissociation *# 
40 REM # d'une solution # 
0 FREM # #: 
LO RENM capote eat op ote me ef 
7 KEM 

30 CLS ? MODE 

0 E=iE-14 


109 INPUT" ph de La solutions"sph 
119 INFUT'pEa de la solutions ":pk 
1E0 KHa=10"{-ph)ih=10" (ph) 

1350 a=( Ea/(ka+th)) 

140 a0=(h-k/h) /a 
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150 
140 
170 
150 
170 
FO 
219 
Zz0 
#50 
40 
T0 
40 
#70 
E50 
z70 
300 
310 


320 


FRINT'la concentration 


PRINT LISINGUH.HH#H "ras: 


a=INEEYS:IF a$="" THEN 
FLOT 30:30: DRAUR 400,0: 
FLOT 30:30: TIRAUWR 0,310: 
LOCATE 57:2£:FRINT"7" 
FÜR v=l TO € 

FPLOT 25, 1204#n+30: DRAWK 
NEXT n 

LOCATE 4:14: FRINT"1/E" 
FÜR nzi TO 13 


PLOT S0+40%: 253 TIRAUR 0: 


NEXT n 

FÜR ph=i TO 13 STEF 0.1 
h=10%(-ph) 

a=ka/(kEa+h) 

FLOT 30+40 #ph:30+300%a 
NEXT ph 


pr de Ja solution:? 2 
Ï a de la solution:? 3 
a concentration est:0.110 nol/1. 


a A Rd 


1/8 


esti"s 


FRINT" mol/1." 


170 


LOCATE 75: 2% 


LOCATE 4:53 


4:00 


4 


: FRINT"ph" 
FRINT"a0" 
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h de la solution:°? 4 
a de la solution:? 
la concentration est: 4. 198 nol/1,. 





ph de la solution:? 6.5 
a de la solution:? 18 
a concentration est:9.891 nol/1. 





4. Couches électroniques 


Le cauchemar des étudiants débutant en chimie ! En effet, de 
manière générale, si la répartition en ‘‘couches”’ des électrons est sim- 
ple (il y a sept couches possibles, distinguées par les lettres K à Q), 


153 


l’ordre de remplissage est plutôt complexe. Ce programme propose 
de visualiser la manière dont les couches se remplissent pour un atome 
quelconque. Cet ordre est en fait imposé par la ligne des DATA (ligne 
210) qui contient des groupes de trois données étant respectivement 
le numéro de la couche (1 à 7), le type de la sous-couche (s, p, d, en 
f), et le nombre maximal d'électrons qu’elle peut contenir. 


Le nombre d'électrons (soit le numéro atomique) de l’atome est 
appelé S. On cherche d’abord si la sous-couche en cours sera com- 
plètement remplie ou non (lignes 140 - 160). Dans un cas (ligne 140), 
on retourne en 110 pour chercher à remplir d’autres sous-couches. 
Dans l’autre cas, il s’agit de la dernière, qui sera alors complète. Si 
enfin la dernière sous-couche ne doit être qu’incomplètement remplie, 
on recherche le nombre d'électrons à placer. 

L'intérêt pédagogique du programme est certain : l'étudiant ‘‘voit” 
les couches électroniques se remplir et admet ainsi qu’il y a un ordre 
à respecter. Nous présentons les résultats obtenus pour l'hélium 
(N=2), le carbone (N = 6), le fluor (N =9), l'aluminium (N = 13), le fer 
(N = 26), l'argent (N=47), le plomb (N = 82), l'uranium (N = 92). 

Un second programme, qualifié de variante, permet simplement 
de générer les couches électroniques d’un atome. 


10 REM XX EHHEHOHEOCHEHEOUNE REED EEE 


ZOO RKEM * * 
50 REM * couches electroniques * 
40 KEM * * 
SO REM #66 6 DE DE DE DEEE EDEN DE DEEE DEEE 
60 RENM 
70 u=0 


80 INFUT "nombre d'electron:'";ss 

90 CLS: LOCATE 2,2:FRINT'"'n=";s 

100 FOR z=i TO 7:x=256:y=192:r=24%#z: GOSUR 280: NEXT 
110 READ a$,b#f,c# 

120 u=u+VAL(c#) 

150 IF uis THEN GOSUR 180:GO0T0 110 

140 IF u=s THEN GOSUER 180 

150 IF uïs THEN b=VAL(c#)-(u-s)-1:c#+=STR# (b) : GOSUE 
210 

169  STOF 

170 DATA 1,5,2,2,5s,2,2,p,6,2,s,2,5,p.6,4,5,22d3i 

0,4,p:6,5,5:24,d,10,5,p,6,6,5,2,4,f,14,5,d,10,6,p 
.6.:6,d,10,7,S,2 
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180 IF c$="2" THEN c$="1" 
190 IF c$="6" THEN c$="Ss" 
200 IF c#t="10" THEN c#="9" 
210 IF a$="1" THEN a=k:b=Kk+VAL (c#):c=2: GOSUR 290:k 
=k+VAL (c$)+1:RETURN 
220 IF a$="2" THEN a=l:b=1+VAL(c#):c=8:GOSUB 290:1 
=] +VAL (c$) +1: RETURN 
220 IF a$="3" THEN a=m:b=m+VAL(c#):c=18: GOSUB 290: 
m=m+VAL (c$)+1:RETURN 
240 IF a$="4" THEN a=n:b=n+VAL(c#):c=52: GOSUB 290: 
n=n+VAL (c$)+1:RETURN 
250 IF a$f="S" THEN a=o:b=o+VAL(c#):c=22: GOSUB 290: 
o=o+VAL(c#)+1:RETURN 
260 IF a$="6" THEN a=p:b=p+VAL(c#):c=32: GOSUR 290: 
p=p+VAL(c$)+1:RETURN 
270 IF a$="7" THEN a=q:b=q+VAL(c#):c=32: GOSUR 290: 
q=q+VAL (c#)+1:RETURN 
280 FOR t=0 TO 2*FI STEP 2*F1/120:FLOT x+rxCOS(T), 
ÿY+rA#SIN(T):NEXT t: RETURN 
290 FOR z=a TO b:x=256-VAL(a$) *24#xC0S(2*xFI#z/c):y= 
192-VAL (a+) x244xSIN(2xFI#xz/c):r=6: GOSUR 280: NEXT 23: 
RETURN 

n= 2 





7 one, 
s 
Pa ,. 
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96 


32e 


ee 
set 


=u 


€T 


RP 
# \ . 
W , L 


\ : 
CE 
7 


n= 47 
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10 
20 
30 
40 
50 
60 
70 
80 
110 
120 
1350 
140 
10 
150 
160 
190 
200 
210 
0,4 
,63 


220 


hel 


15% 


REM HOMME DE DEEE DEEE DEEE DEEE 


REM * * 

REM + couches electroniques * 

REM % * 

REM #66 EOMOHMOHEDEDEDEODEDEDEEEENMONEE 

REM 

INFUT “nom de l''element :"3;2#: PRINT #8 , z# 


INPUT "nombre d'electron:'"ss 

READ af,b#,c# 

dt=a#$+b#t+c$ 

u=u+VAL(c#) 

IF uis THEN FRINT d# : FRINT #8 , d# : GOTO 1 


IF u=s THEN PRINT d£ : FRINT #8 , d# : RUN 

IF uïs THEN b=VAL(c#)-(u-s) 

d#f=at+b#+RIGHTS(STRS(B) ,LEN(STR# (D) )—-1) 

PRINT d£ : FRINT #8 , d# 

DATA 1,5,2,2,8s,2,2,p.6,5,s,2,5,p.6,4,s,2,5,d,i 
1Pr6:5:S:24,d,10,5,p,6,6,5,2,4,f,14,5,d,10,6,p 
6,d,10,7,s,2 

RUN 


itim 


cartonne 


is 
£s 


Zp£ 


fluor 
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Zs* 


2pS 


aluminium 


1s= 
2s%? 
Zp4é 
3s2 
3pi 
fer 
1S* 
Zs* 
£pé 
3s2 
3p#4 
4Ss 
544 
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argent 
15% 


at 
HS 


zp& 
3s# 
5p6 
as 
3410 
4pé 
sx 
447 
plomt 
is? 
2s* 
2p4é 
3s2 
‘3p4 
4S2 
3410 
4apé 
3S% 
4410 
Epé 
&SE 
at1i4 
S5d10 


épi 


ur an i LM 
1s£ 
is 
£p4 
35% 
35p4 
4s< 
5410 
4p4 
as 
4410 
Spé 
LS 
4f14 
5410 
£pé 
&dé 
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5 
Programmes généraux 


1. Equation de Van der Waals 


Pour traduire le comportement des gaz, Van der Waals introduit 
en 1873, l'équation (relative à une môle de gaz). 


P+-2)0 -b) = AT 


où la pression est corrigée par le terme a/V? dû aux forces 
d'attraction entre les molécules, et où le volume V est corrigé par le 
volume propre des molécules (covolume) b. 

Les constantes a et b dépendent de la nature du gaz. 

En calculant les dérivées première et seconde de P par rapport 
à V, le long d’une isotherme, et en les annulant, soit : 


RT 2a 





or) S. = + = 0 
aV T  (V-b} VW 

0°P 2RT 6a 
( ) = = 0 








+ _ 
a T (V-b} V* 
On obtient les relations suivantes : 
Vc=3b ; Tc=8a/27Rb et Pc=a/27b°? 
qui permettent, les constantes a et b étant obtenues expérimen- 
talement, le calcul des valeurs critiques Vc, Tc et Pc et vice versa. 


Remarques : 
— La température critique Tc est celle au-dessous de laquelle 
il faut amener un gaz pour qu'il soit liquéfiable par compression. 
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PcVc 3 : ; 
— AT = — alors que pour le gaz parfait on aurait un rapport 
C 
égal à l'unité. 
— Si on cherche les constantes a et b à partir des valeurs criti- 
ques, on trouve : 
2T 2 
64 Pc 8Pc 
— Pour l'hydrogène a=0,244 et b=0,266 1.atm.mole-' avec 
ces unités la constante du gaz parfait R vaut 0,082 1.atm.k-'mole-". 





Parabole de Mariotte 


Reprenons l’équation de Van de Waals et posons Pv = y, il vient : 


2 

ESP SAT 

y y 

C'est l'équation implicite f(y,P) de l’isotherme. Pour que la tan- 

gente soit horizontale, il faut LA = 0 
dP y 

c'est-à-dire : —by? + ay - 2abP = 0. 

C'est une parabole d’axe horizontal passant par l’origine et par 
y = a/b pour P = O. 

Le long de cette parabole la fonction Pv est constante et on peut 
écrire, pour de petites variations de la pression, que 

Pv = Cte 

On l'appelle donc la ‘‘parabole de Mariotte”’. 


Bien que le gaz suive la loi de Mariotte, le produit Pv est dif- 
férent de RT, ce n’est pas un gaz parfait, sauf au point P=O où 
Pv= ab. 


La ‘température de Mariotte”’ est celle de l’isotherme au-dessus 
de laquelle il n’y a plus de minimum (Tm = a/bR). 


Le programme 1 correspond à la représentation de Clapeyron 
(P=f(V) ) pour l'hydrogène. 
Le programme 2 met en évidence la parabole de Mariotte. 
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Liste des variables 


A, B constantes a et b de l’équation de Van der Waals 
pour l'hydrogène 
R constante du gaz parfait 


P,VetT pression, volume et température 
Pc, Vc et Tc valeurs critiques 


Déroulement du programme 1 


lignes 10 à 70 introduction des données 
lignes 80 à 130 tracé des axes 
lignes 140 à 220 tracés des différents isothermes 


Déroulement du programme 2 


10 
20 
30 
40 
50 
60 
70 
80 
90 
100 
110 
120 
150 
140 
150 
160 
170 
180 
190 
200 
210 


220 


lignes 10 à 70 introduction des données 

lignes 70 à 130 tracé des axes 

lignes 130 à 340 tracés des différents isothermes 
lignes 350 à 500 tracé de la parabole de Mariotte 


REM EQUATION DE VAN DER WAALS 
REM HEC 
REM 

REM INTRODUCTION DES CONSTANTES 
REM HOME EONEOHEMEDE DEEE DE DE DE DEN DE DEN DEEE 
A=0.244:5H=0,0266:KR=0,082 
VEC=S34B: TC=8*A/ (27#R4ER) : FC=<A/ (27*+EB*B) 
CLS 

MODE 2 

FLOT O0,0:DRAWK 640,0 

FLOT 0,0: DRAUR 0,400 

LOCATE 2,3:FRINT "F" 

LOCATE 78,24:FRINT "V" 

FOR T=100 TO 20 STEF -10 


FOR V=0.5 TO 0.008 STEF -0.001 

F= (K#T) / (V-H) -A/ (V2) 

IF F£=0 OR 20#F:=400 THEN GOTO 200 
FLOT S500+%V/0.5,20%F 

NEXT V 

IF T=TC THEN STOF 

NEXT T 


TETC :GOTO 150 
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10 

Z20 

30 

40 

50 

60 

70 

80 

90 

100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
300 
#10 
520 
3350 
340 
250 
360 
400 
410 
420 
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REM EQUATION DE VAN DER WAALS 
REM HEURE EEE DEEE DE DE DEEE DE SEX 
REM 
REM ENTREE DES CONSTANTES 
REM #6 De De DE DEEE DEEE DEEE EEE 
A=0.244:B=0.02766:R=0. 082 
VC= ZE: TC=B#A/ (27#RX%EB) : FC=A/ (27*E%H) 
CLS 
MODE 2 

FLOT O,0:DRAWR 640,0 

FLOT 0,0: DRAWR 0,400 

LOCATE 2,23:PRINT "FV" 

LOCATE 78,24:FPRINT "PF" 

FOR T=100 TO 20 STEF -10 

AF=0: AV=O 

FOR V=0.001 TO 7 STEF 0.02 

= (RAT) / (V-B)-A/ (V2) 

IF F£O OR AF“O THEN GOTO 200 

IF F#3:=600 OR AF#3>=600 THEN GOTO 300 
IF F#V#40>=400 OR  AF#AV#x40:=400 THEN GOTO 300 
FLOT F#3,PxVx40: DRAWR APFx3-F#T, AP#AV#40-FxV#40 
AV=V: AF=F 

NEXT V 

IF T=TC THEN GOTO 400 

NEXT T 

IF T=TC THEN GOTO 160 

REM FARABOLE DE MARIOTTE 

REM #46 M DEEE DE DE DE DE DEEE EE 

FOR F=0 TO 50 

KH=A"2-8+% (E"2)+0%F 

IF KH£O THEN STOF 


30 X1=(A+SOR (EH) ) / (2#B) : XZ= (A-SOR (EH) ) / (2*H) 
440 IF X1#40%=400 OK X2#40:=400 THEN GOTO 500 
450 FLOT F#z,Xix40:FLOT F#5,X2%40 
900 NEXT F 





2. Méthode des moindres carrés 
Applications aux équations non linéaires 


En physique, les graphiques expérimentaux présentent en géné- 
ral des irrégularités. Si les points représentés se répartissent suivant 
une configuration assez régulière, le physicien ajustera la ou les gran- 
deurs étudiés ; c’est-à-dire qu'il substituera aux grandeurs observés 
des grandeurs calculées. 


La première étape de l'ajustement porte sur le choix de la courbe 
d'ajustement, la seconde étape porte sur la détermination de l’équa- 
tion de la courbe d’ajustement. Nous nous intéresserons aux équa- 
tions de la forme : Y =f(X). 


Méthode des moindres carrés 


C'est une méthode analytique qui permet de calculer les différents 
paramètres Wi d’une équation de la forme Y =f(X). 
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Représentons graphiquement les points expérimentaux de 
coordonnées : 

(X1, Y1) , (X2, Y2) (XJ, YJ) 

et les points de coordonnées qui se trouvent sur la courbe d’ajus- 
tement, soit : 

(X1, Y’1) , (X2, Y’2) si (XJ, Y'J) 





La fonction Y =f(X) dépendant de plusieurs paramètres Wi des 
valeurs différentes de ces paramètres conduisent à des fonctions 
Y =f(X) différentes. 

La méthode des moindres carrés consiste à retenir la fonction 
Y =f(X) qui rend la plus faible possible la somme des carrés des écarts : 
Yi Y'i. 

Il est hors de question pour nous de donner, dans cet ouvrage, 
la démonstration mathématique rigoureuse du calcul des différents 
paramètres Wi de la fonction Y = f(X). 


Nous utiliserons la méthode de Deming (1), reprise par Wentworth 
(2), le lecteur intéressé pourra éventuellement consulter les articles 
des auteurs. 

Supposons que la fonction d'ajustement contienne la variable X 
et un certain nombre de paramètres : W1, W2..WJ : soit Y’ (X, 
W1...WJ). 


Il s’agit de trouver à l’aide d’un calcul itératif les valeurs optimi- 
sées des différents paramètres Wi (somme des carrés des écarts 
Yi-Y'i minimale). 
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Soit AWi la différence observée entre deux itérations successi- 
ves pour le paramètre Wi. 
Ecrivons que : AWi = Wi — Wiet Y’i = Y (Xi, Wi, W2... Wi). 


Posons F (Wi) = ue Xi, Wi, W3..Wj 


Deming (1) a démontré que si les Xi sont connus avec précision, 
les AWi se calculent en résolvant le système d’équations suivant. 


En posant Fi = Yi - v'i 


J J J 
F(W1).F(W1).AW1 + LF(W1).F(W2).AW2 ie LF(Wi ).F(WJ).AWJ = F(W1).FF 


Eve svt 


J J J 
F(W2).F(W1).AW1 + LF(W2).F(W2).aW2 UE LF(W2).F(WJ).AWJ = LF(W2).Fi 


J J J J 
LF(WJ).F(W1).AW1 + LF(WJ).F(W2).AW2 PSE LF(WJ).F(WJ).AWJ = F(WU).FF 


Les AWi sont calculés à partir d’un sous-programme d’inversion 


de matrice. 

On arrête le calcul lorsque le coefficient de corrélation est con- 
vergent, c'est-à-dire que la somme des carrés des écarts est minimum 
pour la fonction d'ajustement choisie. 


Liste des variables 


N nombre de points expérimentaux 

J nombre de coefficients de l’équation d'ajustement 
X(N) tableau pour les différents Xi 

Y(N) di  R “ Yi 


V(K) tableau pour les dérivés F (Wi) 
Remarque : V(1) = Y(I) — Y’ 
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Z(J,J) tableau pour les différents éléments de la matrice 
principale 
COR coefficient de corrélation 


Déroulement du programme 


lignes 1 à 20 entrée des données (N et J) 
lignes 40 à 130 entrée de chaque couple de points expéri- 

mentaux (Xi, Yi) 
Remarque : à la ligne 70, on rentre les diffé- 
rentes valeurs des coefficients Wi de départ 
(choix de l'utilisateur) 

lignes 140 et 150 on repère les valeurs maximales en valeur 
absolue des Xi et des Yi 

lignes 200 à 250  initialisation des différents éléments de la 


matrice 
ligne 340 on pose V (1) = Y (1) — Y’ 
ligne 370 on pose DEF FNZ (W) = Y (1) - Y’ 


lignes 390 à 410 les différentes dérivées F (Wi) = V(L+1) sont 
calculées numériquement 


ligne 412 on calcule le coefficient de corrélation COR 
après chaque itération 

ligne 500 saut au sous-programme d'inversion de 
matrice 


Les calculs se déroulent sous vos yeux jusqu’à la convergence 
du coefficient de corrélation (Noter les différents Wi après avoir arrêté 
le programme en appuyant 2 fois sur ESC à la convergence). 


lignes 1000 à 1200 l'ordinateur efface l'écran. Il demande à nou- 
veau les différents Wi calculés à la 
convergence. 

ligne 1025 on trace les axes, les points expérimentaux 
à l’échelle et la courbe d'ajustement. 
Remarque : en ligne 1080, n'oubliez pas de 
rentrer la fonction d'ajustement Y’. 


Nous donnons 5 exemples d'utilisation de ce programme. 

Exemple 1 : la fonction d'ajustement est une droite. 

Exemple 2 : la fonction d'ajustement est une fonction du troisième 
degré. 

Pour les 3 derniers exemples, les fonctions d'ajustement sont des 
fonctions d’exponentielle (exemples les plus difficiles à traiter). 
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En effet, Leydet (3) a montré que les fonctions de mélange en ther- 
modynamique pouvaient être représentées par des fonctions d’expo- 
nentielle. Nous donnons donc 3 exemples de mélange. 


Remarque : ces fonctions s’annulent pour x = 0 et x = 1. 
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1 CLS 

9 REM MOINDRES CARRES 

6 REM HOME EU ENENEN 

7 KREM 

10 INFUT "NE DE FOINTS EXF.N="3;N 

20 INPUT "NB DE COEFFTS.J="3;1J 

35 DIM X(N):DIM Y(N):DIM Z(J,J):DIM D(J) 
40 M=0:F=0: DIM W(J+1):DIM K(J) 

42 K=J+1:DIM V(H) 

44 DIM U(J,J):DIM E(J):DIM F(J) 

90 DW=0.001: FR=0Q. 00000001: MI=1E-50 

993 MODE 2 

60 FOR I=i TO J 

70 INPUT "W="3WCI):FRINT "W"s 15; "="; WCT) 
80 NEXT I 

90 CLS 

100 FOR I=i TO N 

110 PRINT “"I="sI 

120 INFUT "X="5;XCI): INFUT “"Y="sY(I) 

150 FRINT "X"sls"ae"s XCI)s "ets "Ya ls "=" s Y CI) SPRINT 
140 IF ABS(Y(I)):F THEN F=ABS(Y(I)) 

150 IF ABS(X(I)):M THEN M=ABS(X(I)) 

160 NEXT I 

199 S=0 

200 FOR C=i TO J 

210 FOR L=i TO J 

220 Z(L,C)=0 

230 NEXT L 

240 NEXT C 
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250 
300 
310 
516 
320 
340 
1—X) 
341 
350 
360 
370 
2) #X 
390 
400 
410 
aii 
412 
415 
420 
440 
445 
450 
ASS 
460 
470 
490 
491 
492 
493 
500 
800 
610 
620 
60 
640 
641 
642 
645 
660 


1000 
1010 
1020 
1025 
1050 
1040 
1041 
1042 
1043 
1050 
1060 
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FOR L=1 TO J:D(L)=O:NEXT L 
REM CALCUL ITERATIF 
REM He EME EME DEEE NE 
FOR 1=i TO N 
X=X (I) 
V(1L)I=Y(I)-(W(1)#(1-X)-WC1)#(1-X)HEXF(-W(2) #X/ ( 
)+W(3)#X-WC3)#XKEXF(-W(4)#C1—-X)/X)) 
D=V(1)#V(1):S=S+D 
FOR L=1 TO J 
W=W(L) : A=W(L) 
DEF FNZ{(W)=Y(I)-(W(1)#(C1-X)-W(1)#0C1-X)#EXF (—-W< 
/ (1-X))+WC3)#X-WCZT)AXAEXP(-W(4)#(C1—-X) /X)) 
W(L)=A+DW: Zi=FN Z(W) 
W(L)=A-DW:Z2=FN Z (W) 
V(L+1)=(21-22) / (2#DW) 
B=V(L+1)#V(1) 
D(L)=D(L)+B 
W(L)=A 
NEXT L 
FOR C=1 TO J 
FOR L=1i TO J 
H=V (L+1)#V(C+1) 
Z(L,C)=2(L,C)+H 
NEXT L 
NEXT C 
NEXT I 
CR=S/N 
FRINT "COR="3;SQOR(CK) 
FRINT 
GOSUE 2021 
FOR Ls=i TO J:R(L)=0 
FOR C=1 TO J 
R(L)=R(L)+Z(L,C)+#D(C) 
NEXT C 
NEXT L 
FOK F=i TO J 
W(FI)=W(F)-R(F):FRINT W(F) 
NEXT F 
GOTO 190 
CES 
FLOT 06,200:DRAWR 626,0 
FLOT O,0:DRKAWR 0,400 
FOR F=1i TO J:INFUT W(F):NEXT F 
FOR I=1 TON 
FLOT X(1)/M#600,Y(1)/F#180+200:DRAWR 2 
FLOT X(1)/M#600,Y(1)/F#180+200:DRAWR © 
FLOT X(1)/M#600,Y (I) /F#180+200:DRAUR — 
FLOT X(1)/M#600,Y(1)/F#180+2700:DRAWR O©,-2 
NEXT I 
FOR KE=i TO 600 


1070 
1080 


ZT) KX— 


1085 
1090 
1100 
1200 
2000 
2010 
2021 
2022 
2025 
2029 
2026 
2070 
2051 
2055 
2056 
2037 
2040 
2045 
2046 
2050 
2060 
2070 
2080 
2100 
2110 
2120 
2150 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
Les Lun Lan] 


2220 





à A 
22:50 


2240 
2250 
2260 
2270 
2290 
2:00 
2400 
2410 
2420 


2430 


X=(M/600)%F 
Y=W(1)#C1-X)-W(1)#(1-X)XEXF(-W(2)#X/ (1—-X))+W 
WCI)KX#EXF(-W(4)#C1-—-X) /X) 

IF ABS(Y/F#180) 200 THEN GOTO 1100 

FLOT E,Y/F#1B80+200 

NEXT K 

STOF 

REM SF. INVERSION DE MATRICE 


REM HEC ME DE DEEE NE 
FOR C=i TO J 

FOR Lzi TO J 

U(L,C)=Z(L,C) 

NEXT L 

NEXT C 

DET=1 

FOK 6=1i TO J 

MA=MI 

FOR L=G TO J 

FOR C=G TO J 

A=Z(L,C) 

IF ABS(A):=MA THEN GOTO 2050 
IM=L : JM=C: MA=AES (A) 

NEXT C 

NEXT L 

IF MAS:MI THEN GOTO 2100 
STOF 

FT=2Z(IM,9JM) : DET=DET#FT 
F(G)=IM:E(G)=JM 

IF IM=G THEN GOTO 160 

FOR C=i TO J 
A=Z(IM,C):2Z(I1M,0)=2Z(6,C):2Z(6G,C)=A 
NEXT € 

ÎF JM=G THEN GOTO 2200 

FOR L=i TO J 
A=Z(L,G):Z(L,G)=2Z(L,JM):Z(L,JM)=A 
NEXT L 

FOK Lei TO J 

IF LzG THEN GOTQO 2500 

ReZ (L,G)/FT 

FOR C=i TO J 


IF C=G THEN GOTO 2290 

DD=Z (L,C)-R#Z(G,C) 

IF ABS(DD)<#ARS(Z(L,C))xFR THEN DD=0 
Z(L,C) =DD 

NEXT C 

NEXT L 


FOR Z=1 TO J 
Z(Z,G6)=2Z{(Z,G)/FT 
Z(6,2)=-27(6,2)/FT 
NEXT 2Z 
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2500 
2510 
2520 
2530 
25340 
25350 
2560 
2570 
2580 
2600 
2610 
2620 
26:50 
28640 
3000 


3010 


Z(6,G)=1/FT 

NEXT G 

FOR G=1 TO J 

Z=J-G+1 

FL=F (2) 

IF FL=Z THEN GOTO 24600 

FOR L=1i TO J 
A=Z(L,2):Z(L,2Z)=2Z(L,FL)s:Z(L,FL)=A 
NEXT L 

FC=E (2) 

IF FC=Z THEN GOTO 3000 

FOR C=i TO J 
A=Z(Z,C):2Z(2Z,C)=Z(FC,C):2(FC,C)=A 
NEXT C 

NEXT G 

RETURN 


EXEMPLE 1 


cOk= 6.m48528138 


di 0.040529138 


de 


CÔR= 0.048528138 


L'ANPE 
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On part de 1 = 1 ; W2= 1 








Y1 = W1.x + W2 


Fig. 1 


dE EXEMPLE 2 


Ya 


F L. 


in 22% 
Va0008 


On part de W1 = 1 ; W2 = 1 ; W3 = 1 et Wk = 
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Y' = W1.x? + H2.x° + W3.Xx + Wk 


EXEMPLE 3 
Tableau : 


MELANGE : 1 : Tri-n-dodecylamine, C36NH75 
2 : Chlorobenzene, C6C1H5 


x(C36NH75) Hm(J.mol- ‘) 
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COR= 18.2947852 


Er 


ET 
COR= 18.2771143 


EL 


SE 
LICE LR 
Es 





Y'=W1.(41-x)-W1.(1-x).Exp (—-W2.x/(1 - x) ) 
On part de W1 = 1000 
W2 = 4 
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Tableau : 
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EXEMPLE 4 


MELANGE : 1 : Benzene, C6H6 
2 : Cyclohexane, C6H12 


x(C6H6) Hm(J.mol-!) 
0.02055 62.71 
0.04847 143.94 
0.10063 283.30 
0.14769 394.02 
0.18974 471.94 
0.22787 552.61 
0.26264 608.74 
0.29442 654.19 
0.32360 690.16 
0.35041 718.98 
0.37522 741.81 
0.41956 773.21 
0.45795 790.71 
0.49163 798.72 
0.51860 800.12 
0.57886 786.04 
0.60473 773.39 
0.63650 751.99 
0.67183 720.49 
0.71129 675.57 
0.75592 611.98 
0.78027 571.43 
0.80625 523.32 
0.83401 466.66 
0.86376 399.34 
0.89570 319.37 
0.93046 223.37 
0.96744 109.53 
0.98720 44.10 


765.297645 


che 4 libère 


LME 


3.38265383 
cÔk= 6.225633226 


FF. 


En 


FES 





Fig. 4 


Y'=W1.(1-x)-W1(1-x).Exp(-W2.x/(1 — x) ) 
+ W3.x — W3.x. Exp (— W4.(1 — x)/x) 


On part de W1 = 1000 
W2 = 4 
W3 = 1000 
Wa = 4 


177 


Tableau : 
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EXEMPLE 5 


MELANGE : 1 : Propanol 


0.0059 
0.0177 
0.0291 
0.0403 
0.051 
0.071 
0.105 
0.149 
0.189 
0.203 
0.229 
0.280 
0.337 
0.384 
0.478 
0.592 
0.714 
0.833 
0.909 


2 : Eau 


Hm 


— 12.9 
— 36.9 
— 56.5 
— 70.5 
— 78.4 
— 81.3 
— 73.4 
— 59.5 
— 48.8 
— 46.6 
— 39.8 
— 22.6 
— 4,1 

+ 9.6 

+ 32.6 
+ 49.7 
+ 51.9 
+ 35.7 
+ 18.9 


222.993277 


CIEL TT 7 
ER 


ALU 
co 2.69838542 


He 


Ji 
dil dr 63530492 


Hu. 


LT CN 


nl 





Fig. 5 


=W1.(-x)-W1(1-x).Exp (— W2.x/(1 — x) ) 
+W3.x — W3.x.Exp (-— W4.(1 — x)/x) 


On part de W1 = —100 
W2 = 30 
W3 = 300 
Wa = 3 
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Vous n'avez pas encore eu le coup de foudre pour la physique ? 


Votre ordinateur et ses programmes, en vous libérant des longs 
calculs, des formules «bombardées», de l'impatience du professeur, 
vous offrent une nouvelle chance. 


De plus, ces programmes sont soigneusement commentés afin 
que cette nouvelle approche de la physique vous paraisse transparente 
et que vous vous sentiez «maître». Ils vous expliqueront également 
par l'exemple les possibilités graphiques de votre ordinateur. 


La mécanique, l'électricité, la thermodynamique, l'optique et 
la chimie sont abordées dans cet ouvrage au niveau des classes de lycée. 


Les auteurs mettent ici à votre service leur expérience de pro- 
fesseurs de physique. Ils ont déjà écrit de nombreux ouvrages d'élec- 
tronique et d'informatique. 


Les programmes qui suivent correspondent à l'enseignement des 
Sciences Physiques en classes de lycée. Îls ont été écrits pour l'ordinateur 
AMSTRAD. Les copies d'écran pour l'illustration de ce livre ont été 
obtenues au moyen de l'imprimante AMSTRAD DMP 2000. 


Nous avons quelquefois sacrifié quelque peu la présentation 
sur l'écran pour que le programme soit moins long à recopier, nous 
n'avons pas recherché systématiquement les astuces de programmation, 
afin que la construction du programme vous apparaisse plus clairement. 
Nous espérons que vous apprécierez ces compromis. 


couverture JÉRÔME LO MONACO 
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